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Thermal Voltage Converters and Comparator for Very 
Accurate AC Voltage Measurements 


E. S. Williams 


Institute for Basic Standards, National Bureau of Standards, Washington, D.C. 20234 


(July 22, 1971) 


A new fourteen-range set of thermal voltage converters and a thermoelement comparator are used 
to measure ac-de difference, and a-c voltages relative to external d-c standards, with 20 ppm (parts- 
per-million) accuracy at audio frequencies. The imprecision is less than 2 ppm. Corrections relative to 
the very stable middle ranges can be redetermined for every range by a seven-step intercomparison of 


certain adjacent ranges. 


Key words: AC-DC difference; comparator; thermoelement; transfer voltmeter; voltage measurements. 


1. Introduction 


An a-c voltage at audio frequencies is measured 
most accurately at the present time by using a thermal 
voltage converter (TVC) to compare it with a stable 
and accurately measured d-c voltage, which is nom- 
inally equal to it. The basic d-c standards are then, 
in effect, “transferred” to the a-c measurement. The 
TYC may be simply a thermoelement (TE) in series 
with an appropriate multiplier resistor. The output 
emf of the TE is ordinarily monitored with a nul 
detector and a balancing circuit, which may be a 
Lindeck potentiometer. A balancing circuit and a null 
detector are included in most commercial multi- 
range models. 

Before the TVC is used for a-c voltage measurements 
it must be tested for ac-de difference or frequency in- 
fluence so that corrections may be applied. The set 
of TVCs and the new TE comparator described in 
this paper were developed primarily for making these 
tests. However they may also be used for a-c voltage 
measurements as explained in section 9. For ac-de 
difference measurements the TE comparator, whose 
read out instrument may be either a nanovoltmeter or a 
galvanometer, provides a considerable advantage over 
other methods in overcoming the difficulties caused 
by power supply instability and inexact voltage control. 

The 14-range set of TVCs (designated No. 7) consists 
of six resistor units and two TEs, and extends from 1 to 
1000 V. The TEs are rated at 2.5 and 5.0 mA, and each 
one may be attached to any one of the resistors by a 
coaxial connector (see fig. 1). This permits each resistor 
to be used for two voltage ranges, as shown in table 
1. Certain adjacent ranges of the set can be inter- 
compared, and the ac-de differences of all the ranges 


can be determined relative to any one range. This 
set is otherwise similar to an earlier set [1]'!, which 
consisted of two 5-mA TEs and 12 series resistors, 
and covered a range from 0.5 to 500 V. The earlier set, 
designated No. 1, has been extended to 1000 V with 
one additional resistor as part of the present project. 
It is evaluated in a 14-step intercomparison in which 
each range is compared, at reduced voltage, with the 
next lower one. 

In a well designed TVC the ac-de difference is not 
affected by changes in the applied voltage. Therefore, 
an ac-de difference determined for the 100-V range at 
60 V can be applied as a correction when this resistor 
is used at 200 V to test the nominal 300-V range of the 
next higher resistor. Experience has shown this step- 
up procedure to be feasible. However if the higher volt- 
age ranges are not carefully built the ac-de differences 
can change due to self-heating effects to be described 
later. 


2. Design and Construction 


The two TEs are mounted in 2-in brass tubes 
(see figs. 1. and 2, upper right) with coaxial connectors 
for attaching to a resistor. The 2.5 mA-TE has a 400 2 
heater and is used alone as the 1-V range. The 5 mA-TE 
has a 125 C2 heater and a 275 2 resistor is added inside 
its enclosure to make a 400 2, 2-V range. 

The middle ranges (six ranges, 3 to 60 V) make use 
of three resistors—0.8 kf, 3.6 k§2, and 12 k{) (see 
table 1). Each of these is a 2W metal-film resistor 
mounted coaxially in a 3-in brass cylinder 4 1/2 in 
long with coaxial connectors at each end. 


‘ Figures in brackets indicate the literature references on page 153. 


145 











FIGURE 1. New set of thermal voltage converters. 


The simple symmetrical geometry of these units 
permits approximate calculations of their reactances 
[2]. The calculations are inexact because of the nec- 
essary assumptions as to end-effects, and they neglect 
any small residual reactances of the resistors them- 
selves. They do indicate however, that the frequency 
error of each resistor unit (without the TE) should be 
less than 1 ppm at 50 kHz, even for tubes smaller than 
2 in in diameter. (The 3-in tube was chosen for the 
higher ranges where more space is necessary and was 
therefore used for the middle ranges also.) 

The higher voltage ranges have larger resistor as- 
semblies to avoid excessive temperature rise. The 
100-200 V unit has a 40 kO spiraled metal-film resistor 
five inches long.? The 300-600 and 500-1000 V units 
have resistors of tin oxide deposited on a 0.81-in glass 
tube 6 in long. A small voltage coefficient was antici- 
pated in these resistors, but measurements have shown 
it to be negligible. 








*This experimental resistor is not ordinarily available commercially. Four 40 kf. metal 
film resistors in series-parallel can also be satisfactory. 


TABLE 1. NBS thermal voltage converter set No. 7 


Fourteen voltage ranges (column 2) are formed with six series 
resistors and two TEs (F 1 and F 2). 














Series Voltage TE Total 

resistor range resistance 
i 
kQ mA kQ 

Sighibeseiesdandl l 2.5 (F 1) 0.4 

vileaaanines 2 5.0 (F 2) 0.4 

3 2.5 1.2 

a | 6 5.0 1.2 

10 2.5 4.0 

eae 20 5.0 4.0 

30 2.5 12.4 

aa | 60 5.0 12.4 

| 100 2.5 40.4 

ee 200 5.0 } 40.4 

300 2.5 120.4 

i sasincksel 600 5.0 120.4 

500 2.5 200.4 

Di sss aseon 1000 5.0 200.4 

eis 3 A. 





Errors in these ranges are caused mainly by capac- 
itance between the resistor assembly and the outer 
casing, which permits alternating current to bypass 
the TE to ground. Therefore more a-c than d-c voltage 
is required for a given TE output. Frequency compensa- 
tion could be provided by placing relatively small 
capacitors in parallel with part of the resistance. 
However, such capacitors probably would not be suf- 
ficiently stable over a long period of time and might 
be affected by temperature changes which occur in the 
resistor enclosure. 

Compensation was therefore provided (as in the 
earlier set No. 1) with an inner shield which is con- 
nected to the input and surrounds the high- or input- 
end of the resistor (see “S” in fig. 2). The shield is 
positioned, relative to the resistor, to control the capac- 
itance currents and provide optimum high-frequency 
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FIGURE 2. Essentials of thermal voltage converters. 
Thermoelement (TE) enclosure, low voltage unit with resistor (R), and high voltage 
unit (100 to 1000 V) with resistor (R) and frequency compensation shield (S). 
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compensation; perfect compensation is possible only 
at one frequency. 


3. 1000-V TVC 


As stated earlier the corrections for the high and low 
voltage ranges are determined relative to the middle 
range by intercomparison tests. It is necessary there- 
fore that the ac-de difference of all ranges be unaf- 
fected by voltage level. The voltage effect is rare in low 
and middle ranges, but it can be troublesome at higher 
ranges where the heat generated by the resistor is 
appreciable. All intercomparison tests were initially 
made at two voltage levels in order to evaluate this 
effect—or to prove its absence. If, for example, the 
nominal 600-V range is compared with the 1000-V 
range at both 400 and 600 V, and the same results are 
obtained, it is unlikely that the 600-V range is affected 
by the voltage increase (the 1000-V range is very 
probably unaffected at these voltages). However both 
ranges are in question until we are assured that the 
1000-V range is also free of this effect. 

The highest range, then presents a special problem. 
It is more likely to change due to increased voltage 
(heating) and there is no higher range to compare it 
with. If it is compared with another 1000-V TVC at 
two or more voltages and no change is observed it is 
possible that both have the same voltage effect. 
However this is unlikely if TVCs are of somewhat 
different design. If intercomparisons of several units 
of difference design show no change between rated 
voltage and a lower voltage we may safely assume 
this effect is negligible in each one. 

Five 1000-V TVCs (table 2) have therefore been 
built, each with differing components, and a network 
of intercomparison tests was made as described in the 
next section. The original set of TVCs [1] was included 
in these tests. 


TABLE 2. Five 1000-V TVCs 











Unit | Resistance| Voltage range with TE 


Fl F2 


Resistor assembly 








z ¥ 


kQ 
A 400 





1000} Four 100 AQ wire- 

wound resistors in 

series, each 1] 3/8 in 
long. 

1000 Four 100 £2 wire- 
wound resistors in 
series, each | 3/8 in 
long. 

1000 Four 200 AQ metal- 
film resistors in 
series-parallel, each 
2 in long. 

1000 One 200 £2 tin oxide 
resistor 6 in long. 

E 400 1000 Four 400 £2 metal- 
film resistors in 
series-parallel, each 
2 in long. 
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The five 1000-V TVCs are designated alphabetically 
(A through E)$ and all except “‘A” have the inner shield, 
as described below, for frequency compensation. The 
resistor in unit “A” is mounted concentrically in a set 
of six brass rings whose potentials are maintained by a 
capacitance divider. The whole assembly is mounted in 
a 4-in brass cylinder 7 in long. 

Self-heating of a TVC may affect either the fre- 
quency compensation or the resistor itself. Apparently 
the dielectric losses in the insulation between con- 
ductors can be affected by heating and change the 
impedance of the resistor. This change occurred within 
a few minutes after a voltage increase or decrease in 
several resistor types, which were therefore discarded. 

The frequency compensating shield may be moved, 
relative to the resistor, by thermal expansion of mount- 
ing parts and cause a change in ac-de difference. This 
occurred in some of our early units and the change 
was slow, taking up to an hour. Later shields were 
therefore mounted very rigidly since a small displace- 
ment will have a large effect on the frequency influence. 

The shield now in use is in two parts. A brass 
cylinder, with one end closed, is mounted firmly 
against the outer end piece (input end) with a poly- 
styrene insulator one-half inch thick (fig. 3). A movable 
tube fits tightly inside this cylinder and the shield 
length adjustment is made by pushing this piece for- 
ward with a small rod inserted through the end piece 
and the insulator. The rod is removed after each 
adjustment. 

The shield is most effective at the leading edge 
where the voltage difference is greatest between the 
resistor and the shield. Therefore the movable part of 
the shield is cut at an angle, so that only part of it 
(see fig. 3) extends outside the fixed cylinder. Thereby 
the effect of small position changes is reduced. Even 
so, the compensation is affected by approximately 0.01 
percent per millimeter of position change in some 
typical units. 

The shield length adjustment is made in small steps, 
and tests are made after each adjustment. When 
optimum compensation is achieved, the unit is opened 
and the shield parts are locked together with three 
screws through their overlapping portions. 


4. TE and TVC Tests 


As stated earlier the middle ranges are the most 
accurate, and they are the base from which higher 
and lower ranges are tested. However they are no 
more accurate than the TE used with them. Before 
being installed in a TVC the TEs are therefore com- 
pared with one of a group of carefully selected TE stand- 
ards which were made according to NBS specifications. 
This group contains twelve elements made by four 
manufacturers, and are rated from 5 to 20 mA. [1] They 
have either Karma or Evanohm‘ heaters to reduce 


3 Unit “D” is assigned to this set (No. 7). 

*Certain commercial products and instruments are identified in this paper in order to 
specify adequately the experimental procedure. In no case does such identification imply 
recommendation or endorsement by the National Bureau of Standards, nor does it imply 
that the products or equipment identified are necessarily the best available for the purpose. 


















































FiGURE 3. Cross section of high voltage TVC. 
Input end is at the left. Fixed and adjustable shield parts are supported by a polystyrene insulator one-half inch 


thick. 


thermoelectric effects in the heaters. They also have 
small reverse d-c differences (less than 200 ppm), 
high bead resistances (over 1000 MQ), and small ac-de 
differences. An additional two are multi-junction 
elements of different design. All of these TEs have 
been intercompared as current converters with agree- 
ment to better than 3 ppm. Theoretical calculations, 
taking into consideration all known sources of error, 
also indicate that the error should not exceed a few 
ppm at audio frequencies. Because of this and because 
of the excellent agreement between these thermo- 
elements of different design and construction, we 
may safely assume that the average ac-de difference 
of the group is 0 to +2 ppm. 

The low-frequency performance of a TVC depends 
on the TE since the reactance of the resistor is entirely 
negligible even at a few kilohertz. The accuracy of 
a TE at low frequency is mainly dependent on the 
length of the heater. A very short heater permits heat 
to flow more readily to the support stems and cooling 
to occur between current peaks. However most com- 
mercial TEs have errors less than 10 ppm at 20 Hz 
and some have errors less than 10 ppm even at 5 Hz. 

In testing selected commercial TEs for the TVCs 
at audio frequencies we often find a small ac-dc 
difference (up to 4 ppm) which is independent of 
frequency but dependent on heater current. It or- 
dinarily decreases at lower heater currents. When 
these TEs are used in a TVC the ac-de difference 
will therefore be slightly voltage dependent. In the 
step-up test process, where the TVC is alternately 
used at 50 and 100 percent of rated current, this rela- 
tively small thermoelement error is introduced at 
every step of the process. If the process is a long one, 
as in the 12-resistor set, there could be a significant 
accumulation of error. TVCs in the new set are inter- 
compared by testing each resistor once with one TE 
at reduced current (60 or 67%), and again with rated 
current on the other TE. If these TEs differ in ac-de 
difference at these currents, there will again be an 
accumulation of error. 

This is one of the advantages of the set with only 
six resistors. If the 20-V range is the starting point 
(it is probably the most stable), we have only four 
steps up to the 1000-V range, and three steps down 
to the 1-V range. The other advantages are that less 
time is required for the intercomparisons and there 
are fewer units to construct. 


Intercomparison data for this set of TVCs are listed 
in table 3 as well as determinations made using TVC 
set No. 1. (Test methods are discussed in later sections 
of the paper). The middle ranges differ by 2 ppm or less 
even though the resistance varies from 1200 to 12000 2. 
Since the resistance enters into the theoreticai valcula- 
tion of ac-de difference, this agreement is an addi- 
tional indication of accuracy. 

Step-down tests show that the two lowest voltage 
ranges usually have a positive correction (i.e., more 
a-c than d-c voltage is required for a given TE output) 
of 5 ppm or less at 50 kHz. There is no resistor in 
series with the TE at the 1 V range, and at the 2-V 
range the TE is a large part of the total impedance. 
These ranges are therefore affected more by the small 
reactance of the TE. The corrections for these ranges 
are more likely to change if a TE is replaced than for 
the higher ranges. A step-down test from the middle 
range is therefore advisable when TE replacements 
are made. 

The five 1000-V TVCs were compared with the 
300- and 500-V ranges of TVC set No. 1] and the 600-V 
range of the new set No. 7. The corrections are listed 
in table 4 for 20 and 50 kHz and all values are averaged. 


TABLE 3. AC-DC difference 








AC-DC difference (ppm) 























Applied 
TVC-a | TVC-b voltage 20 kHz 50 kHz 

| if T T 

|bu—50 | 3. | & |8u—8y 8 | 5 

= 1 ; 
Fl | 1 | |i-3|—3 +1} —3 
F2 | Fl 1 | +1 0 
F2 “a |—2 |-2 +1) —2 
F1-3 | 2 | +1 | —] 
F1-3 3 |-1 |+1 0} -2 
F2-6 6 —1}] 0 0| -3 
F1-10 | F2-6 6 | +1 +] 
F1-10 10 | 0|+1 +1] +1 
F2-20 | 20 | 0| 0 +1| +1 
F1-30 | F2-20 20 | -1 —2 
F1-30 30. | —1|+1 —1| +2 
F2-60 60 | I-1| 0 —1} 0 
F1-100 F2-60 60 —1 | —4 
F1-100 100 |-2| 0 —5| —2 
F2-200 200 | '-2|—] —5| —6 
F1-300 F2-200 200 | —2 | +6 
F1-300 300 | —4|-3 +1] +4 
F2-600 600 | —4|—4 +1} +1 
F1-500 F2-600 500 | —2 +26 
F1-500 500 | —6|—3 27|+30 
i i - 
F2-1000 1000 | |-6| +27) 
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Table 5 lists intercomparison data for the same fre- 
quencies for seven pairings of the five 1000-V TVCs. 
Values for 600 and 1000 V do not differ significantly and 
so they are all averaged together. 

Corrections (5) for each range of the new TVC set 
are listed as determined in comparison with TVC Set 
No. 1. Intercomparison data(5q— 5») are also listed and 
a computed correction (5;) based on the intercompari- 
son, starting with the average of the measured ac-dc 
differences of the 10- and 20-V ranges. The TVCs are 
identified as the combination of a TE (Fl or F2) with 
the appropriate resistor; i.e., F2-600 forms the 600-V 
range. 

Averages of the corrections for the five 1000-V 
TVCs listed in table 4 are shown in the circles in figure 
4, and the connecting lines show the intercomparison 
values from table 5. The agreement between the vari- 
ous units is illustrated by this figure. 

Each of the values listed in tables 4 and 5 is the 
average of at least four separate determinations. 
An analysis of the data taken in a series of 26 inter- 
comparison tests (an average of 5 determinations 
with each one), involving all ranges at two frequen- 
cies and usually two voltages, showed an imprecision 
(3 oa, where Gq is the standard deviation of the aver- 
age of five determinations) of only 1.4 ppm. 


TABLE 4. Ac-de difference corrections for five 1000-volt TVC (A-E) 
The standard instruments are the 300- and 500-V units of Set No. 1 and the 600-V range of 
Set No. 7. 





Standard | Volts | A | B C | p | E 


AC-DC difference at 20 kHz (ppm) 























500-1 500 — 20 —9 —4 —12 

300-1 300 anit aie 
600-7 600 —3 —22 —2 | —12,-—10 

—8 —3 —10 

Ave. | <6 21 —9 —3 —i] 

AC-DC difference at 50 kHz (ppm) 

500-1 500 +21 | -6 +11 +30 =—3 

300-1 300 | | +30 —4 

600-7 600 +30 | -—9 +30 —6 

+23 | +32 ~6 

Ave. | +25 | —7 | +11 | +30 -5 





TABLE 5. Relative ac-dc differences for seven pairings of members of 
the five 1000-V TVCs 
































Volts 5-53 ‘bn-8c | 5--5p | Sp—Se | Se-S, | 5a-Se | 84-8 
AC-DC difference at 20 kHz (ppm) 
600 —11 | + 14| +2 —7| +7| +3 0 
1000 —17 | +13 0 —6 0| +9] +5 
600 —16 \+1 —10 —2 
1000 —17 +5 
Ave -15 |+15| +1] —7| +4] +6] +2 
AC-DC difference at 50 kHz (ppm) 
600 | 99 \+14 +10] —33| +28] -15| +3 
1000 | -22 | +6] +2| —33| +35| —14| +4 
600 | —25 | —37 | —15| +2 
1000 —24 | | +1 
Ave. —23 | +10 | +6| —34| +32] -—15| +2 
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5. Thermal Compensation 


The coaxial connector between the TE and the series 
resistor provides a low-reactance connection for mini- 
mum frequency error. Tests and calculations show the 
error is less than 1 ppm at 50 kHz. However, it also 
permits heat to flow freely from the resistor, which 
dissipates up to five watts on the 1000-V range, to the 
TE which has a temperature coefficient of emf of ap- 
proximately 0.2 percent/° C (at constant input current). 

Heat conduction through the center conductor of 
this connector was reduced, at 300 V and higher, by 
replacing the two polystyrene insulators in the co- 
axial connectors (one in the resistor output and one 
in the TE input) with two made of boron nitride. This 
material is a good electrical insulator and yet it con- 
ducts heat well enough to provide an effective heat 
sink for this conductor. It is used in both ends of the 
high voltage resistors to improve the heat flow from the 
resistor to the relatively heavy brass casing. This also 
helps to reduce the temperature rise in the resistor. 


Heat flow through the outer part of these connectors 
will gradually raise the temperature of the TE casing, 
and this also causes a heating drift. The increase in 
temperature of the TE bulb is slowed by mounting it 
in a short section of brass tubing to increase the 
thermal mass and this tube is thermally insulated from 
the TE casing. Thermal compensation is also added by 
attaching a 10-0 thermistor to the tube containing 
the TE and connecting it in series with the TE output. 
A resistor appropriately chosen to complete the com- 
pensation circuit (see Fig. 7 and appendix I) will 
draw current through the thermistor sufficient to 
compensate for the temperature increase. 

This resistor value may be computed for TEs with 
average temperature coefficients as in appendix I. 
However, temperature coefficients differ considerably 
and, for best results, it is usually necessary to match 
the resistor to the TE and thermistor combination. A 
suitable value can be found by substituting a variable 
resistor outside the TE casing and adjusting it for mini- 
mum warmup drift on a 1000-V range. However, if the 
TVC is to be used with a TE comparator similar to the 
one to be described, it will frequently be connected 





FiGuRE 4. Intercomparisons of 1000-volt TVCs. 

Ac-de difference corrections (ppm) from table 4 for five 1000-V TVCs are shown in 
the circles and the relative differences from table 5 are shown on the interconnecting 
lines. The algebraic difference between any pair should be equal to the figure on the 
connecting line within experimental error. 





+. 


across a 1000-20 divider, and this will affect the com- 
pensation by about 20 percent. A value about 10 per- 
cent higher than the one chosen experimentally should 
provide a satisfactory compromise and reduce the 
drift by a factor of five or more. With the compensation 
and the large thermal time constant of the TE in its 
casing, the drift can be reduced to less than 5 ppm/ 
minute under typical laboratory condition. The error 
from this is almost completely eliminated by taking a 
sequence of readings with a-c, d-c, reversed d-c, and a-c 
voltage applied, as described in the next section. 


6. TVC Test Methods 


A brief review of other comparison methods may be 
useful before describing the TE comparator. Figure 5 
shows two TVCs connected in parallel. The separate 
balancing circuits (or Lindeck Potentiometers) (B) and 
null detectors (D) may or may not be built in. Stable a-c 
and d-c voltage supplies are connected to the two TVCs 
by the ac-de switch. The switching arrangement shown 
is typical and convenient, although others, such as re- 
versing the d-c supply by a separate switch, are satis- 
factory. The d-c supply must have an ungrounded 
output and it should be possible to ground either ter- 
minal as the voltage is reversed. The reversed d-c dif- 
ferences of the TVC can rarely, if ever, be neglected 
in these tests. 

For convenience and clarity, we shall call one TVC 
the “standard” and the other the “test.”” Normally the 
one with known ac-de difference corrections is the 
standard, and the test TVC is usually treated differ- 
ently in the test procedures, as explained below. 

There are two procedures for this intercomparison. 
One might be called a null-balancing method. With a-c 
voltage applied to both TVCs, the balancing circuits of 
the test and standard instruments are adjusted to null 
the detectors. The d-c supply is then switched to the 
TVCs, and the voltage is adjusted to obtain a null bal- 
ance first on the detector of the test instrument and 
then on the standard. Each d-c voltage is carefully 
measured, and the small difference between the two 
should equal the difference between the a-c voltage and 
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FIGURE 5. Two-Potentiometer method for TVC comparisons. 
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one direction of d-c voltage. The measurement is then 
repeated with reversed d-c voltage, and an average is 
taken. 

This method is tedious and time consuming, an ac- 
curate d-c voltage standard is required, and if there is 
any drift in either instrument output due to warmup or 
other reasons, as there very often is, the data may lack 
both accuracy and precision. 

A more satisfactory procedure and one that will mini- 
mize the effect of drift may be called a deflection 
method. Both balancing circuits are adjusted to near 
null, and then a-c and d-c voltages are applied in the 
succession indicated in figure 5 (AC, DC+, DC—, AC). 
For each input the power supply is adjusted to set the 
test instrument detector to null. The deflection of the 
standard TVC detector is read for each setting, and the 
difference between the average reading on a-c voltage 
and d-c voltage is an indication of relative ac-dc differ- 
ence. If the readings are made at nearly equal time 
intervals, and if the drift rate is reasonably constant 
(even though fairly large), the determinations can be 
made with satisfactory precision. 

If the detector of the standard TVC is a galvanometer 
the reading may be in divisions or millimeters, but they 
can be converted to volts or percent ac-dc differences 
by calibrating the detector scale in one of several ways. 
A calibration factor can be determined simply and di- 
rectly by making a small measured change in the d-c 
input and observing the detector response. 

If the detector is a microvoltmeter and if the balanc- 
ing circuit is a Lindeck potentiometer indicating the 
millivolt output of the thermoelement, the calibration 
step can be eliminated. With the four deflection read- 
ings on the detector the difference between the average 
reading on a-c voltage (E,) and those on d-c voltage 
(Eq) is computed from the differences of the correspond- 
ing readings of the microvoltmeter. The ac-dc differ- 
ence of the test instrument is 


Esa roa Esa 


5:=6,+ nE sa 


where the subscript s indicates that these values come 
from the standard instrument, and 6, is the ac-dc differ- 
ence correction for the standard.® 

The factor n relates small changes in thermoelement 
heater current (A/) (at a fixed frequency) to correspond- 
ing changes in output emt (AE/E = nAl/I approximate- 
ly). For a TVC we may substitute nAV/V for the right 
hand expression, so that n=VAE/EAV. The value of 
n is 2 if the thermoelement has a square law response, 
possibly at very low heater currents, but it is usually 
1.7 to 1.9 at rated current. Determinations of n should 
therefore be made in a d-c test at five or more current 
levels by measuring the values indicated. A plot of n 
against E can be made so that values corresponding to 
any emf can be found for substitution in the equation 
above (see sec. 8). 


5 See eq 3 of appendix Il. This general formula is applicable here also. 
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7. TE Comparator 


There are several variations of basic comparator 
circuits which are used to minimize the difficulty 
caused by power supply instability [1, 3]. They employ 
a voltage divider circuit to which two emfs are con- 
nected. The divider is adjusted to null a detector, 
and at this point the divider setting corresponds to 
the ratio of the emfs. If one emf (E;) is held constant 
as a-c and then d-c voltage is applied to the TVCs, 
the other emf (£,) will change if there is a relative 
ac-de difference in the TVCs at the frequency (5; # 5s). 
The resulting unbalance in the divider will produce a 
change in the detector deflection proportional to the 
ac-de difference. However, small fluctuations in the 
power supply will produce nearly equal proportional 
changes in the emfs and therefore the detector will 
not be affected appreciably. The stabilizing effect 
depends on how well matched the time constants 
and response characteristics of the TEs are, but the 
effect usually affords a significant advantage. Also, 
the monitored emf need not be held constant so 
exactly as with a simple balancing circuit. 

In the new comparator (fig. 6) the higher of the 
two emfs is connected across a 1000-0 Kelvin-Varley 
divider, and the lower one is connected to the variable 
tap. The detector is brought to null by adjusting the 
divider. A Lindeck potentiometer and the same de- 
tector (with key k,) is used to monitor the test TVC 
output (E,). The detector is labeled “N” in the dia- 
gram because the use of a nanovoltmeter is suggested. 

The test procedure is similar to that already men- 
tioned. Preliminary settings of the potentiometer 
and divider are first made so the readings of the 
nanovoltmeter will fall near mid-scale. Then with a-c 
voltage applied, the key k; (labeled “SET” on the 
instrument, figs. 8 and 9) is depressed and the appro- 
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FIGURE 6. New TE comparator. 
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priate power supply is adjusted for a null on the 
detector. Key kz, (READ) is then depressed and the 
detector indication is read. The two-step procedure 
is repeated for two directions of d-c voltage and for 
a second a-c input. The average indication in nano- 
volts ® (N,) for a-c input and the average for d-c input 
(Na) is computed and the ac-de difference of the 
test instrument is 


Na—Na 


5: = 5s + Ns E. 


(appendix II eq 5) 


Figure 6 shows the test TVC output (E,) connected 
to the high emf terminal (Ey in fig. 8). If E; is smaller 
than E£, it will be necessary to interchange them, but 
this is a minor inconvenience. If E; is connected to 
the E, terminal the potentiometer is switched to that 
input with switch S» (fig. 8) and E; is monitored as 
before. The same switch also reverses the polarity 
of the detector, and therefore the sign of the indica- 
tion. The slightly modified equation for 5; is eq 11 in 
appendix II. 

A galvanometer may also be used as a detector but 
since the circuit resistance is relatively high a photo- 
electric amplifier is necessary for sufficient sensitivity. 
The scale must be linear over the portion to be used 
and a scale calibration is also necessary to relate the 
readings in divisions or millimeters to a voltage change. 
The calibration can be made quite easily however by 
moving switch S,; from READ to CAL with ke closed 
and observing the resulting deflection change. This 
increases the resistance of the 1000-1 divider by 0.5 2 
and changes the divider current by 0.05 percent. The 
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FIGURE 8. New TE comparator. 


® 8, — 5, is conveniently in parts-per-million if E, is in millivolts and Na — Na is in nanovolts. 








FIGURE 9. Photograph of TE comparator. 


voltage applied to the galvanometer circuit is changed 
by the same amount and the deflection change D, 
is used in eq 10 and 12 in appendix 2 to determine 
5,. A measurement of E, is made only to find the cor- 
responding value of n, from a plot as before. 


8. n Tests 


Measurements of n are conveniently made with the 
comparator. The thermoelement to be tested is con- 
nected to the E,, input and the E, input is unused. A 
key labeled “n TEST” (fig. 6) is locked down and 
switch S. is moved to Ey, connecting the potentiometer 
voltage across the divider. The divider is set arbi- 
trarily but preferably near full range and the potenti- 
ometer is adjusted to null the detector. The input 
voltage is then changed by a small measured amount 
AV/V, and the resulting change in TE output (AF) 
is measured by observing the change in indication of 
the detector (nanovolts or microvolts).? Switch S» 
is then moved to E, and a measurement of TE output 
(E) is made with the potentiometer. The values for 
each factor of the equation given above (n=VAE/ 
EAV) are then known and n is calculated from this 
equation. 

If a galvanometer is used as the detector a measure- 
ment of D,; is made as described above, and the deflec- 
tion change resulting from the measured change in 
input voltage (AV) is designated D.. Then n=0.05 
D./qD,, where q is the percentage change in input 
voltage (100 AV/V). 

Determinations of n need be accurate only to a few 
percent, and, since they are quite stable, retesting is 
not usually necessary. 


9. AC Voltage Measurements 


As suggested earlier, the comparator may be used 
with one TVC to measure an a-c voltage such as the 


7 A small change in the input to a TVC can be measured with a volt box and potenti- 
ometer or by a setting with a calibrated d-c voltage supply. It is also convenient to use a 
resistor with a shorting switch in series with a TVC to introduce a small current change 
(Al) in the thermoelement heater (i.e., 20 1 with a 40-kM, 200-V TVC will give 100 A//] = 
0.05 percent which is equivalent to 100 AV/V). 


input to a digital voltmeter or the output of a calibrated 
a-c voltage supply. Switch S,; has a third position 
marked “POT’R” which disconnects the divider from 
the E, input. If a TVC output is connected to Ey, and 
switch S»2 is moved to that input, the potentiometer and 
detector can be used to measure the TVC output and 
changes in it. The circuit is as figure 5 without the test 
TVC. The d-c voltage is accurately measured and nom- 
inally equal to the a-c voltage. 

It is suggested that d-c voltage be applied to the TVC 
first and a detector reading taken. Then, at nearly equal 
time intervals, a second reading is taken with a-c volt- 
age and a third with reversed d-c voltage. The differ- 
ence between the average of the readings with a d-c 
voltage (Eq) and the one with a-c voltage (Eq) is related 
to the input voltage difference very much as before. 


——— Ee es Ea 
} ac” } acl 1+ nE = 6.) 


With a nanovoltmeter the emfs are read directly, but 
if the detector is a galvanometer a scale sensitivity 
factor must be determined. 


10. Design Details 


The case of the TE comparator has thermal insula- 
tion and electrostatic shielding, and the comparator is 
carefully built to minimize thermal emfs and contact 
resistances. The switch decks of the Kelvin-Varley 
divider and S,; and S» have enclosed silver contacts and 
the keys k; and ke are of a low-thermal type. The back 
contacts on these keys connect a 100-0 shunt across 
the detector to reduce sensitivity and noise until one 
key or the other is depressed. The fine control for the 
divider is a ten-turn helical resistor with a special low- 
thermal sliding contact. The resolution of the uncali- 
brated divider is not quite sufficient for setting an exact 
null on the detector, but setting within a few divisions 
is adequate in the ordinary use of the instrument. 
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11. Conclusion 


ranges has reduced by about half the time required, 
over previous TVC sets, to make a complete intercom- 
parison test. Changes in any member of the set are 
readily detected since frequent intercomparison tests 
are more feasible. 

The TE comparator combines the stabilizing fea- 
ture of other comparators with much of the simplicity 
of the single Lindeck potentiometer or balancing 
circuit. It is relatively inexpensive to build and calcu- 
lations are simplified and direct, especially where the 
detector is a nanovoltmeter. Repeated tests have shown 
that the imprecision is less than 2 ppm and ac-dc trans- 
fer accuracies, with the TVC set, are 10 ppm at audio 
frequencies and 20 ppm up to 50 kHz. 
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13. Appendix | 


The temperature coefficient of emf, B, of the 
average TE is about —0.2 percent/° C, with constant 
input current. Thus a corresponding decrease in the 
resistance, R=R.+R,+R, (see fig. 7) is required 
to maintain a constant voltage across the resistor 
R;,. For a 1 degree increase in temperature, the 


Ry 











FIGURE 7. Thermal compensation for TE. 
R. TE output resistance (about 10). 
R, Thermistor resistance (about 101). 
R, Series resistor to complete compensation circuit. (R- and R, have negligible 
temperature coefhcients.) 


change in thermistor resistance, AR;=ak,/100 
where a is the temperature coefficient of the thermistor 


The use of one resistor unit for two TVC voltage -in percent/° C —typically about — 4 percent/° C. There- 


fore, for thermal compensation, 
B= 100AR,/R=aR,/R 
and R=akR,/B 


With the typical values listed we find R=200 Q, 
and, subtracting R, and R;, the value for R; is 180 2. 


14. Appendix Il 


Development of Equations * 
The ac-de difference of a TVC is defined as 


Ve a V~ 
6=—-__—_4 (1) 
J d 


where lV’, is the a-c voltage and Vq the average of the 
two directions of de-voltage required to produce the 
same output emf. In the circuit of figure 6 the a-c 
and d-c voltages are adjusted to give the same emf 
E,, of the TVC under test, as indicated by a null on 
the detector NV with k; closed. The same voltages are 
applied to the standard TVC. Then, 


Va=Va(l+ &) =VG + 8) (2) 


where Vj, is the d-c voltage required to produce the 
same output emf of the standard as V,, and the sub- 
scripts, ¢ and s, refer to the test and standard instru- 
ments respectively. If Vi;—Va<<1, then, closely 
enough 


Vi—-Va 


aiid Va 


¥rom the definition of n given in the paper 


= Esa - Esa ‘ 
5: — 5s = a (3) 
With the polarities as shown in figure 6, 
N=mE,—E, (4) 


where JN is the detector voltage with kz closed and 
ky open, and m is the divider ratio (with the detector 
resistance Rq >> R;). 

Thus. 


Na—Na 


= § +——_—* 
mean NsE sa 


(5) 


where the subscript a and d have the same meaning 
as before. 


* These equations were developed by F. L. Hermach. 
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If a galvanometer is used as the detector, the current 
through it with ke closed is, by Thevenin’s Theorem, 


mE, ei E, 


l»=———— 


R. (6) 


where R, is the resistance of the galvanometer circuit 
with E, and E,=0. With E; and E, held constant and 
pR, inserted by means of Sy, 


R, 





(7) 


I ge 


where R; is the resistance of the galvanometer circuit. 
Since mE,—E,<<E, and p<<1, we have to a 
sufficient degree of approximation, 





Esp _(Dy—De) _Di 


where D, and D, are the resulting galvanometer deflec- 
tions and S is the galvanometer current sensitivity. 
In the ac-dc test, with mE; constant, 


Esa— Esa __Da— Da 





Thus, from (3), (8) and (9), 
D —Da 
au 2 (10) 
st 1 


The test and standard TVCs may be interchanged 
if E;< Es. If E; is applied to the E,, input of figure 8 and 


held constant, with the detector and potentiometer con- 
nected as shown in the figure, a similar analysis leads 
to the following equations: 





i= 5,4 TeoNs (11) 
and 
8.= 5,42 Ce Pe) (12) 
The characteristic n is determined from 
AE/E 
n=AVIV (13) 


by applying known changes in input voltage and 
observing the changes in output emf with a high resist- 
ance voltmeter, as described in the text. If a galva- 
nometer is used as the detector instead of the voltmeter 
its sensitivity is determined by inserting pR, as de- 
scribed. From equations similar to (8) and (9) 


AE, _pAD 
gE dD, 





(14) 


If the thermocouple resistance of the test TVC, 
Ric, is significant, p in equation (7) should be replaced 
by p’, the fraction of Ri + Ric inserted by S;. Similarly 
if the potentiometer resistance R, is significant, p 
in equation (14) should be replaced by p’’, the fraction 
of R,+ Ry, inserted by S. 


(Paper 75C3&4-321) 
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Formulation and Experimental Verification of an 
Axisymmetric Finite-Element Structural Analysis 
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A finite-element structural analysis is described for application to problems in which geometry 
and loading are axisymmetric and material properties are isotropic elastic. An attempt to minimize 
restrictions imposed on the shape and orientation of the triangular finite-elements has been largely 
successful. This facilitates use of the analysis, with automatic finite-element mesh generation, in 
parameter or optimization studies. A series of laboratory tests to verify the analysis are described 
in which the magnitude and distribution of boundary loading was known within narrow limits. 


Key words: Axisymmetric; elastic; experiment; finite-element analysis; force transducer; stiffness 


matrix; structural analysis. 


1. Introduction 


In the structural optimization of a force transducer 
of the elastic column type, one is concerned with 
relationships between surface strains and such 
parameters as boundary loading, dimensions, and 
material properties. Finite-element analysis is prob- 
ably the most versatile structural analysis method 
currently available for studying these relationships. 
The finite-element formulation described here was 
developed for application to such optimization prob- 
lems in which both the geometry and loading are axi- 
symmetric and the materials are isotropic elastic. 
This class of problems is analytically two-dimensional. 
Ring-shaped finite elements of triangular cross section 
and a global Cartesian coordinate system are used in 
the analysis. This approach generally follows the out- 
lines given by Clough [1, 2], Rashid [2, 3, 4], and 
Wilson [5]. 

For efficient structural optimization it is necessary 
to program the computer to generate automatically 
the mesh that subdivides a structure into a network 
of finite elements. To facilitate this, an attempt has 
been made in this formulation to keep to a minimum 
the restrictions imposed on shape and orientation of 
the triangular elements. Explicit integration, used 
here to develop general formulas for the stiffness of 
finite elements, yields logarithmic terms which re- 
quire special treatment to avoid excessive errors 
when one side of a triangular element is oriented at 


*Present address: Naval Scientific and Technical Intelligence Center, Washington, 
D.C. 20390 
' Figures in brackets indicate literature references at end of paper. 
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a small non-zero angle with respect to the axis of 
symmetry. The special treatment adopted here is to 
expand the logarithmic integrals into rapidly converg- 
ing infinite series. Although the stiffness integrals 
are rather intricate in both the logarithmic and the 
series forms, they do avoid the need for numerical 
integration. 

An experimental verification of the analysis was 
obtained by a series of tests conducted on the struc- 
tural body of a force transducer (load cell) made of 
18 percent nickel 250 grade maraging steel. Axial 
compressive load was applied to the load-cell body. 
The load was distributed over a small circular area 
at the top and over a narrow ring-shaped area at the 
bottom. Thus, the location of applied load was known 
within narrow limits. Strains on the surface of the 
load-cell body, known to vary significantly with load 
location or distribution, were measured and found 
to be in good agreement with the finite-element 
analysis. 


2. Finite-Element Formulation 


In a finite-element analysis the continuum structure 
is subdivided into a network of elements that are 
connected to adjacent elements only at common nodal 
points. Elastic displacements within the individual 
elements are assumed to be defined by generalized 
functions that assure displacement compatibility 
along common boundaries of adjacent elements. 
The stiffness matrix of each element, relating nodal 
point forces and displacements, is then computed in 
terms of the assumed displacement functions, dimen- 
sions, and material properties of the element. The 
stiffness matrix for the entire structure, relating 
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FIGURE 1. Longitudinal cross section of load-cell body used in labo- 


ratory tests 


applied external forces to nodal point displacements 
throughout the structure, is formed by superposing 
the element stiffness matrices. 

A longitudinal cross sec ion of the load-cell body 
used in the laboratory tes 1s shown in figure 1. The 
figure illustrates a triangu.ar finite element mesh. A 
similar but finer mesh was used for the analyses that 
were correlated with the laboratory tests. 

The triangular cross section of a general element is 
shown in figure 2. Each nodal point (actually, a nodal 
circle in the axisymmetric case) has independent 
displacement components in the r and z directions. 
The generalized displacement functions assumed for 
this formulation are given by 
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u wi lr-nz—2390 O 0 Qs (1) 


v 0 0 0 lr-nz—2Z] |& 








in which 

u=displacement in the r direction of any point in 
the element, 

v=displacement in the z direction of any point in 
the element, 

and 
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FIGURE 2. Cross section of a general axisymmetric finite element. 


Q@», = a generalized coordinate. 


A generalized coordinate, a,, represents either the 
r or z component of displacement of nodal point i 
or a rate of change of a displacement component 
with respect to r or z. The assumed displacement 
functions give linear variations in displacement along 
element boundaries and, therefore, complete dis- 
placement compatibility between adjacent elements. 
Stresses within adjacent elements are not, in general, 
in equilibrium along common boundaries, but the 
resultant forces acting at nodal points are required 
to be in equilibrium. Substitution of nodal point 
coordinates (defined in fig. 2) into eq (1) gives the nodal 
point displacements in terms of the generalized 
coordinates. 


Uj 182¢e¢ 0° 0 0 ay 
uj 1 aj bj 0 0 0 Qe 
Ux lab 0 0 0 Qs ; 
v{=|}0 001 0 04 jas’ (2) 
vj 0001 a b Os 
Vk 0 0 0 l ak by as 

or {w}=[A] {a}- (3) 


Element strains are obtained by applying the defini- 
tions of strain from elasticity theory [6] to eq (1) as 
follows: 
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{e} =[ B (r,z)] {a} - 


The assumed stress-strain relationship for isotropic 
elastic materials is 


(5) 


or 
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or {o} =[D] {e}, (7) 


in which E = Modulus of elasticity 
and v= Poisson’s ratio. 

By a process of equating internal to external virtual 
work, using eqs (3), (5), and (7), the general equation 
for nodal point stiffness of an element can be shown 
to be (see ref 1, 2, 3 and 7) 


[k]=[4- 1" fo [B]"[D}(B)dr ) [A-"], (8) 


or [A] = [4-*]"(k} [44]. (9) 
In eq (8), dV is a differential volume within the element. 
Integration is performed over the volume of the 
element, and the superscripts 7 and —1 indicate 
matrix transposition and inversion, respectively. 
The stiffness matrix [k] relates the six components 
of force {f} acting at the three nodal points of a 
triangular element to the six components of nodal 
point displacement {w} by 


{ft=[k] {u}. 
Using r, z—z;, and @ as the variables of integration, 


and defining the integrals of the functions of these 
variables by 


(10) 


I(f(r,z—zi, 0))= 2m | | f(r, z—z) rd(z—z)dr,(11) 


the nonzero elements of [k], eq (9), are 
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(1+v)(1—2p) 
For explicit integration of eq (11) in the r—z plane, 
the notation shown in figure 2 is used with the con- 
straints 
nS rr, 
Tj < Tj, 


if r;>= rz then 2 < Zr. 


Since r; may be less than, equal to, or greater than 
rk, four possible sets of integration limits, cases 1 
through 4 in figure 3, are considered. These are: 


Case 1: 


h(ftr,2—2)) =20 |" | ftr,2— ard e—ziddr (12) 
rj 2ij 


Case 2: 


L(f(r,2—2)) = 20 [? [P feyz—a)rdle— aide (13) 
Tr 2ij 











FIGURE 3. Four different integration limit cases. 


Case 3: 
Is(ftr.2—2i)) =20 [” [ fir.2—z)rd(e— adr (14) 
Case 4: 
L(f(r,z—zi)) = 2m |” jin f(r,z—z)rd(z—z)dr (15) 


in which the limit notation zj represents the value 
of z—z along line i—j. Three of the functions in 
[k] integrate directly, for either orientation of side 


j—k, to 
2 
1)=3 TA (ri t+rj+TrK), 
1(=)=204, 
r 
and 


“25 2 
i )=5 1A (bj) + bx), 


in which A= $ (ajby — axb;) =area of the element 


cross section. 
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The integrals of the other three functions in [k], for the 
four cases shown in figure 3, are 
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The logarithmic terms in the above integrals re- 
quire special treatment to avoid excessive errors in 
numerical computations. By a limiting process employ- 
ing L’Hospital’s rule it can be shown that, as rj or both 
r; and rx; tend to zero, the limit of the logarithmic term 
is zero. Therefore, if one or two nodal points of an 
element lie on the z axis, the logarithmic terms are 
omitted from the computation. A more difficult com- 
putational problem is presented when one side of a 
finite element is oriented at a small nonzero angle 
relative to the z axis. To permit the use of elements with 
a side so oriented, the logarithmic terms can be ex- 
panded in the series (Peirce’s formula 768 [8]) 


ee fe ee —})aiX 
In (1+x)=x 9 +3 ati: tl 1) 4 


which converges in the interval —1<y< 1. This 
gives the following alternate forms of the integrals: 
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Notice that use of the series forms of two of the inte- 
grals, I: [(z—zi)?/r?] and [j[(z—z;)?/r2], is limited to 
elements for which either 5; or b, equals zero; that is, 
one side of the element must be oriented in the r direc- 
tion. If an element has one side oriented at a small non- 
zero angle with respect to the z direction, no nodal point 
on the z axis, and no side oriented in the r direction, 
some other method is required to evaluate these two 
integrals. In such an event, the two integrals could be 
evaluated by the application of numerical integration 
to eqs (13) and (15) for the function (z—z;)?/r°. Com- 
putational experience indicates that the series forms 
of the integrals are efficient for numerical computations 
within the bounds: 


: l ' ; : 
Case 1: => => 10-® and 4 < 107! 
pA aj i 
: ] ik : ik 
Case 2: => > 10-6 and * < 10-1 
2 aj rj 
7 ] : 
Case 3 ~>> 10-6 and “ < 10-1 
v4 ar Tj 
i | jk : A jk 
Case 4: = ee 10-§ and — < 10-'. 
é ar rk 


For cases falling below the lower bound it is satisfac- 
tory to set the integrals equal to zero. For cases falling 
above either upper bound the general integrals con- 


taining the logarithmic terms can be used without 
difficulty. 


The stiffness matrices of the various finite elements. 
[k] in eqs (8) and (10), are superposed to form the 
stiffness matrix of the entire structure. This is done 
by adding the stiffness matrix elements that relate 
displacements to resulting forces for common or ad- 
jacent nodal points of adjacent elements. The resulting 
stiffness matrix of the entire structure [K] relates the 


forces applied to the structure {F} to the resulting 


nodal point displacements {w} 


equation (F} = [K] {u}. 


according to the 
(16) 


For problems involving large numbers of unknown 
displacements, eq (16) can be solved efficiently by 
iteration. For the work reported here, Gauss-Seidel 
iteration was used along with overrelaxation and group 
relaxation as described by Wilson [5]. For the analysis 
mesh shown in figure 1, 250 cycles of iteration resulted 
in convergence to a state having an absolute sum of 
unbalanced residual forces of one-half percent of 
the applied load. This iteration required about 48.8 s. 
running time on a UNIVAC 1108 computer. Evaluation 
of the stiffness integrals, eqs (12) through (15), and 
formation of the stiffness array for the entire structure 
required only about 6.1 s. computer running time. 


3. Experimental Verification 


To determine how well the finite-element analysis 
described the structural response of a load cell, 
a series of laboratory tests were conducted on the load- 
cell body dimensioned in figure 1. The load-cell body 
was made of 18 percent nickel 250 grade maraging 
steel. For the principal verification tests, an axisym- 
metric compressive load of 100,000 lbf (444,800 N) 
was applied through a small circular area on the top 
surface and through a circular ring-shaped area on 
the bottom surface. In a preliminary test, an axisym- 
metric compressive load of 800,000 lbf (3,558,400 N) 
was applied to the load-cell body through two solid 
cylindrical mild steel blocks. In all tests, the loads were 
applied by dead weight testing machines. 

Strains on the surface of the load-cell body were 
monitored by fifty-four metal foil strain gages of 0.062 
in (0.16 cm) gage length. Forty-six gages were oriented 
in the longitudinal direction and located on three 
lines defined by the intersections of the surface and 
three radial planes spaced 120 degrees apart. Eight 
gages were oriented in the circumferential direction 
and located beside longitudinal gages. Twenty-four 
longitudinal gages were located in sets of three at 
common heights and at the three angular locations to 
detect bending due to load eccentricity. 


In the preliminary test, figure 4, the upper mild 
steel block was plastically indented by the spherical 
top surface of the load cell, giving a contact surface 
of about 1.41 in (3.58 cm) radius at 800,000 Ilbf 
(3,558,400 N) applied load. The lower mild steel block 
was used to permit plastic indentation by a small 
region that protruded about 0.002 in (0.005 cm) 
at the center of the bottom surface of the load-cell 
body. This resulted in continuous contact over the 
bottom surface upon application of the first 50,000 
Ibf (222,400 N) load increment. Recorded longitudinal 
surface strains due to a load of 800,000 lbf (3,558,400 N) 
are plotted in figure 5. The curves in figure 5 were 
obtained by finite-element analysis. One curve is for 
800,000 Ibf (3,558,400 N) load applied at the centers 
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FIGURE 4. Setup for applying 800,000 bf (3,558,400 N) compressive 
load ina preliminary test. 


of the top and bottom surface. The other curve is for 
800,000 Ibf (3,558,400 N) load applied at the radii of 
the edges of the contact surfaces in the laboratory test, 
figure 4. Although most of the test data points lie near 
or are bracketed by the analytical curves, this experi- 
ment is clearly insufficient for a precise verification 
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FIGURE 6. Setup for applying 100,000 lbf (444,800 N) compressive 
load in one of the verification tests. 


of the analysis. Figure 5 shows that surface strains 
can be very sensitive to variations in the distribution 
of boundary loading. 


Figure 6 shows the loading arrangement for one of 
the principal verification tests. Axisymmetric compres- 
sive loads of 100,000 Ibf (444,800 N) were applied to the 
top surface of the load-cell body through a mild steel 
block and to the bottom surface through a 0.2 in (0.5 cm) 
wide by 0.035 in (0.088 cm) thick steel ring. Rings with 
mean radii of 0.40, 0.79, 1.58 and 2.27 in (1.02, 2.01, 4.01 
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Ficure 5. Longitudinal surface strains for an applied load of 800,000 
lbf (3,558,400 N). 
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and 5.77 cm) were used in different tests. Figure 7 
shows the recorded longitudinal surface strains for 
four of these tests (the plotted data points) along with 
corresponding curves obtained by finite-element 
analysis. For the finite-element analysis the bottom 
loading radius corresponds to the mean loading ring 
radius and the top loading radius corresponds to one- 
half the contact surface radius. Much of the scatter in 
the data is probably due to bending. The value of the 
modulus of elasticity used in the finite-element 
analysis was 26.8 < 10® lb/in? (1.85 x 10!° N/m2). This 
value was computed by averaging, for several tests, the 
data obtained frbm the three strain gages located at 
z=4.31 in (10.95 cm). This particular data was very 
repeatable and bending errors were essentially can- 
celled by the averaging. Admittedly, this procedure for 
determining the elastic modulus imposes good agree- 
ment between analysis and experiment in the region 
near z= 4.31 in (10.95 cm). A Poisson’s ratio of 9.3, 
a value listed in the literature of steel producers for 18 
percent nickel 250 grade maraging steel, was used in 
the finite-element analysis. The generally good agree- 
ment between the analysis and this test data for the 
entire length of the load-cell body (fig. 7) indicates 
that the analysis would be suitable for the study of 
such subjects as stress concentration and the optimum 
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FIGURE 7. Longitudinal surface strains for an applied load of 100,000 
lbf (444.800 N). 





geometry of material test specimens and fixtures. 

For use as a force transducer, a load-cell body of the 
type tested is ordinarily instrumented with strain 
gages over only a small region near mid-height. An 
important consideration in the design of such a load 
cell is the change in surface strain sensed by these 
strain gages due to a change in the distribution of 
boundary loading. A portion of the curves and test 
data in figure 7 that bear directly on this question are 
replotted to a much greater scale in figure 8. The test 
data points in figure 8 have been averaged from the 
strain readings of sets of three gages at different 
angular locations to minimize bending errors. Figure 9 
is a plot of the differences between surface strain for a 
particular radius of load at the bottom surface and the 
corresponding strain for a 2.27 in (5.77 cm) radius of 
load. Figures 8 and 9 indicate very good agreement 
between the finite-element analysis and the test re- 
sults at the three strain gage locations represented. 
The greatest difference between the analytical and aver- 
aged experimental values represented in figure 9 is 
2x 10-*. This difference approaches the resolution of 
the digital indicator (1 <X 10~*) used in the laboratory 
tests. This data verifies the accuracy of the analysis in 
predicting the sensitivity of such a load cell to changes 
in the distribution of boundary loading. 
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FIGURE 8. Longitudinal surface strains averaged for sets of three 
gages at common heights and different angular locations. 
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FIGURE 9. Difference between surface strain for a particular loading 
radius and the corresponding strain for a loading radius of 


2.27 in (5.77 cm). 


4. Conclusion 


A finite-element structural analysis has been for- 
mulated for application to problems in which geometry 
and loading are axisymmetric and materials are iso- 
tropic elastic. An attempt to minimize restrictions im- 
posed on the shape and orientation of the triangular 
elements employed has been largely successful. This 
facilitates use of the analysis, with automatic finite- 


element mesh generation, in parameter or optimization 
studies. 

An experimental verification of the analysis has been 
obtained by a series of tests conducted on a load- 
cell body. In these tests the magnitude and distribution 
of boundary loading was known within narrow limits. 
Strains measured on the surface of the load-cell body 
were in good agreement with the finite-element 
analysis. 


C. H. Melton was largely responsible for instrument- 
ing the load cell for the measurement of surface strains. 
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A Burnett Apparatus for the Accurate Determination of Gas 
Compressibility Factors and Second Virial Coefficients, 
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Results for Argon and Carbon Dioxide 
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A Burnett apparatus has been developed for the determination of the compressibility factor of 
gases, including corrosive halogenic gases, to an accuracy of 0.1 percent to 0.01 percent at temper- 
atures from 0 to 225 °C and for pressures up to 250 bars. The apparatus has been used throughout this 
temperature and pressure range and at temperatures within 15°C of the critical temperature of the 
sample gas. The apparatus features a rugged yet highly sensitive and reproducible null-type pressure 
transducer to isolate the sample gas without introducing a pressure uncertainty greater than the pre- 
cision of the pressure measurements, which is better than | part in 50,000. The pressure measurements 
are accurate to at least | part in 20,000. The reduction of the data is discussed extensively. The ca- 
pability of the apparatus has been evaluated in terms of the compressibility factor represented by a finite 
pressure or density virial expansion and derived from argon data at 25 °C and pressures up to 250 bars 
and in terms of the second virial coefficient of carbon dioxide based on data at 50°C and pressures 
below 35 bars. The value obtained for this second virial coefhiicient. — 102.2 + 0.2 cm*/mol, compares 
favorably with the value of —102.5 cm*/mol obtained from a recent correlation and is in substantial 
disagreement with the value of —100.7 +0.4 cm*/mol determined by Dadson and co-workers from the 
Burnett method and the piezometer method. 


Key words: Adsorption, gas on metal; argon, compressibility factor of; Burnett data reduction; carbon 
dioxide, second virial coefficient of; compressibility factor; gas, compressibility factor of; gas, PVT 





properties of; measurements, gas; second virial coefficient. 


1. Introduction 


The accurate determination of the equation of state 
of a gas over a wide range of temperature and pressure 
provides important information about the gas both 
from a fundamental and an applied point of view. 
For example, the thermodynamic properties of a gas 
can be reliably extended beyond its ideal gas proper- 
ties if its equation of state is known with sufficient 
accuracy. Equation of state data are conveniently 
represented by means of the compressibility factor. 
The compressibility faetor Z is defined in terms of the 
gas constant R and the measurable gas variables pres- 
sure P, volume V, absolute temperature 7, and quan- 
tity, expressed as the number of moles n, by the rela- 
tionship 

Z=PV/nRT. (1) 
The accuracy of Z obtained from a straightforward 
measurement of each variable is limited by the 
accuracy with which the difficult measurements of 
n and V can be made. The direct or indirect measure- 
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ment of one (n) or both (n and V) of the difficult vari- 
ables can be eliminated by the use of the virial expan- 
sions of Z as a function of P or of the density p and the 
fact that the dependence of Z on P is dominantly 
linear at sufficiently low pressures and approaches the 
ideal gas law 


Z=PV/nRT=1 (2) 


in the limit as P goes to zero [1,2].! The Burnett PVT 
method [3-7] is an experimentally simplified method 
in which neither n nor Vis measured; instead, accurate 
isothermal pressure measurements are made before 
and after stepwise expansions between two volumes. 
The decrease in the overall measurement error 
achieved by the elimination of n and V measurements 
can in principle lead to more accurate values of Z 
than can be obtained from a method requiring the meas- 
urement of each variable in eq (1). For this reason we 
have developed a Burnett apparatus for the accurate 
determination of the compressibility factor of gases, 


! Figures in brackets indicate the literature references at the end of this paper. 





including corrosive halogenic gases [8], at tempera- 
tures from 0 to 225 °C and for pressures up to 250 bars 
(1 bar= 105 N - m~*). The apparatus has been used by 
us throughout this temperature and pressure range 
including temperatures within 15°C of the critical 
temperature of the sample gas. The apparatus is 
equally well-suited and has been used for the accurate 
determination of second virial coefficients. 

This paper is concerned with both the apparatus 
for obtaining Burnett data and the methods for re- 
ducing the data. These subjects are described briefly 
in the remainder of this section. In the next section the 
apparatus is discussed in more detail, and its dis- 
tinctive features are emphasized. Since the reduction 
of the data obtained with the apparatus is an essential 
element in the evaluation of its capability, the data 
reduction is discussed extensively in the following 
section. The discussion is illustrated with results for 
the compressibility factor of argon at 25°C and pres- 
sures up to 250 bars and for the second virial coef- 
ficient of carbon dioxide at 50 °C. 

The apparatus, shown schematically in figure 1, 
consists essentially of a sample volume V, and an ex- 
pansion volume V,, as well as components for (1) 
isolating, evacuating, or filling both volumes, (2) 
measuring the pressure in the sample volume, and (3) 
thermostating and measuring the temperature of the 
apparatus. The gas in the sample volume is isolated 
from the expansion volume by an expansion valve and 
from the external system by a metal diaphragm which 
is the pressure-responsive element of a_null-type 
transducer. A counterbalancing gas acts on the other 
side of the diaphragm; the diaphragm is in its null 
position when the pressure on each side of the dia- 
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phragm is the same. The pressure of the sample gas 
is obtained indirectly from the measurement of the 
pressure of the counterbalancing gas with a piston gage 
and a barometer. The sample and the expansion vol- 
umes include the appropriate volumes of the trans- 
ducer, expansion and filling-evacuation valves, and the 
interconnecting pressure lines; in our apparatus all 
these volumes are thermostated at a temperature which 
is regulated to about 0.001 °C and measured with a 
platinum resistance thermometer. Because no correc- 
tion needs to be made for temperature differences 
among these volumes, a correction that is often a 
significant source of error in PVT measurements is 
thereby eliminated. The transducer produces negli- 
gible variations of the sample volume and has a negli- 
gible effect on the high precision and accuracy of the 
pressure measurements. This is due to the excellent 
reproducibility of the null position of the diaphragm 
and the high sensitivity of the transducer. 

The experimental procedure is begun by filling 
the sample volume to the initial pressure, thermostat- 
ing the apparatus at the isotherm temperature T,, 
isolating the volumes from each other, and then 
measuring the gas pressure Py in the sample volume. 
The corresponding compressibility factor Zo is 


Zo=PWV(Po)/noRT-, (3) 


where no represents the initial number of moles and 
V\(Po) indicates that the sample volume is pressure 
dependent as a result of its distortion under pressure. 

In the next step, which we call the expansion step, 
the expansion volume is evacuated and _ isolated, 
and the gas in the sample volume is allowed to expand 
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into it. After temperature and pressure equilibration 
has occurred, the pressure P,; of the gas that now 
occupies both volumes is measured. For this first 
expansion step the compressibility factor Z, is 


Zi:=P,[V(Pi1) +Vu(P1) |/noRT-, (4) 


where no and 7, are the same as in eq (3). By taking the 
ratio of eqs (3) and (4) and defining N, as the ratio of 
the combined sample and expansion volumes to the 
sample volume, i.e., Ni; = [V,(P:) + Vig(P1)]/Vi (Po), 


one obtains 
PoP: =NiZolZi, (5) 
in which N, depends slightly on both Po and P,;. 


The ratio of the volumes can be written as 
N, = N(1 + AN, (Po, P:)], where N, known as the cell 
constant, is the value of the volume ratio in the limit as 
the pressure goes to zero and AN,(Po, P;) is a cell 
constant correction for the dependence of the cell 
volumes on pressure. 

Generally the sample volume is isolated again and 
the expansion step is repeated for the reduced quan- 
tity of gas n, that is now in the sample volume. This 
leads to a new pressure ratio, corresponding to eq (5), 


P,/P2=N2Z,/Z2, (6) 


where N2=N[1+AN.2(P;, P2)]. This procedure is 
usually repeated to successively lower pressures with 
the successively smaller quantities of gas that remain 
in the sample volume after each expansion until the 
lowest pressures that can be measured both accurately 
and precisely are reached. The generalized form of eqs 
(5) and (6) for the rth expansion is 


P, JP N,Z;~i1Z-. (7) 


The combination of eq (7) and the relationship 
P=pZRT, obtained from eq (1), leads to a correspond- 
ing Burnett equation in terms of the density 


Pr= Pr- i/Ny. (8) 


Equations (7) and (8) are basic equations in the re- 
duction of Burnett data and may be transformed into 
different relationships in which the compressibility 
factor is represented by a finite virial expansion in either 
pressure or density. Any one of the relationships may 
be applied to Burnett data to form an overdetermined 
set of equations for evaluating the virial coefficients 
and generally the cell constant by means of nonlinear 
statistics [9]. Regardless of the relationship used, the 
cell constant is not determined from volume measure- 
ments. In the data reduction the pressure corrections 
AN,, AN2, . . . are not considered as unknown quan- 


tities but are estimated from theoretical and empirical 
distortion equations for pressure vessels or are deter- 
mined from a separate distortion experiment. 
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2. Apparatus and Experimental Procedure 


2.1. Burnett Vessel and Gas Manifold System 


An overall view of the apparatus is shown in figure 2. 
In the Burnett vessel, as shown in figure 3, the sample 
chamber, 0.022 m i.d. by 0.076 m long, and the expan- 
sion chamber, 0.022 m i.d. by 0.051 m long, are bored 
from opposite ends of 0.064 m o.d. nickel rod to within 
0.051 m of each other. Nonrotating end closures are 
used for coupling to the transducer and to the filling- 
evacuation line. The top gland nut is split to accommo- 
date the transducer coupling. The coupling to the 
expansion valve is made through 0.0032 m o.d. by 
0.00080 m i.d. stainless steel capillary tubing silver- 
brazed into the sample chamber inlet. The coupling 
from the expansion chamber is made through similar 
capillary tubing silver-brazed into a T-connection which 
in turn is coupled to the expansion valve and to the 
filling-evacuation valve. A well for the platinum resist- 
ance thermometer is bored radially into the Burnett 
vessel between the sample and the expansion chambers. 

Since halogenic gases react readily with most met- 
als, particularly if a trace of moisture is present, the 
components exposed to the sample gas are fabricated 
from metals such as nickel and stainless steel that are 
resistant to corrosion. We have found that gold provides 
still better resistance to corrosion than nickel and 
stainless steel do, and have therefore gold-plated every 
surface exposed to the sample gas in the Burnett vol- 
umes except the surfaces in the Burnett valves. Ex- 
treme care is taken to keep the system free of moisture 
at all times. The sample gas is not reclaimed; if neces- 
sary, neutralizing columns are used for its disposal. 
Cold traps prevent corrosive gases from entering the 
vacuum system, which includes low and high vacuum 
pumps, a cold cathode vacuum gage, and several ther- 
mocouple vacuum gages. Commercial diaphragm com- 
pressors boost the gas pressure from the supply cylinder 
pressure to about 600 bars for storage in high pressure 
receivers. 





FIGURE 2. Burnett PVT apparatus. 






































FIGURE 3. Cross section of the Burnett volumes and valves. 


Carefully selected commercial nonrotating-stem 
valves are used; our experience shows that for a gas 
system the sealing surfaces of rotating-stem valves 
tend to gall and consequently to leak. The condition 
of the expansion and the filling-evacuation valves is 
particularly important; minimum deformation of the 
valve seat, long periods of trouble-free use, and 
certainty of valve closure are essential characteristics. 
A small valve with a block of 316 stainless steel, 
ports to accommodate 0.0032 m o.d. capillary tubing, 
and a stem of 17-4 stainless steel hardened to 40 
R/C is used for each of these two valves. Valve blocks 
for all the valves are selected on the basis of the 
concentricity of the seat and the stem-packing section. 
The shoulder of the valve seat is beveled slightly 
and carefully lapped to the same angle as that of the 
conica! end of the valve stem. The sealing surface is 
kept small to minimize the torque required for valve 
closure. To prolong valve life, all valves are closed 
with a torque wrench, and the required torque is 
recorded on each valve block. 








2.2. Pressure Measurement System 


a. Transducer 


A rugged yet highly sensitive and reproducible 
null-type pressure transducer separates the sample 
gas from the pressure measurement system. The 
transducer is shown in figure 4; an earlier model has 
been described previously [10]. The transducer 
consists basically of a prestressed metal diaphragm 
as the pressure responsive element, a capacitance 
sensor for detecting deflections of the diaphragm as 
small as 0.025 wm or less, and diaphragm backing 
surfaces incorporated as part of both the base and 
the capacitance sensor. 

The diaphragm is plane and circular, 25 wm thick, 
and prestressed to reduce inelastic behavior. The 
backing surfaces, spherically concave with a chord of 
0.0238 m that determines the effective diameter of 
the diaphragm and a sagitta of 25 um, are necessary 
because the stepwise expansions in the Burnett 
method result in temporary pressure imbalances 
across the diaphragm many times greater than the 
elastic range of a diaphragm with the required 
sensitivity. 

The earlier model has been modified to improve its 
performance, especially above 125 °C. A metallized 
ceramic disk replaces the former metal-epoxy capac- 
itance sensor assembly. This has enabled us to operate 
the apparatus at temperatures up to 225 °C, and the 
disk could even be used at temperatures as high as 
500 or 600 °C. The diaphragm is now fabricated from 
rolled nickel, instead of Monel, sheet stock; its thermal 
expansion thus closely matches that of the nickel 
base, and the diaphragm tension does not depend 
significantly on the temperature. A _polytetrafluoro- 
ethylené ring, prestressed by the clamping of the top 
and the base of the transducer and supported by 
inner and outer metal extrusion rings, replaces the 
former V-ring as the outer pressure seal. The dia- 
phragm sealing ring has been redesigned so that it 
stretches the diaphragm radially as it forms a pressure 
seal and thus does not distort the diaphragm when 
the transducer is assembled. Consequently the sample 
volume is more precisely defined and the transducer 
behaves more reproducibly. An important source of 
error is thereby reduced significantly. 

The transducer capacitance, consisting of the capac- 
itance contributed by the capacitor plate, diaphragm, 
and intervening gas plus the residual capacitance, 
constitutes one arm of a capacitance bridge. For 
each temperature and pressure the bridge is initially 
adjusted so that it is balanced when the pressure on 
each side of the diaphragm is the same. A subsequent 
pressure imbalance deflects the diaphragm from its 
null position and causes a bridge imbalance. The 


corresponding signal is amplified, phase detected, 
and indicated on a 100—-0-100 ywA output meter. 


The counterbalancing gas pressure is then adjusted 
until the diaphragm returns to its null position as 
indicated by the bridge balance. 
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FIGURE 4. Cross section of null-type pressure transducer. 
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A. Follower, steel 

B. Extrusion rings, lava stone 

C. Top, nickel 

D. Clamping bolts for pressure seal K 

E. Ring, ceramic 

F. Capacitor lead stop, silver-brazed bead 

G. Clamping ring for capacitance disk 

H. Inner clamping screws for clamping ring G 

I Outer clamping screws for clamping ring J 

J Clamping ring 

K. —_ Polytetrafluoroethylene ring with brass extrusion rings 
L. Capacitance disk 

M. Capacitor plate, metallized 

N Diaphragm, nickel 

O. Capacitor lead, drill rod 

P. Electrical and oil shield, Monel 

Q. Counterbalancing gas inlet tube, nickel 

R. Packing gland nut, steel 

S. Electrical insulating sleeve, polytetrafluoroethylene 
T. Extrusion rings, ceramic 

U. Packing rings. polytetrafluoroethylene 

V. Metallized surface to capacitor plate 

W. Thrust screws for clamping capacitance disk 
xX Thrust screws for diaphragm pressure seal A-A 
Y. Thrust screws for stressing diaphragm N 

Z.  Stressing ring, nickel 

4-—A. Diaphragm seal 

B-B. Diaphragm stressing-sealing ring 

C-—C, Diaphragm backing surfaces 

D-D. Base 

E-E. Sample gas port 
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Thé temperature dependence of the transducer 
null capacitance is negligible for temperature changes 
even two orders of magnitude greater than the fluc- 
tuations of 0.00] °C about the isotherm temperature 
that represent the resolution of the temperature con- 
trol system. For each isotherm the transducer is 
calibrated before Burnett measurements are made: 
The counterbalancing gas, which has been argon in 
all of our work, is admitted to both sides of the dia- 
phragm, and the null capacitance is measured as a 
function of pressure from atmospheric pressure up to 
the highest Burnett pressure to be measured. The 
pressure dependence found does not represent any 
significant change in the diaphragm null position but 
is caused primarily by the change in the dielectric 
constant of the counterbalancing gas. The pressure 
sensitivity of the transducer, as used by us in the 
Burnett apparatus for pressures up to 250 bars, ex- 
ceeds the sensitivity of the piston gages and the 
barometer used to measure the pressure. 


b. Piston Gages and Barometer 


Although Bourdon gages are used, for convenience, 
as pressure indicators at various locations in the ap- 
paratus, the basic pressure measurements are made 
with piston gages and a standard barometer — instru- 
ments which are primary pressure standards. Two 
controlled-clearance piston gages with an accuracy 
of better than 1 part in 20,000 cover the range from 
about 2 to 3 bars gage pressure up to 250 bars and 
beyond. A gage with gas as the pressure fluid covers 
the lower part of the range, and a gage with oil as the 
pressure fluid covers the higher part. Both gages cover 
the region from about 30 to 45 bars and in this overlap 
region they agree with each other to within | part in 
40,000 to 50,000 of the measured pressure. 

The nominal effective area of the oil gage piston is 
26 X 10-* m? and of the gas gage piston, about 
340 X10-* m?. The pistons are fabricated from tungsten 
carbide; the cylinders, from stainless steel. The top 
of the oil gage piston is hemispherical and fits into a 
conical recess in the coupling that connects the piston 
to the weight holder on which the weights are loaded. 
In the gas gage a lubricated ball bearing fits into a 
conical recess in both the piston and the coupling. In 
each case the purpose of the spherical contact, which 
represents our modification of the gages, is to reduce 
the amount of side thrust applied to the piston. Simi- 
larly, the weight holders and the weights that nest on 
them are machined carefully to avoid eccentric loading. 
The prevention of side thrust is particularly important 
for the gas gage because its piston receives little lateral 
support. The gas gage piston and cylinder must also 
be kept scrupulously clean. For each gage the load 
applied to the piston is spun manually and continues 
to spin for up to an hour or more before it stops. 

Even though the pistons and cylinders were care- 
fully lapped and cleaned before they were delivered 
to us, we were able to improve the performance of the 
gages by using a simple motor drive assembly to rotate 
the piston and move it up and down within its cylinder 





for several hours while the piston and cylinder were 
flooded with mineral spirits. The resultant burnishing 
increased the spin time, and the flooding helped to 
remove any fine particles that might have been present. 
The burnishing procedure can be used to improve 
performance at any time, particularly if the cylinder 
has become slightly scratched. 

A simple piston gage without the controlled-clear- 
ance capability is used to make less accurate meas- 
urements over the range from 2 to 3 bars gage pressure 
down to 0.2 bars absolute pressure or less: the abso- 
lute measurements are made under an evacuated 
bell jar. The gage has three stainless steel piston and 
cylinder assemblies with nominal effective areas of 
800 x 10-* m?, 86 x 10-® m?, and 8 X 10~* m?, and uses 
gas as the pressure fluid. 

The atmospheric pressure is measured with a stand- 
ard barometer of modified U-tube design consisting of 
a column and a cistern. The level of the mercury in the 
column is detected photoelectrically and is measured 
with a 31 in (0.79 m) scale and a vernier graduated to 
0.001 in (25 xm); the level of the mercury in the cistern 
is detected electrically and is measured with a microm- 
eter graduated to 0.0001 in (2.5 wm). The back pressure 
above the mercury in the column is measured with a 
McLeod gage and is generally about 4 <10-® bars or 
less. “he accuracy of the barometer scale is easily 
checkec vith gage blocks. Repeated tests have shown 
that the oarometric pressure measurements are ac- 
curate to 33 X10~® bars. 

c. Other Components 


Small adjustments in the counterbalancing gas pres- 
sure are made with a compound gas injector whose 


MERCURY HEIGHT 
SENSOR—__ 





FIGURE 5. 





internal volume is changed by the manual displace- 
ment of a 0.013 m or 0.0032 m diam piston. The coun- 
terbalancing pressure acts directly on the gas gages but 
is separated from the oil gage either by a commercial 
diaphragm transducer or by a mercury U-tube. The 
U-tube is an adaptation of an earlier model by Keyes 
[11] and is shown in figure 5. To reduce capillary effects, 
the lower part of each bore containing the mercury is 
made large in diameter. The bore diameter is reduced 
slightly above the mercury to make the gas volume in 
the bore a small fraction of the total volume of the 
counterbalancing gas. As a result, the bore gas volume 
has only a small effect on the performance of the gas 
injector. The mercury level in each column is indicated 
by electrical contact between the mercury and a fine 
metal probe whose height is initially adjusted to within 
25 um of the reference height and is checked period- 
ically. The modified design permits the valve at the 
base of the U-tube to be replaced conveniently. 


2.3. Temperature Control and Measurement System 


An oil bath and a servo heater system are used to 
control the temperature of the Burnett apparatus. As 
shown in figure 6, the Burnett vessel, the transducer, 
the expansion and the filling-evacuation valves, and the 
interconnecting lines are all immersed in oil in an inner 
chamber located in a thermally insulated bath. The 
inner chamber is double-walled, and the space between 
the walls can be evacuated to reduce radial heat flow. 
An oil-circulating pump at the bottom of the bath is 
coupled mechanically to a motor mounted on the top. 
Oil drawn into the manifold and the pump is ejected 
upward into the inner chamber and flows by the Burnett 








Vercury U-tube. 
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FIGURE 6. Oil bath. 


apparatus, over the top of the inner chamber, and then 
downward past the heat exchanger to complete the 
cycle. Sufficient mixing of the oil occurs in the mani- 
fold and the pump to eliminate temperature gradients 
acquired during the flow of the oil past the heat 
exchanger. 

The heat exchanger consists of a cooling coil and a 
series-parallel heater for low and high power output. 
A coolant is circulated through the coil and a refrig- 
erated bath, which is thermostated to within 0.003 °C 
from —20 to 5 °C. The cooling system is used for oil 
bath temperatures from 0 °C to slightly above room 
temperature. 

The servo heater system consists of an adjustable 
Wheatstone bridge with a resistance thermometer con- 
stituting one arm, a d-c amplifier and recorder, a three- 
mode controller for proportional, derivative, and inte- 
gral response to the error signal, and a power amplifier 
connected to a heater circuit. The resolution of the tem- 
perature control is 0.001 °C; slight differences from the 
desired temperature are readily corrected by the use 
of a variable resistor in the bridge circuit. 

The temperature as defined by the International 
Practical Temperature Scale of 1968 is measured with 
a capsule platinum resistance thermometer located in 
the well between the sample and the expansion cham- 
bers in the Burnett vessel. The thermometer leads are 
brought out through a protection tube filled with oil to 
ensure good thermal contact between the thermometer 
and the Burnett vessel. A Mueller bridge with a photo- 
electric galvanometer as the detector is used to measure 
the resistance. 


Differential five-junction copper-constantan thermo- 
couples are used to measure temperature differences 
over the Burnett vessel and transducer as well as near 
the top of the oil bath. A small temperature difference 
can exist between the Burnett vessel and the top of 
the bath due to heat loss through the apparatus supports 
and oil-vapor loss to the surroundings. This difference 
is limited to 0.01 °C by means of a manually-controlled 
auxiliary heater at the top of the bath. Temperature 
differences over the Burnett vessel are less than 
0.002 °C. Electrical noise in the thermocouple circuits, 
which include the switching and measurement instru- 


mentation, corresponds to a temperature uncertainty 
of less than 0.0001 °C. 


2.4. Experimental Procedure 

An experimental procedure has been evolved which 
ensures the reproducibility of the pressure measure- 
ments, the precise definition of the sample volume, 
and the detection of any sample gas leakage. A series 
of runs is made at each isotherm temperature. Each 
run may include up to approximately eight expansions; 
the initial pressure for each run is chosen so as to 
give a well-distributed set of pressures for the iso- 
therm. Before any Burnett measurements are made, 
the oil bath is thermostated at the desired tempera- 
ture and the transducer is calibrated over the desired 
pressure range. The transducer null capacitance at 
atmospheric pressure is checked before and after 
each day’s measurements. 

The first step in a run is to evacuate the sample gas 
part of the apparatus and flush it thoroughly with the 
sample gas. After this part of the apparatus has been 
reevacuated, the sample and the expansion volumes 
are filled with fresh sample gas to the highest pressure 
desired for the run. In the flushing and the filling proc- 
ess and throughout the run, any pressure imbalance 
across the transducer diaphragm greater than a few 
bars is directed toward the base, which is more rugged 
than the capacitance sensor and its clamping arrange- 
ment and hence provides better diaphragm support. 
This precaution helps to ensure the stability of the 
diaphragm null position and thus aids in the precise 
definition of the sample volume. Similarly, as a pres- 
sure measurement is made, the diaphragm is slightly 
overloaded toward the base prior to the final balancing 
of the counterbalancing and the sample gas pressures. 

After the filling process is completed and the filling- 
evacuation valve closed, sufficient time is allowed 
for temperature and pressure equilibration. During 
the course of a run the temperature of the Burnett 
vessel may vary over a range of about 0.002 °C due to 
instabilities primarily in the temperature control 
bridge; equilibration times are chosen to keep short- 
term temperature changes well within this range and 
temperature differences within the ranges discussed 
earlier. 

While the equilibration is taking place, preparations 
are made for the first set of Burnett measurements. 
The counterbalancing gas, which was maintained at 
a higher pressure than the sample gas during the filling 
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process, is gradually vented until the transducer null 
indicator returns to its Zero position. The approximate 
counterbalancing gas pressure is then read from a 
Bourdon gage, and a variable capacitor in the capaci- 
tance bridge is adjusted in accordance with the 
transducer calibration so that the bridge balance will 
truly correspond to the diaphragm null position. The 
piston gage to be used for the precise pressure meas- 
urement is loaded with approximately the right 
combination of weights that will be needed. 

When the resistance thermometer and the differen- 
tial thermocouple readings show that the equilibration 
is essentially complete, the expansion valve is closed 
to isolate the sample and the expansion volumes. While 
the valve is slowly closed, the diaphragm is maintained 
in its null position by adjustment of the counterbalanc- 
ing gas pressure with the compound gas injector. An 
amplifier and supplementary null indicator about ten 
times as sensitive as the transducer output meter are 
used in this process, partly because of their greater 
sensitivity and also because the location of the supple- 
mentary null indicator is more convenient to the person 
closing the valve and operating the compound gas 
injector. 

After the expansion valve is closed, the first set of 
Burnett measurements can be made. In the operation 
of the piston gage a precise zeroing of the transducer 
output meter would require an unnecessarily long time 
to be achieved. Instead the gage operator calibrates 
the output meter over its central range by noting the 
difference in the scale readings that correspond to two 
nearly equal loads on the gage. The slight difference in 
the diaphragm position represented by the different 
scale readings changes the sample volume, which is 
bounded by the diaphragm, by a negligible amount — 
less than 2 ppm. The scale readings are made when 
the piston is at its reference height as shown by an elec- 
tronic proximity sensor mounted underneath the weights 
loaded on the gage. The height of a gas gage piston is 
adjusted by means of the compound gas injector: the 
height of the oil gage piston is adjusted by means of a 
separate oil injector. For either type of gage the piston 
is raised above its reference height by means of the ap- 
propriate injector and is then allowed to fall slowly down 
to its reference height while the piston and its load are 
spinning steadily. The gage pressure and barometric 
pressure measurements are made almost simultane- 
ously, and the resistance thermometer and differential 
thermocouple measurements are made immediately 
afterward. These measurements constitute the first 
set of Burnett measurements. 

Next the expansion valve is opened and after about 
a minute is closed as before. A second set of Burnett 
measurements is made. Since the expansion step has 
not been performed, this second set of measurements 
should be in close agreement with the first set. If the 
second set, in fact, does not differ significantly from the 
first, the filling-evacuation valve is opened and the ex- 
pansion volume is evacuated. After isothermal condi- 
tions are reestablished and a thermocouple vacuum 
gage, located as close to the expansion volume as is 
convenient, shows a pressure of about 4 < 10~-® bars, a 





third set of measurements is made. If this set does not 
differ significantly from the first two, the _filling- 
evacuation valve is closed and the expansion valve is 
opened. The gas initially in the sample volume ex- 
pands to fill the expansion volume as well. After this 
first expansion the measurement procedure is repeated 
to successively lower pressures. 

To be considered satisfactory, the three sets of meas- 
urements corresponding to the initial filling or a sub- 
sequent expansion must show a pressure reproducibility 
of 1 part in 50,000 or better. Occasionally a set of meas- 
urements will be repeated if some temporary disturb- 
ance, such as a temperature or barometric fluctuation, 
has caused a poorer reproducibility. If a serious ex- 
perimental failure, such as sample gas leakage, were 
to occur, the measurements for that run would not be 
used in the data reduction. 


3. Data Reduction and Discussion of 
Results 


3.1. Background 


In contrast to the simplicity of the experimental 
measurements, the reduction of Burnett data involves 
rather formidable nonlinear statistical analysis of 
relationships containing the compressibility factor, 
defined by a pressure or a density virial expansion, 
and optional constants. The authors have previously 
presented a general discussion of the analysis of 
different Burnett relationships using two nonlinear 
statistical methods, due to Gauss and to Deming [9]. 


Related papers have also been published by others 
{12, 13]. 


Our experience indicates that compressibility 
factors and second virial coefficients calculated from 
Burnett data of! high precision and accuracy are as 
accurate as those calculated from the best data ob- 
tained with any other PVT method. Various treatments 
of relationships, that is, various models, based on eqs 
(7) and (8) have been used in our analysis. In these 
models the cell constant, where it appears, is regarded 
as a single parameter common to all the experimental 
runs, as a separate parameter for each run, or as a 
constant to be determined from Burnett measurements 
for a nearly ideal gas. Meaningful comparisons of 
the results are difficult for at least two reasons: (1) 
the same systematic errors may propagate differently in 
each model and (2) the number of parameters to be 
evaluated varies from one model to another. The 
effect of the latter is subtle; it may introduce an 
analytical bias, as the magnitude of the correlation 
between the parameters, which is basically high, 
increases significantly with an increase in the number 
of parameters. It is not surprising, therefore, that slight 
discrepancies exist among the results of different 
models, even when applied to the same data. For 
example, a comparison of the compressibility factors 
derived from argon data taken with our Burnett 
apparatus shows discrepancies among the values of 
approximately 0.0] percent to 0.04 percent, which, 
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although small, are nevertheless greater than the 
estimated statistical uncertainties of less than 0.01 
percent. For the six models considered, the results of 
four show comparative discrepancies of about 0.0] 
percent. Another example of this type of discrepancy 
occurred in an informal round-robin analysis of helium 
data by participants at a two-day Burnett conference 
held at the National Bureau of Standards in February 
1969. Initial analyses indicated discrepancies of 0.04 
percent among the results of the different models as 
compared to statistical uncertainties of 0.0] percent 
in the results themselves. 


We shall not repeat our previous detailed presenta- 
tion of the nonlinear statistical analysis of Burnett 
data [9], but shall present a brief discussion which 
emphasizes the important differences between the 
various models in order to help one choose the models 
that introduce the least bias. The fact that the data can 
be reduced in several ways is an attractive feature of 
the Burnett method and enables one to study the dif- 
ferent effects of systematic errors. 


The basic Burnett equations may be transformed into 
different relationships among the variables of interest. 
In terms of a finite virial expansion in pressure, 


I 
Z=1 a 2 biP', eq (7) may be expressed as 


(+> bP.) 














Pru i=1 
P =N, ; : (9) 
4 i+y bP) 
i=1 
P,(NiNo . . . Nr) = A (1+ bP), 0) 
r wV¥2- +--+ 4¥r pH Pan a 
P,-2P, aaa l + AN,-1(Pr-2, P,_1) 
i 1+ AN,(Py-1, Pr) 
(1+ 2 bPt_.)(I +> Pt) 
x — > (11) 
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In terms of a finite virial expansion in density, 


J 
Z=1+ Ss Bip), eq (8) may be expressed as 


j=1 


—— "ie /___...___¥ 
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It should be noted that in these equations the second 
virial coefficient is represented by 6; or B;. For both 
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the pressure and the density relationships 
N,=N{(l + AN,(P r-l> P,)] 
and 


N,N, ee N,=N'[(1+ AN, (Po, P;)| 


x [1+ AN2(P1,P2)]... (1+ AN,(P,y-1, Pr)) - 
Any one of these relationships may be applied to 
Burnett data to form an overdetermined set of equa- 
tions for evaluating the virial coefficients and the 
other parameters. The latter consist of the cell con- 
stant (except in eq (11) where the cell constant has 
been eliminated), the initial density po in eq (12), 
and sometimes the quantity Po/Zo in eq (10). The latter 
two parameters are characteristic of an individual 
run. As mentioned in the introduction, the pressure 
corrections (AN,, ANo,...) are not considered as 
parameters but are determined independently of the 
parameter evaluation. 


3.2. Compressibility Factor of Argon at 25 °C 


Compressibility factor results derived from the 
different relationships applied to our argon data are 
given in table 1 [9}| The models are designated briefly 
at the top of each column. For each model the param- 
eters were evaluated simultaneously for all the data. 
Where a fixed cell constant was used, it was pre- 
determined from eq (9) applied to helium data taken 
at the same temperature for pressures below 35 bars. 
For these conditions helium behaves nearly as an 
ideal gas, a characteristic which simplifies the data 
reduction. 

The results at each pressure are essentially the same 
except for those derived from the density relationship 
with an adjustable cell constant included as either 
a separate parameter for each run or a single param- 
eter common to all the runs. The discrepancies among 
the first four columns, which are about 0.01 percent, 
could be due to the presence ot systematic error which 
propagates differently in each model. A likely source 
of such error is the correction for the volume dis- 
tortion of the sample and expansion chambers with 
pressure. A study of simulated Burnett problems 
based on our argon data indicated that other probable 
sources of error had a negligible effect on the com- 
pressibility factor. In this study a simulated syste- 
matic error of 1 part in 30,000 was superimposed on 
the argon pressure data. The density reduction of 
these data, with the cell constant either fixed or re- 
garded as a parameter, resulted in a change in the 
compressibility factor of less than 0.001 percent for 
all the pressures. An effect of less than 0.001 percent 
was also obtained when a simulated systematic error 
of 1 part in 20,000 was superimposed on the RT product 
in a density reduction of the data. 

The volume distortion correction was determined 
from distortion equations and strain gage measure- 
ments of the external deformation of the Burnett 


TABLE 1. Comparison of results derived from Rurnett argon data at 25 °C and to 250 bars*: 























3 Runs (7, 7, & 6 expansions/run) 
Pressure series formulation Density series formulation 
. ’ | _ Me 
Eq (9) Eq (9) Eq (11) Eq (12) Eg (12) Eq (12) 
Cell constant One cell constant Cell constant Cell constant One cell constant All cell constants 
fixed adjusted eliminated fixed | adjusted adjusted 
N 1.74505 1.74502 1.74505 | 1.74511 1.74511 
1.74509 
| 1.74511 
Zon208 wens 0.99384 0.99387 0.99388 0.99380 0.99375 0.99376 
Z101.33 pars .95356 .95365 .95367 .95352 | .95333 .95334 
abaen .95202 .95209 .95212 .95194 .95169 .95173 

















* The compressibility factor is shown to five digits to indicate the differences in the numerical results and not to indicate the overall 


accurecy. 
b 1 bar= 105 N- m=? 


vessel due to internal pressure [14, 15]. The effect 
of the distortion correction and its estimated uncer- 
tainty on the compressibility factor is shown in figure 
7 for argon data analyzed with the use of the pressure 
equation (10). The effect of the uncertainty is small, 
about 0.01 percent at the maximum pressure of 250 
bars. The results of simulated problems indicate that 
the effect is the same in the density formulation. In 
the pressure series formulation the effect of the dis- 
tortion correction uncertainty is consistent with the 
small differences between the results derived with the 
cell constant fixed and those derived with the cell 
constant evaluated or eliminated. 

We now consider the basis for the differences, 
particularly at the higher pressures, which are sig- 
nificantly greater than 0.01 percent. These are between 
the results derived from the pressure series formula- 
tion and the results from the density series formula- 
tion, particularly with the cell constant included as 
either a separate parameter for each run or a single 
parameter common to all the runs. If the exact model 
were known, that is, if all the systematic errors 
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FIGURE 7. Relative compressibility factor correction due to pres- 


sure distortion of Burnett cell volumes, argon run. 


were included in each equation, then the uncertainties 
in the values of the compressibility factor, the cell 
constant, and the initial densities would be consistent 
with the magnitude of any random errors in the pres- 
sure. This has been illustrated with a simulated Burnett 
problem for a particular model in which random errors 
were superimposed on the pressure data [16]. In 
the real world one cannot account exactly for all of 
the systematic errors. Consequently, a bias is intro- 
duced in the statistical analysis, and the results are 
not as reliable as the statistical uncertainties indicate. 
In addition, the number of parameters varies from one 
model to another. This variation adds to the bias since 
the absolute value of the correlation between the 
parameters increases significantly with an increase 
in the number of parameters. 

In general the correlation between virial coefficients 
is high regardless of which PVT method is used to 
determine them. In the Burnett method additional 
highly correlated parameters are also present. In the 
density reduction of the argon data with a fixed cell 
constant, the absolute value of the computed correla- 
tion between each of the virial coefficients is about 
0.85 on a normalized basis of 0 to 1, whereas it is about 
0.6 between each of the initial densities and about 
0.75 between any one of the virial coefficients and any 
one of the initial densities. In contrast, if the cell 
constant is included as a parameter, the statistical 
fit is improved over that for the fixed cell constant, 
in the sense that the sum of squares is reduced. 
However, the correlation between any two of the 
parameters increases significantly! The correlation 
between each of the virial coefficients becomes about 
0.93; between each of the initial densities, about 
0.994; and between any one of the virial coefficients 
and any one of the initial densities, about 0.87. The 
correlation between the cell constant and any one of 
the virial coefficients is about 0.88 and between the 
cell constant and any one of the initial densities, about 
0.995. Exceptionally high correlations such as 0.994 
and 0.995 indicate that these parameters are not truly 
independent and thus cannot be accurately evaluated 


[17]. 
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Additional data, not highly correlated to the data 
already in the analysis, would reduce the high param- 
eter correlation. As a result, the evaluation of the 
parameters would be enhanced and any analytical 
bias in the results would be reduced. In one approach 
the analysis would consist of data from the same ap- 
paratus for two or more sample gases with different 
compressibility characteristics. The data would be 
taken at the same temperature, and a single cell 
constant parameter, common to all the data, would be 
evaluated. Another approach would be to use a fixed 
cell constant obtained from helium data. This approach 
is less desirable because the cell constant would be 
entered without any statistical uncertainty, that is, as an 
exact quantity, and hence the computed statistical 
uncertainties would be too small. The significance of 
high parameter correlation should be considered care- 
fully in estimating the reliability of results obtained 
from the density reduction methods. For the pressure 
reduction methods the situation is similar, although the 
parameter correlation is not as high. 

Our values of the compressibility factor shown in 
table 2 represent the average of the results of the three 
pressure models and the density model with a fixed 
cell constant; the results of the density models with one 
or three separate cell constants were not included in 
the average because they were considered to be less 
reliable due to the high correlation between the 
parameters. For pressures up to 250 bars the overall 
accuracy of the averaged compressibility factor is 
estimated to be within + 0.0002 on the basis of our error 
analysis. The statistical uncertainty of each value of the 
compressibility factor is no greater than one sixth 
of the overall accuracy of the averaged compressibility 
factor. Table 2 also includes comparative values of 
the argon compressibility factor derived from the 
piezometer results of Michels [18] for pressures up 
to 202.7 bars. The maximum difference is small, only 
0.05 percent at 202.7 bars; PVT results are generally 
considered to be accurate if the overall error is less 
than 0.1 percent. 























TABLE 2. Comparative values of the compressibility 
factor for argon at 25°C* 
Pressure This paper” Michels‘ 
bars Compressibility factor Compressibility factor 
10.133 0.9938 0.9938 
101.33 .9536 .9534 
202.65 .9520 9515 
*] bar= 10 N- m*?. 
® The overall accuracy is estimated to be + 0.0002. 


{18} 


3.3. Second Virial coefficient of Carbon Dioxide 


at 50 °C 


In addition to accurate values of the compressibil- 
ity factor, reliable second virial coefficients are also 
obtainable from the Burnett method. Our own ex- 
perience with argon at 25°C and with carbon dioxide 
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at 50°C, 19 degrees above its critical temperature, 
shows that this is true over a wide range of reduced 
temperatures greater than unity.2 The argon results 
are for a reduced temperature of 1.98 and have been 
discussed previously [9]; the carbon dioxide results, 
which are discussed in the following paragraphs, 
are new and are for a reduced temperature of 1.06. 

Two groups of carbon dioxide measurements were 
made: one group over the pressure range from 35 to 
4.5 bars and another group over the range from 4.5 
to 0.5 bars. The higher pressure data are the more 
accurate data and were reduced by the density model 
with a fixed cell constant determined from helium data 
at the same temperature. From the more accurate 
data a second virial coefficient for carbon dioxide 
of —102.25+0.04 cm*/mol was obtained, where the 
uncertainty is only the statistical uncertainty and is 
necessarily too small because it does not include the 
uncertainty in the cell constant. Our estimate of the 
overall systematic uncertainty, based on our error 
analysis, is 0.2 cm*/mol. This value of the second 
virial coefficient was supported by a graphical evalu- 
ation which included the less accurate carbon dioxide 
data for the lower pressures. It compares favorably 
with the value of —102.5 cm*/mol obtained from a 
recent correlation [19] which included an evaluation 
of all existing data but is in substantial disagreement 
with the value of —100.7+0.4 cm?/mol determined 
by Dadson and his co-workers [20, 21]. It should be not- 
ed that Dadson et al. carried out their measurements 
with both the Burnett method and the piezometer 
method. The fact that their results from both these 
methods agreed with each other and yet disagreed 
with the results of the correlation prompted our own 
measurements. Our results and more detailed compari- 
sons will be the subject of a later publication. 

Various potential sources of systematic error were 
investigated in our own work. In simulated problems 
based on the higher pressure data, superimposed 
systematic errors of +25 percent in the distortion 
correction and | part in 20,000 in the RT product had 
only a negligible effect on the second virial coefficient. 
The uncertainty in the distortion correction would be 
expected to have a negligible effect since the distortion 
correction itself is very small even at the highest 
pressures in the carbon dioxide data. A cell constant 
was computed for the carbon dioxide data in a separate 
calculation by treating it as an adjustable parameter 
in the density series formulation. The value obtained 
agreed with the predetermined value obtained from 
helium data to within | part in 20,000, which is excellent 
agreement in view of the large magnitude of the non- 
ideality of carbon dioxide at 50°C (T/T criticai= 1-06) and 
of the extremely high correlation between the param- 
eters in the analysis. 

The possibility of adsorption of carbon dioxide on the 
surfaces of the sample and the expansion volumes was 
investigated by means of a graph of the low pressure 
data in which the pressure ratio P,_,/P, was plotted 


2 The reduced temperature is the ratio of the sample gas temperature to its critical 
temperature, where both these temperatures are expressed on the absolute scale. 





as a function of the pressure P,_,. The significance 
of the graph is based on the equation 


Pee 
p, N+ h(N— DPA, 


(13) 
which is the low pressure form of eq (9). The extrapo- 
lation of the curve drawn through the points to the verti- 
cal axis yielded an intercept whose value, the cell con- 
stant, was essentially the same as the predetermined 
value obtained from the helium data. This agreement 
would not be expected if adsorption were present. 
For a nonpolar gas such as carbon dioxide and for our 
range of pressures, adsorption is expected to be sat- 
urated [22]; that is, above a relatively low pressure 
very little more adsorption occurs as the pressure is 
increased further. In the sample volume the adsorbed 
gas would remain essentially unchanged upon an ex- 
pansion, whereas in the expansion volume gas would be 
desorbed in the evacuation process and adsorbed 
in the expansion process. As a result, the number 
of free molecules would be reduced upon an ex- 
pansion, and the ratio P,_,/P,; would therefore in- 
crease rapidly at low pressures. The extrapolation 
to zero pressure would yield a value for the cell 
constant that is greater than the value predetermined 
from helium data for which adsorption is not signifi- 
cant. Furthermore, the graph of this ratio versus pres- 
sure would show increasing nonlinearity at low pres- 
sures instead of the decreasing nonlinearity character- 
istic of real gases as the effect of the third and higher 
order virial coefficients not shown in eq (13) becomes 
insignificant. Neither of these signs of adsorption is 
present in our data. Hence we feel confident that any 
adsorption in our apparatus is so slight as to have an 
insignificant effect on our results. 


We gratefully acknowledge the benefit we have 
received from general discussions with M. Klein and 
from comments regarding adsorption made by P. T. 


Eubank. 
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Sliding short methods represent a measurement tool of substantial potential for the measurement 
of small losses such as are associated with waveguide connectors or adaptors. Until recently, however, 
the use of these methods has been inhibited by the uncertainty of the error contribution due to non- 


ideal short behavior. 


A recent analysis by Almassy has shown that by the use of proper techniques, the error contribu- 
tion from this source is usually negligible, provided that the adaptor (or connector) is “well matched.” 
It is the purpose of this paper to eliminate this latter restriction, develop additional measurement 


methods, and describe further applications. 


Key words: Adaptor; connector: efficiency; sliding short. 


1. Introduction 


A measurement problem, of continuing interest in 
the microwave art, is that of adaptor evaluation. 
With the advent of the power equation [{1|' methods 
and use of terminal invariant parameters [2], the 
dependence upon impedance properties is suppressed, 
and the dissipative characteristic emerges as the 
parameter of major interest. A similar observation 
may be made with regard to connectors. Here the dissi- 
pation at the connector interface, or more specifically, 
the lack of loss repeatability, represents a basic limi- 
tation to the attainable accuracy in many microwave 
measurements. 

The application of sliding short methods, to this 
measurement problem, has been known for some time 
(3, 4]. These techniques were generalized in conjunc- 
tion with the development of the power equation 
methods [1, 2]. Despite this revival of interest, how- 
ever, the technique has been subject to a major limi- 
tation: an ideal sliding short is assumed. 

In practice, of course, no such device exists. 
Moreover, the losses in the short may be of the same 
order (or larger) than those in the joint or adaptor to 
be evaluated. In the absence of further information, 
this loss represents a potential source of substantial 
error. 

A recent analysis by Almassy [5] has shown that, 
with proper techniques, the sliding short losses can 
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be substantially reduced, provided that the adaptor 
(or joint) reflection is small (i.e., |Si:|<1, |S22|<1). 

It is the purpose of this paper to eliminate this 
restriction, develop additional measurement methods, 
and describe further applications. 





2. Background 


The major thrust of this paper will be directed 
towards the measurement of the efficiency of an arbi- 
trary two-port. In particular this two-port may be a 
connector, waveguide-coax adaptor, etc. By definition, 
the efficiency is the ratio of the (net) power output, 
to the (net) power input, and is a function of the load 
(but not of the source) impedance. Throughout this 
paper, the efficiency is assumed to be large (small 
losses). (In general the accuracy of this method de- 
creases rapidly with losses greater than 10 dB.) 

It will prove convenient to briefly review the existing 
theory. The recommended instrumentation takes the 
form shown in figure 1. The basic configuration will 
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Generalized reflectometer for use in efficiency 
measurement. 


FIGURE 1. 
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be recognized as a generalized reflectometer (g-reflec- 
tometer), where the junction parameters are arbitrary 
except as subsequently noted. In particular, the 
tuning transformers 7,, 7, permit one to impose 
certain conditions on the measurement system. The 
device to be evaluated is connected to port 2, while 
arms 3 and 4 are terminated by power meters. 

Although power meters have been specified for arms 
3 and 4, it will prove convenient, momentarily, to as- 
sume an alternative detection scheme such that the 
complex ratio b3/bs is measured, where 63, b, are the 
emergent wave amplitudes. Let arm 2 be terminated 
by a moving short (T2=e/®), and let the values 


(3) be plotted in the complex plane as @ is permitted 
4 


to vary. The resultant locus is a circle, an example of 

which is shown in figure 2. The parameters of this 

locus, in particular the radius of the circle, R, and the 

(absolute) distance, R,, between its center and the 

origin, play a major role in the efficiency measurement. 
Inspection of figure 2 indicates that ” 

















SP] =R+Re, (1) 
bs max 
Fa =R—-R.. (2) 
| bs | min 
and solving for R, R- yields: 
‘ =3( bs max bs min ) (3) 
_1(|bs| | bs 
~ 2 ( bs max bs min ) (4) 




















FIGURE 2. 


Locus of b3/bs as a function of 8. 


2 It is assumed that the origin is within the circle. As long as the efhciency is high, this 
condition is assured. As the losses increase, however, it is possible that R- > R. For a 
more complete discussion see [2]. 





The determination of R, R. thus involves amplitudes 
only, the phase detection capability is not required. 
Finally, the power meters respond to the square of the 
amplitude such that P3 = | b3 |? ete. 

In order to make an efficiency measurement, the 
two-port is connected as shown in figure 3. Note that 
the designations have been chosen such that port 2 
of the four arm junction mates with port 2 of the two- 
port. The “load” impedance, for which the efficiency 
is measured, is actually provided by the four arm 
junction. 

To be more specific, it is convenient to postulate 
that P, is canstant. (In practice this is often done by a 
“leveling” or feedback arrangement.) The source im- 
pedance, I,, for the “equivalent” generator, which 
now obtains at port 2, depends upon 7, and the ad- 
jacent coupler [6], but is independent of 7, and the 
remaining coupler. The efficiency, which is obtained 
in the measurement described in the following para- 
graph, is that for an assumed power flow from terminal 
1 to terminal 2, with the termination Ig. It will be 
denoted by yi (T,). 

The efficiency measurement calls for connecting a 
moving short to terminal 1 and determining the radius, 
R,, of the resulting circular locus using (3). The two- 
port is then removed, and the operation repeated at 
terminal 2 to obtain Ro. It has been shown [2] that 





Ne (Ty) = Ri/Re. (5) 


In an alternative method, it is convenient (but not 
essential) to use the configuration shown in figure 4. 
Here the efficiency, n:2(T1), is measured for a power 
flow from terminal 2 to terminal 1, and for the termi- 
nating load I’). 

The measurement procedure now calls for the ad- 
justment of 7, such that P; vanishes. The parameters, 
R, R-, of the circular locus are then measured with the 
moving short connected to terminals | and 2. It will be 
shown (see appendix) that 72(I;) is given by: 





FIGURE 3. G-reflectometer with two-port connected. 





FIGURE 4. Alternative configuration for g-reflectometer. 
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It should be noted that the adjustment of 7,, in 
both cases, is dictated by the value of load impedance 
for which the efficiency is measured. Provided that 
the indicated configurations of couplers and tuners are 
used, the operation, in theory, is independent of the 
adjustment of 7, Tuner T, is useful, however, be- 
cause of certain practical considerations which will 
be explained later. If the suggested coupler positions 
are interchanged, (i.e., the configuration of figure 4 
is used to measure 7e(I,), or the configuration of 
figure 3 is used to measure 72(I;)) the indicated role 
of T, is shared by 7, and the desired adjustments 
become interdependent and more difficult to realize. 


3. Implementation 


The application and implementation of these meth- 
ods is perhaps best described in terms of a specific 
problem—an efficiency measurement of a waveguide 
joint or connector. 

Referring to figure 5, it is convenient to first postulate 
a length of lossless waveguide lead, an ideal moving 
short, and an “ideal” reflectometer. (In terms of the 
description contained in (7), in the next section, an 
“ideal” reflectometer is one for which b=c=0.) The 
object is to measure the loss of the indicated waveguide 
joint. For reasons which will emerge, it is useful to re- 
cord the reflectometer response as a function of short 
position. For the “ideal” system of figure 5 the expected 
response is that shown in the inset. This result is ex- 
plained as follows. 

As long as the short is to the left of the waveguide 
joint, the reflection coefficient presented to the reflec- 
tometer is of unit magnitude and variable phase. The 
system responds only to the magnitude; this accounts 
for the “straight” line section. With the short to the 
right, however, the joint losses will lead to reduced val- 
ues for the reflection coefficient magnitude. In partic- 
ular, this loss is determined by the longitudinal current 
component and thus will be a maximum when the short 
position is a multiple of a half-wavelength, and a mini- 
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FIGURE 5. 


An “ideal” measurement system and its response to the 
efficiency measurement of a waveguide joint. 
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mum (zero*) when the position is an odd multiple of a 
quarter wavelength. In this case, the maximum and 
minimum values to the right of the joint are averaged, 
and the ratio of this average to the value left of the joint 
is the efficiency. 

If this experiment is carried out in reality, a typical 
record is that shown in figure 6. The most obvious 
change is in the downward slope, this is caused by the 
waveguide loss. The measurement procedure now calls 
for a projection of the maxima and minima to the plane 
of the joint and the parameters R, R, evaluated from 
these projections. The mathematical basis for this will 
be given in the following section. 

Although an “ideal” reflectometer was assumed in 
figure 5, in principle it is possible to start with arbitrary 
couplers and provide only for the tuning adjustments 
specified in the preceding section. In general, depend- 
ing upon the adjustment of 7,, the curve will show an 
oscillatory behavior on both sides of the joint. As a 
practical matter, it is generally desirable to adjust 7, 
such as to keep the amplitude of this oscillatory com- 
ponent within nominal limits. This point will be con- 
sidered in greater detail in a following paragraph. 

For an arbitrary two-port, it is not possible (in gen- 
eral) to make a continuous recording as shown here; 
in this case it is necessary to choose separate reference 
planes (usually coinciding with the waveguide flange) 
in the input and output arms; the determination of 
R, R- follows from the projections of the maxima and 
minima to the reference planes as already described. 
At this point it is necessary to decide whether the 
associated joint losses are included as part of the two- 
port. In many cases the joint loss, or at least its re- 
producibility, limits the accuracy to which the 
two-port efficiency may be specified. 


4. Analysis 


The validity of the foregoing procedure was estab- 
lished by Almassy [5] under the condition that the 
two-port is “well-matched,” (i.e., both |Si:| and |S22| 
are much smaller than unity). The purpose of the 
present analysis is to demonstrate that the method 
loses very little of its accuracy, even when this con- 
dition is not satisfied. 

In general the functional relationship between bs/b,, 
and a termination of reflection coefficient I’, at either 
terminal ] or terminal 2, is in the form: 


63 _ al'+b (7) 
by cl+1 
Here the complex parameters a, b, c depend upon 
both the reflectometer and the two-port when the 
termination is at terminal 1, but only on the former 
when the termination is at terminal 2. For an ideal 
short, [=e where @ is a function of short position. 


3 This assumes the equivalent circuit of the joint is that of a small series resistance. 





In order to account for non-ideal behavior, the model 
here assumed is: 


P= reve Xai (8) 


where r, & are the magnitude and phase of the plunger 
reflection, a, 8 are the loss and phase constants of 
the line in which the “short” moves, and / is its posi- 
tion with respect to the reference flange. 

If (8) is substituted in (7), the resulting relation can 
be written as follows: 


bs b—ac* re 
by 1—|cl?re-4*! 
(a— bc) re%e 2(a+jBN ( 1+ c*re ibe- 2a—jBN ) 
[1 —|cP re dal} ( ] +crej¥e-2(a+jB)!) 





where (*) denotes the complex conjugate. 

Although this expression is a complicated one, a 
great deal can be learned about the problem by in- 
spection of this result. First, if the short is ideal, 
r=1, JW=7, a=0, and the expression becomes: 


bs - bac” (a— bc) e728" 1 — c* eb!) 
bs 1—|el? (1 — |e|?) (1 —ce-258!) 





(10) 


Here it is noted that 63/6, is the sum of two terms. 
The first is a constant (as / varies), while the second is 
of constant magnitude but variable phase. (Note that 
the last factors in the numerator and denominator of 
the second term are conjugates of each other.) This 
explicitly demonstrates the previously stated result: 
the locus of b3/bs is a circle of radius 








R= =e (11) 
og 
and distance to the center 
iil 
R.= Tce : (12) 





Next, it is desirable to return to the more general 
relationship contained in (9) and to consider the 
functional dependence of |b3/bs| upon /. It is this 
relationship which is plotted in figure 6. 

















SHORT POSITION——————o> 


FIGURE 6. “Typical” response in an actual measurement. 





In the discussion which follows, it will prove useful 
to consider not only the functional relationship of 
(9), but also the family of functions which is obtained 
from (9) tor different values of Ww. Conceptually, one 
assumes the measurement is repeated a number of 
times, where the moving “shorts”’ have the same value 
for r, but differing values of w. 

For a given /, the value of |b3/bs| will depend upon 
Ww, and for some particular choice of W will have the 
maximum value: 


bs _ ja— bel re“ + 1b—ac*re = | 





(13) 








bs max i~ joyfrte** 


For some other choice of W, |b3/bs| will have the 
minimum.?: 


bs _ ja—bel|re-?! — |b — ac*re~™ | 


a 1 — |c|2r2e-4a! (14) 











bs\min 
For each value of /, there will be some member of the 
family (value of w) for which these maximum and 
minimum values are realized, but these limits are never 
exceeded. Equations (13), (14) are evidently the 
envelope of the family. If this envelope is denoted 
by E, then: 


_ ja—beire“™ + |b— ac* re tat} 


E= (15) 


| ins] ic|?re -4al 





Because of the oscillatory behavior of |b3/b4|, the 
envelope may be easily inferred from a single member 
of the family provided that several cycles, or more, 
have been recorded. From (15), the extrapolation of 
this envelope to the plane where /=0 yields: 





la— be|r+ |b—ac*r*| 


E 








i=o= (16) 


1—|c|?r 


The use of this result to obtain approximate values for 
the efficiency will be evaluated in the following section. 
It is next instructive to consider the general behavior of 
E as | varies. 

Returning to (15), the first term in the numerator 
tends to dominate the entire expression, and provides 
a simple exponential decay. The second term in the 
numerator determines the amplitude associated with 
the oscillatory behavior. If c vanishes, this amplitude 
is constant, while if 6 vanishes, the amplitude decays 
exponentially. If both these conditions obtain the oscil- 
latory behavior is absent. The most distinctive situa- 
tion occurs when 6 and ac* are nominally equal. In this 
case there will be some value of / (which may or may 
not be physically realizable) for which this second term 
will be a minimum (possibly zero). 

Finally, the second term in the denominator is usu- 


* This assumes that the first term exceeds the second which, for the applications en- 
visioned, is usually the case. 
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ally small ® with respect to unity. As a whole the denom- 
inator increases as / increases, and the value of the 
complete expression decreases somewhat more rapidly 
than would otherwise be the case. As a rule, this effect 
is rather small. The most general form of the envelope 
is a distorted hyperbola as shown in figure 7. 


| 











bs 
bal | 
| 
| 
SHORT POSITION —= 
FIGURE 7. General form of system response. 


5. Practical Considerations 


As already noted, the measurement technique calls 
for recording |b3/bs| as a function of / (short displace- 
ment). The envelope of this curve is then projected 
to the plane where /=0, and Eyjax, Emin Substituted for 
bs/bs| in (3), (4). The resulting approximations to R, R- 
are denoted by primes and given by: 





»» . @--bei(i—e) 
i 1—|c|*(1—e)* 





(17) 


Ri= ont Sa a 
l—jc#(1—e)* 





(18) 


where r= 1—e, and € is usually a small quantity. 

The evaluation of the error, e, due to nonzero e, 
in an efficiency determination using (5), begins with 
the definition : 


——_———a (19) 


Substitution of (11), (17) into this result leads to 


e = 2e(|c2|? — |e, |?) (20) 
where the subscripts 1, 2, differentiate between the 
ereflectometer parameters which obtain at ports 1 
and 2 respectively, and only the first order terms have 
been retained. 

If |ei|, |ce|, € were known, their values could be 
substituted into (20) and e would become a “‘correc- 
tion factor” rather than error. As a rule, however, the 








‘ht is not within the scope of this paper to consider the many “pathological” situations 
which could be invented, e.g., couplers without directivity or reversed in direction, extreme 
departures from impedance match conditions, etc. Because the mathematical formulation 
is general enough to include these cases, these additional qualifying statements are required. 


potential accuracy improvement does not warrant 
the extra effort this requires. As a general guideline, 
in the configuration of figure 3, |c| ~ |['y| where IT, 
is a reflection coefficient which is “representative” 
of the two-port (e.g., adaptor) and the load for which 
the efficiency is measured. In the configuration of 
figure 4, |c| is determined by the adjustment of 7,; if 
this is such as to minimize the amplitude of the oscil- 
latory behavior of |b3/b4| in response to the short mo- 
tion, the foregoing guideline (|jc| ~|I,|) is again 
satisfied.® 

For an adaptor of VSWR of 1.5, a typical value of 
\c| is 0.2, while € will ordinarily be less than 0.01. The 
maximum expected error, due to the nonideal short is 
thus 0.1 percent over a wide range of practical oper- 
ating conditions, and the typical error is probably 
0.01 percent or less. 

To complete this discussion, there will also be an 
error in the determination of R,;, Rez due to the non- 
ideal short. Ordinarily the error contribution from 
this source will be one or more magnitudes smaller 
than that already described and thus negligible. 

Thus far, the discussion has implicitly assumed 
that the same short is used at both terminals | and 2. 
In some applications, the evaluation of a waveguide- 
coax adaptor for example, this is obviously impossible. 
In this case, the generalization of (20) is 





ec (€2 —€,) + 2(€2|\co|* —€1\c1|?). (21) 
Here €;, €2 are associated with the shorts used at ports 
1 and 2 respectively. 

It will be immediately recognized that the difference 
between €; and €2 is now a potential source of sub- 
stantial error. Fortunately, however, this error can be 
eliminated. The procedure is to make the measurement 
by each of the two described procedures (using (5) 
and (6)).?7 Assuming that the same pair of shorts are 
used for each of these measurements, it is easily 
shown that the geometric mean of the two efficiency 
determinations gives the desired result. Conversely, 
the square root of the ratio of these measurements 
gives the ratio between the reflection coefficient 
magnitudes of the respective plungers. 

Returning to the expression for efficiency, it is 
convenient to write (5) in the form: 


ewe. 8 
_ 


22 
R (22) 


This shows that the vertical displacement between 
R,z and R;, (fig. 6) is a direct measure of the difference 
between 72; and unity. Typically, this difference is no 
more than a few percent, and the performance require- 


6 Although this guideline should prove adequate for most practical purposes, an “exact” 
evaluation of ¢;, cz is not difficult to obtain. In particular, the generator in figures 3 or 4 is 
replaced by a variable termination which is adjusted such that 5, vanishes when the system 
is excited via port | (or 2). The reflection coefficient observed at port | (or 2) is now —¢, 
(or — C2). 

7 This assumes that in the first of the described methods 7, has been adjusted such that 
“Ty” equals the load reflection coefficient for which yz is desired. Moreover, it will be 
recognized that, during the second measurement, the two ports of the adaptor are reversed, 
in relation to the g-reflectometer, as compared with their positions during the first 
measurement. 
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ments on the associated measurement system are 
primarily for resolution and stability rather than 
“absolute” accuracy. For example, a 10% error in 
either Re or in (Re—R,) will cause an error of only a 
few tenths of a percent in m1. Conversely, a system 
resolution of a part in 10° is required if (R2—R,) 
is to be observed with a resolution of one percent, etc. 


The role of tuner 7, has been referred to a number of 


times. Although the operation is, in theory, inde- 
pendent of this adjustment, certain second order error 
considerations have been noted. In addition, however, 
is the practical requirement of matching the excur- 
sions of |b3/b4| to the dynamic range limitations of the 
x—y recorder. If these excursions are permitted to 
be large relative to Ro — R,, the resolution, with which 
R.—R, is measured, suffers. With a little experience 
this adjustment is not a difficult task. 

Finally, it is of interest to note that the recording pro- 
vides not only an indication of the joint or other losses, 
but also an indication of the loss of the transmission line 
in which the short moves. In many cases, this provides 
a useful confirmation of the system sensitivity and cali- 
bration. Indeed, it is possible to measure the loss in 
terms of the waveguide attenuation/wavelength with- 
out otherwise calibrating the system. 


6. Applications 


This technique represents a powerful tool for certain 
measurement problems. An immediate example, al- 
ready alluded to, is the evaluation of connector loss. In 
order to gain some insight into the behavior of the Type 
N connector, a coaxial line was devised with a contin- 
uous outer conductor and a slip joint on the inner con- 
ductor such as is found in the Type N connector. This 
joint showed a loss of 0.006 dB at 9 GHz. In another test 
piece, the center conductor was continuous, while the 
outer conductor included a replica of the joint found 
in the Type N. This exhibited a loss of 0.05 dB.2 The 
measured loss of a complete “Type N” joint, in which 
the slots had been deleted from the outer sleeve, was 
0.01 dB. 

This technique has also found application in the 
evaluation of certain parameters associated with a 
cryogenic noise source developed for the Comsat Cor- 
poration. In particular the corrected value of the tem- 
perature at the output port required both the loss per 
unit length in the output waveguide, and the loss in a 
waveguide window. Both of these quantities were easily 
measured via these methods. 

Another use of these methods is in conjunction with 
the large antenna gain calibration project at the Jet Pro- 
pulsion Laboratory [6]. Here a loss correction factor is 
obtained for the 30-40 foot waveguide run between the 
gain standard horn and the reference plane. The mis- 
match corrections are also explicitly accounted for with 
this technique. 


* Although the purpose of these measurements was ostensibly that of evaluating the con- 
nector design, the real motivation was that of demonstrating the flexibility of the measure- 
ment technique. The above result is rather surprising, and may not be actually representative 
of the general design. In any case, it appears desirable to further pursue this subject 





7. Summary 


The usefulness of the sliding short methods of “‘atten- 
uation”’ measurement has been substantially enhanced 
by application of the first order correction theory, for 
nonideal short behavior, developed by Almassy [5], 
and the subsequent extension outlined in this paper. 

The microwave art is rapidly approaching the place, 
if indeed it is not already there, where connector imper- 
fections represent a major barrier to further advances 
in the accuracy of microwave measurements. Although 
a great deal of recent effort has been expended in the 
direction of improved VSWR specifications, the more 
important parameter, for many applications at least, is 
the dissipation characteristic. This method should 
prove a particularly valuable tool in assessing the mag- 
nitude of this problem. 

The procedure may be implemented in a variety of 
ways, the major requirement is for a high degree of sta- 
bility in the measuring system. One such system has 
been described in some detail by Almassy ® [5]. At the 
Jet Propulsion Laboratory, a commercial version of the 
“precision insertion loss test set’ developed at JPL 
[7] has been adapted for this use. 

The existence of this first order correction theory 
also makes it possible to envision the extension of slid- 
ing short methods to much lower frequencies than has 
previously been feasible. In particular, because the 
theory explicitly takes account of line loss, it may prove 
useful to construct the inner and outer conductors in 
the form of a helix. Another possibility is that of im- 
mersing the line in a liquid dielectric. 

Finally, the tuning adjustments, called out in this 
procedure, can be eliminated if a phase detection capa- 
bility, such as is found in network analyzers, is assumed. 
Much of the detail for doing this has been worked out, 
but remains to be implemented. 


8. Appendix 


The purpose of this appendix is to derive (6). 

Referring to figure 4, the relationship between b3/bs 
and the reflection coefficient, I’, presented at terminal 
2 by the two-port and its load, I), may be written,” 


b, _AT,+B . 
b, CI,+D (23) 


where A, B, C, D are parameters of the g-reflectometer. 
In a similar way, the relationship between I’, and I’; 
is given by, 


._ali+b ; 
ors. cl; 4 l (24) 


* The “compensation” included in Almassy’s system has not been found necessary in 
these measurements 

© She numerator and denominator could obviously be divided by D to obtain better agree- 
ment with the functional form of (7). The existing form, however, provides better continuity 
with earlier work in this area 
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where for the moment, I; is assumed to be arbitrary, 
and a, b, c are parameters of the two-port whose efh- 
ciency is required. (Note that this represents a change 
in terminology from (7).) 

The required efficiency is given by [2], 


|la—be|(1—|Ti|?) 
1+ cl )2?— lal) + bP 





M2(11) (25) 


By use of (11), (12), 
shown that 


and a little algebra, it can be 


Re) _ lAl— [BP 
ni-(Heh we 
R:) ) \|AD—BC| ats 
To obtain the counterpart expression involving 


R,, Rey it is convenient to substitute (24) into (23) 
which (for an arbitrary |")) yields, 


bs _ (Aa+Be)Vi+ (Ab+B) | 








bs (Ca+Dc)l;+ (Cb+D) 2%) 
This is in the same form as (23) so, 
Rs (1- (Fe) ) 
|Aa+ Bc|?—|Ab+ B|? , 
= (da+Be)(Cb+D)—(Ab+B)(CatDea| 
By hypothesis, however, 6; vanishes when I; 


represents the load reflection for which the efficiency 
is desired. Therefore, 


(Ab+B) =— (Aa+Bc)Ii, (29) 


from which, 
= al';+b 


r cl,+1 





(30) 





Substitution of this result into (28), (26) leads to, 


w=) 


|A| - |a—be| - (1— ||?) 











~ fitch: (Clari+6)+DU+ery)* = 

Re 2 

Re (1-()) 
___|Al- ({l+el|?—|ali+ b)?) ‘i 


~ |l+eP|- |C(ali+6)+DO+cr))| 


Finally, comparison of the ratio of (31) to (32) with 
(25) yields (6). 
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Accurate Microwave High Power Measurements Using a 
Cascaded Coupler Method 


Kenneth E. Bramall* 


Institute for Basic Standards, National Bureau of Standards, Boulder, Colorado 80302 


The use of directional couplers to extend the range of low-level power meters is a well established 
technique. Calibration of bolometer-coupler units up to 20 dB is routine at NBS and other laboratories 
and allows power measurement to the 1 W level. For higher coupling ratios, however, the problems in 
calibration using conventional techniques becomes more difficult. This paper describes a cascade 
method which allows calibration of low power meter-coupler combinations for measurement up to the 


megawatt range. 


The technique is based upon system linearity and the uncertainty limits are shown to be 1 percent 
to 1.5 percent using NBS low power standards. By using the NBS calibration services and commercially 
available equipment, the above power levels can be measured with uncertainty limits less than 2.7 
percent. The method can also be adapted to measurement of large coupling ratios (20-80 dB). 

Results are given for power (200 watts in WR-90 waveguide) and coupling ratio (40 dB coaxial 


coupler) measurements. 


Key words: Accurate: calibration; directional coupler; high power; measurement: microwave. 


1. Introduction 


High level (above 1 W) microwave CW power meas- 
urements are commonly referenced to low-power 
standards by using a low-level power meter in com- 
bination with a power-reducing device such as a direc- 
tional coupler or attenuator. Although better known and 
more economical than other methods for measuring 
high power, use of this method has been limited by 
uncertainties in relating the high power to the power 
indicated by the low-level power meter. An improved 
technique for determining this relationship is the main 
topic of this paper. The use of the directional coupler 
is emphasized, but the basic ideas presented will also 
apply to other attenuating or coupling devices. 


2. Cascaded Feed-Through Power Meter 


Calibration Technique 


The basic principles of the new calibration tech- 
nique are described by discussing the calibration of 
a 40 dB coupler in combination with a low-level side- 
arm power meter having a dynamic range of 20 dB 
(0.1 to 10 mW). 

It has been shown that the best accuracy is obtained 
by calibrating the sidearm power meter and the coupler 
as a unit [1].! This requires a measurement of the 
coupler output power while a reading is obtained on 
its sidearm power meter. If a standard power meter 


*Electromagnetics Division, Power-Current-Voltage Standards, National Bureau of 
Standards, Boulder, Colo. 80302. 
Figures in brackets indicate the literature references at the end of this paper. 


is available for measuring power in the range of 1 to 
100 W (the power needed to obtain a reading within 
the range of the sidearm power meter) the calibration 
can be readily achieved. For the sake of explaining the 
technique, however, the calibration is to be achieved 
using a 10 mW standard power meter. 

The calibration procedure employs an auxiliary 20 
dB coupler and sidearm power meter with a dynamic 
range of 20 dB (0.1 to 10 mW) to transfer the calibra- 
tion from the 10 mW standard power meter to the 40 
dB coupler unit. Before starting the calibration, the 20 
dB coupler is connected to the 40 dB coupler as shown 
in figure 1. At first glance, it appears that two coupler 
type, feed-through powermeters are in series. However, 
an alternate interpretation is also possible. In partic- 
ular, the 40 dB coupler can be regarded as an extension 
of the main arm of the 20 dB coupler so that terminal 
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FIGURE 1. Basic cascaded coupler system. 





surface 2 serves a dual role; it is the output port of both 
the 20 and 40 dB couplers. 

The 20 dB coupler unit (containing the 40 dB coupler 
unit) is calibrated by obtaining simultaneous readings 
on its sidearm power meter (PM1) and on the 10 mW 
standard power meter connected to terminal surface 2. 
If the standard power meter is now removed and re- 
placed by a high power load (and for the present it will 
be assumed that their impedances are equal) the 20 dB 
coupler unit will permit measurements of the load 
power in the range of 10 mW to 1 W. The only require- 
ment is that the component parameters be linear, i.e., 
independent of power level. 

As the 1 W level is approached, the indicator on the 
20 dB coupler reaches the upper limit of its dynamic 
range (10 mW) and the indicator on the 40 dB coupler 
comes into the lower region of its dynamic range (0.1 
mW). Since the load power can be accurately predicted 
from the power meter on the sidearm of the 20 dB 
coupler, the 40 dB coupler unit is calibrated by obtain- 
ing simultaneous readings on the two sidearm power 
meters. The calibration is thus transferred from the 
20 to the 40 dB coupler unit without introducing uncer- 
tainties associated with mismatch, with flange losses, 
or with losses occurring within the couplers. Following 
this calibration, the input power can be further in- 
creased another 20 dB and load powers in the range of 
1 to 100 W determined by the 40 dB coupler-sidearm 
power meter unit. Again, the components must be 
linear. 

As a practical matter, the power meter on the side- 
arm of the 2U dB coupler may be damaged by excessive 
power when the load power is increased beyond | W. 
This power meter overload problem is simply avoided 
by removing the 20 dB coupler after the 40 dB coupler 
unit is calibrated. The calibration of the 40 dB coupler 
unit is not affected by doing this because a directional 
coupler, like other three-port devices, has the property 
that the power split between the main and sidearms is in- 
dependent of the source characteristics (see appendix). 

The method described in the previous paragraphs 
can obviously be extended to systems employing any 
number of cascaded couplers. For example, figure 2 
shows four coupler units where the calibration pro- 
cedure is as follows: 


1. Adjust power input so both PM] and PMS are 
within their dynamic ranges. 

. Obtain simultaneous readings of PM] and PMS. 

. Replace PMS with the high power load. 

Increase power until PM2 is in its operating 

range (but upper limit of PM1 is not exceeded). 

. Obtain simultaneous readings of PM] and PM2. 

. Remove PM] and the 20 dB coupler from the 
system. 

7. Increase power until PM3 is in its operating 
range (but upper limit of PM2 is not exceeded).? 


wh 


nu 


2 Since coupler 1 was removed from the system, PM] may be transferred to coupler 3. 
Likewise PM2 may be transferred to coupler 4. Thus, only two power meters are required 
for any number of cascaded couplers. Although errors may be introduced by this procedure, 
they are negligible for couplers having large coupling coefhcients. For couplers having small 
coupling coefhcients, the error can be made negligibly small by using load isolators on their 
sidearms or by initially terminating the sidearms with loads having refiection coefhicients 
approximately equal to the sidearm power meters. 
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FIGURE 2. Cascaded coupler system used to calibrate an 80 dB 
coupler-power meter combination. 


8. Obtain simultaneous readings of PM2 and PM3. 

9. Remove PM2 and the 40 dB coupler from the 
system. 

10. Increase power until PM4 is in its operating 
range (but upper limit of PM3 is not exceeded). 

11. Obtain simultaneous readings of PM3 and PM4. 

12. Remove PM3 and the 60 dB coupler from the 
system. 

13. Measure power in the range 104 to 10° W. 


By following the steps above, the 80 dB coupler and 
its sidearm power meter are calibrated as a unit and 
the load power is given by the expression 


_ Pi Ps Ps Ps 
Pu= BP, Ps Py (1) 


where P,, is the reading of the nth side arm power 
meter when a calibration was effected and P}, is the 
reading of this power meter when it measured the load 
power. The reading of the standard power meter is 
P,. This expression for the load power, P,, is exact 
when the standard power meter and the load im- 
pedances are equal. The effects of dissimilar imped- 
ances will be considered next. 

When the standard meter and the load impedances 
are not equal, mismatch is introduced when the former 
is exchanged for the latter. This is the only impedance 
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mismatch introduced, however, because the same 
load is used for the remainder of the calibration 
procedure. The impedance mismatch introduced with 
the cascaded system is basically the same as for a 
single coupler unit and can be handled using existing 
techniques. 

Engen [1] has shown that the load power determined 
by a single calibrated coupler unit is given by the 
expression 


tp L—IPi) |1—Voels|? 
P,= —! P, . be ; 
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(2) 








where P, is the reading of the side arm power meter 
when the unit was calibrated with a standard power 
meter whose reading was P,. The reading P| of the side 
arm power meter determines the power delivered to 
the load after it replaces the standard meter. The 
reflection coefficients of the standard meter and the 
load are given, respectively, by I’; and I’,. The factor 
Ice is defined by Engen [2] as the equivalent gen- 
erator reflection coefficient which the coupler provides, 
and its measurement is discussed in the referenced 
paper. This coefficient is given in terms of the coupler’s 
scattering coefficients as 


Pir =S»—-—* (3) 


when the reference planes are chosen so that the 
source is connected to port 1, the output of the coupler 
main arm is port 2 and the sidearm is port 3. 

Similarly, an equivalent generator reflection coefh- 
cient ['¢¢ can be defined for the first feed-through 
power meter in the cascaded coupler system with 
reference planes as indicated in figure 2. By doing 
this and comparing the expression for the load power 
given in eq (1) with eq (2), eq (1) can be modified to 
include the effects of impedance mismatch. The results 
are: 
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Only the [ce of the first feed-through power meter 
appears in eq (4), although an equivalent generator 
reflection coefficient can also be defined for each of 
the other feed-through power meters in the cascaded 
system. Unlike the first feed-through power meter, 
the others measure the power delivered to the same 
load that terminates their output when they are cali- 
brated. As a result the power measurement is not 
affected by their equivalent generator reflection coefh- 
cients. They need not be known unless, as discussed 
below, the load is exchanged for another after the 
calibration is complete. The Ice of the first feed- 
through power meter must be known, however, in 
order to evaluate the mismatch error in eq (4). This 
requires that I~ be measured before the calibration 
begins and the system is disassembled. 

In addition to giving a method for measuring I¢¢ in 
the above reference, Engen showed that it is possible 


to reduce the value of Ics for a single coupler by using 
a tuner in the output arm. A tuner can also be incorpo- 
rated into the cascaded system and used to reduce the 
value of Ie that exists when the load replaces the 
standard meter. The tuner can be added anywhere be- 
tween the first coupler and the output port. There is an 
advantage, however, in locating it between the first and 
second couplers. After the calibration is transferred to 
the second coupler, the tuner becomes part of the 
source for this coupler. Thus the tuner can be removed 
with the first coupler, eliminating the need for a tuner 
that will operate at high power levels. 

An alternative method for measuring and reducing 
[Ge is provided by employing the first coupler as part 
of a reflectometer [3, 4]. When this reflectometer is 
tuned with respect to the output port, gz is made small 
and a method is provided for measuring the reflection 
coefficient magnitudes of the standard meter and the 
load. Again, as with the single tuner, the added reflec- 
tometer components can be removed with the first 
coupler. More details about this method are given 
later with the experimental results. 

The high power load can be exchanged for another 
after the calibration is complete. However, another mis- 
match error will be introduced. Specifically, if after the 
last coupler has been calibrated using a particular load, 
another load with a reflection coefficient I), is substi- 
tuted, the expression for the power delivered to this 
second load is 
2 
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where I, is the equivalent generator reflection coefh- 
cient associated with the last coupler, and the other 
terms are as previously defined. Evaluation of the addi- 
tional mismatch error in eq (5) requires additional 
measurements to determine I%;¢ and Iie. 


3. Application to Attenuation Measurements 


The calibration technique as described provides a 
method for accurately determining the power delivered 
to a particular high power load. The technique does 
this by using several power meters with limited 
dynamic ranges to produce the same results as would 
be obtained using one power meter with a large 
dynamic range. This characteristic also makes the 
technique useful for other applications such as 
evaluating the linearity of detectors with large dynamic 
ranges or for measuring attenuation using the power 
ratio method of measurement. 

The power ratio method for measuring attenuation 
requires measurements of the input and output 
powers of a two-port device. This method has been 
shown to be capable of high accuracy when using 
bolometric power meters for attenuation measurements 
up to 20 or 30 dB [5]. A wider range of attenuation can 
now be measured using a cascaded coupler system, 
since it permits a wider range of input power levels 
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to be measured. For best accuracy, the output power 
can be measured with the same standard power meter 
used to calibrate the cascaded coupler system. The 
maximum attenuation that can be measured is limited, 
of course, by the component’s ability to withstand the 
high power. This limitation, for several measurements, 
is not as severe as it may seen. For example, by using 
power meters with dynamic ranges of 0.1 to 10 mW, 
40 dB of attenuation can be measured with a source of 
1 W, 50 dB with a source of 10 W, and 60 dB with a 
source of 100 W. Several components including direc- 
tional couplers can tolerate these (and higher) power 
levels. 

The measurement of a coupler coupling coefficient is 
an attenuation measurement between the input and 
the side arm when the output is terminated in a non- 
reflecting load. An 80 dB coupler, for example, requires 
an input power level of 10 kw in order to obtain a 
reading of 0.1 milliwatt at the side arm. A measurement 
of this input power is provided by using a cascaded 
system of 20, 40, and 60 dB coupler units. Using this 
system and computing the coupling coefficient yields 


Pi P: Ps Pp |1—Voels|?___ 1 
P,P. P, P| 1—Falc 11 —SssIF'; |? 


where Sm», are the scattering coefficients of the coupler 
and I¢ is the reflection coefficient at the input to the 
coupler with the nonreflecting high power load and 
standard power meter connected. The expression for 
this reflection coefficient is 





K=10 Logio (6) 





: Si3Ssil’s 

[ c=Su + ] we (7) 
Note that since the ratio of P;/P is taken, where Ps. 
is the side arm power reading obtained with the stand- 
ard meter, a calibrated standard power meter is not 
required for these measurements. Consequently, 
this meter can be incorporated with a tuner to obtain 
[,~0 and reduce the mismatch error introduced. 
The advantage of tuning [g¢e¢-- 0 at the low power 
level is apparent. 


4. Alternate Calibration Method 


The system shown in figure 3, is an alternate 
calibration method which uses attenuators between the 
coupler sidearms and the power meters. Compared 
with the purely cascaded system, it is physically 
shorter and the power ratio calibration steps can have 
any value if the attenuators (A; and Aes) are con- 
tinuously variable. The utility of this method is that it 
allows couplers having relatively low coupling ratios to 
be transformed into ones having the equivalent of 
high coupling ratios. For example, a 40 dB coupler can 
be transformed to the equivalent of a 60 or 80 dB 
coupler. 

The calibration procedure for this system is similar 
to that for the purely cascaded system, the difference 
being that the calibration is transferred back and 
forth between power meters | and 2 as the power level 








sR 


FIGURE 3. Variation of the cascaded coupler calibration technique. 
The attenuators are alternately used to reduce the coupling to one power meter while 
the other retains the calibration. 


is increased. The attenuators are alternately used to 
reduce the coupling to one power meter while the 
other retains the calibration. 

The first coupler unit can be considered as part of 
the source for the second unit and the setting of 
attenuator A, can be changed without introducing an 
error when the calibration is referenced to power 
meter 2. When the calibration is referred to power 
meter 1, however. a change in the setting of attenuator 
A» can introduce an error, because the amount of 
power reflected back into the system may change 
which in turn can change the relation between the 
first sidearm power meter reading and the load power. 
This error can be made negligibly small by using an 
isolator as shown in figure 3. The other errors associ- 
ated with this calibration scheme are the same as 
for the strictly cascaded coupler system. 

A tuner for reducing [ge can be inserted between 
the couplers. It can be removed after the calibration 
is first transferred to power meter 2 and an error will 
not be introduced when the calibration is transferred 
back to power meter 1. Similarly, a reflectometer can 
be used to reduce and evaluate the mismatch error 
introduced at arm 2. 


5. Experimental Systems and Results 


A WR-90 waveguide coupler system was assembled 
in the laboratory to obtain experimental data. The 
system incorporated a reflectometer and provided the 
capability of obtaining two separate calibrations of the 
last low power meter-coupler unit. This system is 
shown in figure 4. One calibration for the 40 dB coup- 
ler unit was obtained by using couplers 1 through 4 
to provide 10 dB power ratio calibration steps. The 
other calibration was obtained by using couplers 2 
and 4 to provide power ratios of 20 dB. The two cali- 
brations disagreed by no more than 0.1] percent which 
indicates the precision of the method. 

No unusual problems were encountered or unpre- 
dicted results obtained in performing the double 
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FIGURE 4. A coupler system that incorporates a reflectometer and provides the 
capability of obtaining two separate calibrations of the 40 dB coupler. 


calibration. Intercomparison power measurements 
were made at the 70 to 200 W level between the 40 
dB coupler unit and a stirred-water calorimeter [6]. 
The disagreement was around | percent, well within 
the combined error limit for the two meters. 

In the tuning procedure for this system, tuner A 
was first adjusted so that the coupler sampling the 
reflected power (coupler 5) would have near infinite 
directivity. Tuner B was then adjusted to obtain 
lce =~ 0 for the system using power ratios of 20 dB. 
Finally, tuner C was adjusted to obtain [ge = 0 for 
the system using power ratios of 10 dB. All of the tuning 
was done with respect to the output of the system by 
using well known tuning procedures [3, 4]. The tuning 
operation was performed by sliding high and low cali- 
brated reflecting loads within a waveguide section 
attached as part of the coupler system. The low re- 
flection load was used as the reflection coefficient 
standard. By using this approach, the reflectometer 
measured the combined reflection coefficient of the 
load and flange at the output of the system. 

The reflectometer, the first two couplers, the three 
tuners, and the low-power low-pass filter in the main 
line were removed when the calibration was trans- 
ferred to couplers 3 and 4. 

A second system using coaxial components was 
assembled to measure the output to sidearm power 
ratios of a set of 5 coaxial couplers having nominal 
coupling values of 12, 20, 23, 30, and 40 dB. These 
are slot type, quarter wave precision couplers as de- 
scribed by Hudson [7] and have 50 dB directivity and 
maximum VSWR of 1.02 at all ports. Their center 
frequency is 1.45 GHz and measurements were made 
at frequencies in the range | to 2 GHz. For the lower 
value couplers (up to 20 dB), the output-to-sidearm 
power ratio was measured directly using standard 
bolometric power meters (VSWR < 1.03). The cascade 
system was used to measure the higher value couplers, 
and because VSWR of the ports is low, tuners were 
not required to reduce the value of Te. The limit of 
uncertainty in the measurement of the 40 dB coupler 
was computed to be 0.8 percent (0.035 dB) which 


includes the effects of impedance mismatch, repeat- 
ability, etc. The 40 dB coupler was also measured 
using the NBS coaxial calorimeter [8] and a standard 
bolometric power meter. The two measured power ratio 
values differed less than | percent. 

6. Errors 

The limit of uncertainty in a high power measure- 
ment, using the cascade system, depends primarily 
upon the uncertainties introduced by the reflection 
coefficients and by the low level standard and side 
arm power meters. Since bolometric power meters 
conveniently yield accurate measurements in the 0.1 
to 10 mW range, they are used in the following dis- 
cussion to demonstrate the accuracy obtainable. 

There are at least two methods available for obtain- 
ing less than 0.1 percent for the uncertainty introduced 
by the reflection coefficients. In one, the uncertainty 
is eliminated, in principle, by adjusting the reflection 
coefficients of the standard meter and the high power 
load to be equal. In the other method, the uncertainty 
is reduced by using a tuner to obtain a small, near 
zero, equivalent generator reflection coefhicient and 
then correcting the measured power using the meas- 
ured reflection coeficient magnitudes of the standard 
meter and the load. 

The limit of uncertainty associated with the meas- 
urement of P, results primarily from the uncertainty 
in the bolometer mount effective efficiency and from 
the uncertainty in the measurement of the bolometri- 
cally substituted de power. The National Bureau of 
Standards presently provides routine effective efh- 
ciency calibrations with limit of uncertainty of 1 percent 
for waveguide bolometer units [9]. The uncertainty 
can be as low as 0.2 percent, however, for a WR-90 
waveguide mount calibrated against the NBS micro- 
calorimeter [10]. 

Self-balancing bolometer bridges have been devel- 
oped that enable a measurement of the substituted 
dc power to be made with an uncertainty of 0.09 per- 
cent +0.1 «W [11]. This corresponds to uncertainties 
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of 0.09, 0.1, and 0.19 percent in measurements of 10, 
1 and 0.1 mW, respectively. 

In a power ratio measurement, some of the syste- 
matic errors, such as those resulting from the bolom- 
eter mount reflection coefficient and its effective 
efficiency, cancel. The only uncertainties that need 
to be evaluated are those associated with the measure- 
ments of the bolometrically substituted de powers. 
By using the bolometer bridge mentioned above, power 





ratios of 10 and 100 can be measured with uncertainties 
of 0.19 and 0.28 percent, respectively. 

The uncertainty in a high power measurement is 
the sum of the uncertainties introduced by the standard 
power measurement, by the mismatch error, and by 
the uncertainties introduced in the power ratio meas- 
urements. Table 1 shows an example breakdown of 
these uncertainties. 

The power ratio uncertainties in table 1 can be re- 
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wae” TOTAL UNCERTAINTY IN A HIGH POWER MEASUREMENT (3) 
100W 1.95 1.75 1.15 0.95 
10kW 2.33 2.03 1.53 1.23 
MW 2-71 2.31 1.91 1.5] 














TABLE 1. Examples showing the uncertainty limit obtainable in a 


high power measurement by using the cascaded coupler calibration 


technique and equipment that is available commercially. 
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duced by the use of the same bolometer bridge and 
accessory instrumentation in measuring the numerator 
and denominator powers. The amount of reduction 
will, in part, depend upon the systematic error common 
to both measurements. 

Other sources of error that need to be considered in 
practice include rf power instability and the effects 
of bolometer mount temperature changes. Because 
the calibration is performed by making simultaneous 
power measurements, only short term stabilities are 
important, and the exacting power and temperature 
stability requirements are reduced. 

The output from some high-power sources may 
contain large harmonic components. Since the coupling 
for the harmonic power is generally not the same as 
for the fundamental power, this can be a source of 
error. This error is easily reduced by incorporating 
low-power, low-pass filters with the power meters, as 
shown, for example, in figure 4. 

In general, the calibration technique is based upon 
the assumption that the system is linear. The param- 
eters of all the components used in the system must 
be independent of the rf power level and any other 
variable, such as environmental temperature, that may 
affect the linearity. The linearity of all the components 
except the directional couplers can be easily measured, 
since they operate at low power levels. The linearity of 
a coupler unit, however, must be measured at high 
power levels. Since this measurement requires a high 
power meter, another coupler unit can be used. In 
this way the power linearity of different types of 
couplers can be compared by comparing the side arm 
power readings as the power level is changed. With 
several comparisons any nonlinear components can 
be detected. 


7. Summary 


Accurate measurements of microwave power at 
high CW levels can be accomplished using the cas- 
caded feed-through power meter calibration technique. 
The accuracy is based upon existing low-level stand- 
ards, system linearity, and the ability to accurately 
measure power ratios. The technique can be used at 
any frequency in the rf spectrum where directional 
couplers (or other suitable 3-port devices), low-level 
power meters, and low-level power standards are 
available. (See fig. 5.) 

Alternative systems, which employ directional 
couplers and attenuators, provide an additional 
flexibility in the calibration of coupler-low-power- 
meter units. The use of feed-through devices other 
than directional couplers may be. advantageous, 
especially at the lower frequencies where the wave- 
guide components are large. Since absolute power 
measurements are not required for measuring power 
ratios, detectors other than power meters, such as a 
precision variable attenuator and detector combina- 
tion or an IF detector, are possible alternatives. Both 
of these detectors, however, may be subject to addi- 
tional errors cuased by impedance changes. 
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Large values of attenuation or coupling ratio can 
be accurately measured using the cascade technique, 
provided the components can withstand the high power 
required. 
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8. Appendix 


The calibration technique uses one calibrated 
feed-through power meter to calibrate another, and 
since a feed-through power meter has an input and 
two outputs, the theory for 3-port junctions is appli- 
cable. The ratio of the powers at ports 2 and 3 of any 
linear 3-port junction, with a source connected to 
port 1 as shown in figure 5, is given by [12] 
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where Syn are the scattering coefficients of the 
junction, and I, and [; are the reflection coefficients 
of the loads at ports 2 and 3. 

The expression shows that the ratio of P2 to Ps is 
independent of the source reflection coefficient, I. 
For a linear system, the reflection and scattering co- 
efficients are not a function of the power level and the 
ratio is independent of the source amplitude. In 
general, the ratio, or the calibration of a feed-through 
power meter, is independent of the source character- 
istics, and an error is not introduced by changing the 
source at port 1. 
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FIGURE 5. A generalized 3-port junction. 
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type samples, and have failed to provide a well-characterized widely 
distributed and easily available standard. The NBS Office of Stand- 
ard Reference Materials is in the process of acquiring, analyzing and 
certifying a series of biological standards. These standards, when 
issued, will include six botanical standards, a blood standard and a 
tissue standard. 


Bell, F. C., Wollin, H. F., Editors, Report of the 55th National 
Conference on Weights and Measures 1970, Nat. Bur. Stand. 
(U.S.), Spec. Publ. 342, 315 pages (April 1971) $1.75, SD Catalog 
No. €C13.6/3:342. 


Key words: Conference, weights and measures—laws; weights and 
measures — regulations; weights and measures—technical require- 
ments; weights and measures—history: weights and measures. 


This is a report of the proceedings (edited) of the Fifty-fifth National 
Conference on Weights and Measures, sponsored by the National 
Bureau of Standards, held in Salt Lake City, Utah, July 13-17, 
1970, and attended by state, county, and city weights and measures 
officials and representatives of the Federal Government, business, 
industry, railroads, and associations. 


Bennett, H. S., Absorbing centers in laser materials; J. App/, 
Phys. 42, No. 2, 619-630 (Feb. 1971). 


STEELS GI ELY ; 
Key words: Antimony; heat conduction; laser materials;, Nd-doped 
glass; platinum: stress components; and thermoelastic theory. 
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One of the severe problems encountered in high-power-solid-state 
laser systems is the thermal damage to laser rods and optical 
elements arising from metallic or dielectric inclusions; i.e., impuri- 
ties with physical and optical properties which differ substantially 
from those of the host material. Such inclusions may absorb an 
appreciable amount of the incident radiation and thereby may pro- 
duce major stresses within the host material. In this paper, the 
dependence of the maximum value of the tensile stress upon the 
size of the inclusion and upon the physical properties of the host 
is examined. The feasibility of using optical techniques to detect 
metallic and dielectric inclusions in laser materials before they cause 
damage also is studied. The computations suggest that the use of 
laser pulse widths of the order of microseconds or longer may be 
more promising for the detection of small incipient absorbing cen- 
ters than the use of nanosecond pulse widths. 


Bennett, L. H., Swartzendruber, L. J., Interactions between 
small magnetic clusters in copper-rich Cu-Ni-Fe alloys, /. 
Appl. Phys. 42, No. 4, 1547-1548 (March 15, 1971). 


Key words: Alloys; clustering; copper: iron; local moments; Méss- 
bauer effect: nickel. 


Traces of iron in “nonmagnetic” Cu-rich Cu-Ni alloys produce small 
magnetic clusters largely limited to an Fe site and its Ni nearest 
neighbors, as observed in Mdéssbauer source experiments. The 
present work concerns *7 Fe Mossbauer absorber data for higher 
Fe concentrations (}%-3%), in the 10%-30% Ni range, which 
demonstrates the existence of a magnetic-moment distribution 
centered about two distinct moment values. For the lower Fe con- 
centrations, the smaller and more abundant of these moments is 
identical to that observed in the source experiments and is associ- 
ated with isolated small Fe-Ni clusters. The other moment is about 
twice as large and can be ascribed to coupled pairs of these small 
clusters. These results are almost free of effects of metallurgical 
clustering. By appropriate heat treatment, clustered samples are 
obtained with different moment distributions and higher Curie 
temperatures. 


Boyd, M. E., Mountain, R. D., Methods for determining the 
second virial coefficient of a gas from speed-of-sound data, 
Phys. Rev. A. 2, No. 5, 2164-2167 (November 1970). 


Key words: Acoustic thermometry; equation of state; equilibrium 
gas properties: helium; second virial coefficient; speed of sound in 
gases. 


Two methods of analyzing data on speed of sound in gases to obtain 
the second virial coefficient B (T) are compared. The older method, 
which assumes a form for the temperature dependence of B, is 
shown to correspond to finding an exact solution to an approximate 
differential equation for B, while a method recently proposed by 
Bruch solves the exact equation, but in an approximate manner. 
Examination of the errors in each method indicates that the first 
method is preferable. 


Boyne, H. S., Laser frequency stabilization techniques and 
applications, Proc. 24th Annual Frequency Control Symp., on 
Frequency Control, Atlantic City, N.J., April 27-29, 1970, pp. 
233-339 (U.S. Army Electronics Command, Fort Monmouth, N.J., 
1970), IEEE Trans. Instr. Meas. IM-20, No. 1, 19-22 (Feb. 1971). 


Key words: Frequency; laser; laser frequency measurement; stabi- 
lization; techniques; time standard. 


A review of progress in laser stabilization techniques and laser 
frequency measurement is given. Methods for relating laser frequen- 
cies to the time standard and methods for absolute laser frequency 
stabilization are described. Experimental information on repro- 
ducibility and noise characteristics is reported. Application to fre- 
quency and wavelength standards is discussed. 


Boyne, H. S., Hall, J. L., Barger, R. L., Bender, P. L., Ward, J., 
Levine, J., Faller, J., Absolute strain measurements with a 30 
meter vacuum interferometer, Proc. Conf. Laser Applications 
in the Geoscience, Huntington Beach, Calif., June 1969, pp. 215-225 
(1970). 





Key words: Earth strain; earth tide; geophysics; laser strainmeter: 
seismograph. 


We present details on the design and performance of a 30 meter 
interferometric strain gauge. We also discuss a practical method 
for recording absolute earth strain measurements by comparing 
length changes in the interferometer with an absolute wavelength 
standard. 


Breden, L. H., Technical Standards Coordinator, Fluorinated 
ethylene-propylene (FEP) plastic lined steel pipe and 
fittings, Nat. Bur. Stand. (US.), Prod. Stand. 34-70, 12 pages 
(Mar. 1971) 10 cents, SD Catalog No. C13.20/2:34-70. 


Key words: FEP lined pipe; fluorinated ethylene-propylene plastic 
lined pipe; pipe and fittings, steel, plastic lined; plastic lined steel 
pipe; steel pipe and fittings, plastic lined. 


This Voluntary Product Standard covers requirements and methods 
of test for the material, dimensions, construction, and performance 
of commercially available steel pipe and fittings lined with fluori- 
nated ethylene-propylene (FEP) plastic intended to be used for 
conveying acids, gases, solvents, and other corrosive materials. 


Breden, L. H., Technical Standards Coordinator, Poly(vinyl- 
shrinkable film, Nat. Bur. Stand. (U.S.), Prod. Stand. 29-70, 
16 pages (Jan. 1971) 15 cents, SD Catalog No. C13.20/2:29-70. 


Key words: Balanced and unbalanced film; film, plastic; orientation 
release stress; plastic heat-shrinkable film; retraction ratio; shrink 
film, plastic. 


This Voluntary Product Standard covers plastic film which, within 
its useful temperature range, will shrink at least 20 percent in either 
or both the machine and the transverse directions. Included are 
requirements and methods of test for shrink properties, dimensions, 
yield, weight, and slip. 


Breden, L. H., Technical Standards Coordinator, Poly(vinyl- 
idene fluoride) (PVF.) plastic lined steel pipe and fittings, 
Nat. Bur. Stand. (U.S.), Prod. Stand. 35-70, 12 pages (Mar. 197i) 
10 cents, SD Catalog No. C13.20/2:35-70. 


Key words: Pipe and fittings, steel, plastic lined; plastic lined steel 
pipe; poly(vinylidene fluoride) plastic lined pipe; PVF. lined pipe; 
steel pipe and fittings, plastic lined. 


This Voluntary Product Standard covers requirements and methods 
of test for the material, dimensions, construction, and performance 
of commercially available steel pipe and fittings lined with poly- 
(vinylidene fluoride) (PVF) plastic intended to be used for convey- 
ing acids, gases, solvents, and other corrosive materials. 


Breden, L. H., Technical Standards Coordinator, Polystyrene 
plastic sheet, Nat. Bur. Stand. (U.S.), Prod. Stand. 31-70, 14 
pages (January 1971) 15 cents, SD Catalog No. C13.20/2:31-70. 


Key words: Impact resistance, plastic sheet; oriented sheet; plastic 
sheet; polystyrene plastic sheet; sheet, plastic. 


This Voluntary Product Standard covers requirements and methods 
of test for dimensions, and mechanical and physical properties 
of hot-melt extruded polystyrene plastic sheet. The standard does 
not cover polystyrene foam or biaxially oriented sheet. 


Breden, L. H., Technical Standards Coordinator, Polytetra- 
fluoroethylene (PTFE) plastic lined steel pipe and fittings, 
Nat. Bur. Stand. (U.S.), Prod. Stand. 33-70, 12 pages (Mar. 1971) 
10 cents, SD Catalog No. C13.20/2:33-70. 


Key words: Pipe and fittings, steel, plastic lined; plastic lined steel 
pipe; polytetrafluoroethylene plastic lined pipe; PTFE lined pipe; 
steel pipe and fittings, plastic lined. 


This Voluntary Product Standard covers requirements and methods 
of test for the material, dimensions, construction, and performance 
of commercially available steel pipe and fittings lined with poly- 
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tetrafluoroethylene (PTFE) plastic intended to be used for conveying 
acids, gases, solvents, and other corrosive materials. 


Brenner, A., Note on an organic-electrolyte cell with a high 
voltage, J. Electrochem, Soc. 118, No. 3, 461-462 (March 1971). 


Key words: Graphite electrode; high energy density cell; lithium 
deposition; organic electrolyte. 


An unusually high back EMF of 5.2 volts was obtained after electroly- 
sis of a 10% solution of lithium fluoborate in sulfolane with a graphite 
anode. 


Brenner, A., Note on the electrodeposition of magnesium 
from an organic solution of a magnesium-boron complex, 
J. Electrochem. Soc. 118, No. 1, 99-100 (January 1971). 


Key words: Decaborane; electrodeposition of magnesium; lithium 
deposition; nonaqueous electroplating; organometallic compound. 


A magnesium electroplating bath which yields smooth, white de- 
posits, has been prepared from lithium methyl, decaborane, and 
magnesium chloride. The solvent consisted mostly of tetrahydrofuran. 


Brenner, F. C., Kondo, A., Research for a uniform quality 
grading system for tires. IV. Tread wear, Rubber Chem. Technol. 
44, No. 1, 106-121 (March 1971). 


Key words: Automobile tires; grading; rating; test method; tread 
wear. 


This research had as its objective the development of a test pro- 
cedure for rating tires for tread wear. The results of preliminary 
experiments on roads of different composition and severity are 
discussed. The authors conclude that the course on which tires 
are to be rated should include several different pavements and a 
variety of maneuvers of unequal severity. 

Two courses, one near Bryan, Texas, and the other near San Antonio, 
Texas, were designed on this principle. Eighteen brands of tires 
including bias, radial, and G78 belted bias tires were run on each 
course. The tires ranked in order of decreasing projected mileage on 
each course were highly correlated, although the mileages were 
not eal, 


Brown, D. W., Lowry, R. E., Wall, L. A., The radiation-induced 
polymerization of 3,3,4,4,5,5,5-heptafl pent 1 at 
high pressure, J. Poly. Sci. 8, Part A—1 , 3483-3493 (1970). 





Key words: 3,3,4,4,5,5,5-Heptafluoropentene-1; irradiation; polym- 
erization; pressure. 


A study was made of the radiation-induced polymerization under 
pressure of 3,3,4,4,5,5,5-heptafluoropentene-l. Polymerization rates 
increase with pressure (activation volume equals — 1] cm*/mol) and 
temperature (activation enthalpy equals 6.5 kcal/mol) in liquid phase. 
At 13800 atm and 25°C freezing occurs; the polymerization rate in 
the solid is very small. In liquid phase polymerization can continue 
for many hours after the irradiation is terminated. An active species 
is formed by radiation which initiates polymerization in the dark 
period. 


Brown, R. L., Diffusion of a trace gas into a flowing carrier, 
Intern. J. Chem. Kinetics 2, 475-477 (1970). 


Key words: Diffusion equation; diffusion of gases; flow reactor; gas 
mixing; mixing by diffusion; relaxation methods. 


Numerical methods are used to solve the differential equation for 
diffusion of a trace gas into a flowing carrier gas having a parabolic 
velocity profile in a cylindrical tube. Steady state solutions are given 
in the form of contour diagrams of constant trace gas concentration. 


Bullis, W. M., Baroody, A. J., sr., Editors, Methods of measure- 
ment for semiconductor materials, process control, and 
devices, Quarterly Report, April 1 to June 30, 1970, Nat. Bur. Stand. 


(U.S.), Tech. Note 560, 58 pages (Nov. 1970) 60 cents, SD Catalog No. 


C13.46:560. 


Key words: Alpha-particle detectors; aluminum wire; carrier 
lifetime; die attachment; electrical properties; epitaxial silicon; 
gamma-ray detectors; germanium; gold-doped silicon; metallization; 
methods of measurement; microelectronics; microwave devices; 
nuclear radiation detectors; resistivity; semiconductor devices; 
semiconductor materials; semiconductor process control; silicon; 
thermal resistance; thermographic measurements; ultrasonic bonder; 
wire bonds. 


This quarterly progress report, eighth of a series, describes NBS 
activities directed toward the development of methods of measure- 
ment for semiconductor materials, process control, and devices. 
Principal emphasis is placed on measurement of resistivity, carrier 
lifetime, and electrical inhomogeneities in semiconductor crystals; 
evaluation of wire bonds, metallization,adhesion, and die attachment; 
and measurement of thermal properties 6f semiconductor devices 
and electrical properties of microwave devices. Work on related 
projects on silicon nuclear radiation detectors and specification of 
germanium for gamma-ray detectors is also described. Supplemen- 
tary data concerning staff, standards committee activities, tech- 
nical services, and publications are included as appendixes. 


Bullis, W. M., Editor, Methods of measurement for semicon- 
ductor materials. process control, and devices, Quarterly 
Report, July 1 to September 30, 1970, Nat. Bur. Stand. (U.S.), Tech. 
Note 571, 58 pages (Apr. 1971) 60 cents, SD Catalog No. C13.46:571. 


Key words: Alpha-particle detectors; aluminum wire; base transit 
time; carrier lifetime; die attachment; electrical properties; epi- 
taxial silicon; gamma-ray detectors; germanium; gold-doped silicon; 
metallization; methods of measurement; microelectronics; micro- 
wave devices; nuclear radiation detectors; probe techniques (a-c); 
resistivity; semiconductor devices; semiconductor materials; semi- 
conductor process control; silicon; thermal resistance; thermo- 
graphic measurements; ultrasonic bonder; wire bonds. 


This quarterly progress report, ninth of a series, describes NBS 
activities directed toward the development of methods of measure- 
ment for semiconductor materials, process control, and devices. 
Work is continuing on measurement of resistivity, carrier lifetime, 
and electrical inhomogeneities in semiconductor crystals; specifica- 
tion of germanium for gamma-ray detectors; evaluation of wire 
bonds, metallization adhesion, and die attachment; measurement of 
thermal properties of semiconductor devices and electrical proper- 
ties of microwave devices; and characterization of silicon nuclear 
radiation detectors. New effort is being started on the measurement 
of transit-time and related carrier transport properties in junction 
devices. Supplementary data concerning staff, standards committee 
activities, technical services, and publications are included as 
appendixes. 


Burdick, M. D., An instrumental procedure for evaluating 
adherance of porcelain enamel cover coats direct-to-steel, 
Proc. Porcelain Enamel Institute Technical Forum, University of 
Illinois, Urbana, Ill., Oct. 7-9, 1970, 32, 160-167 (Porcelain Enamel 
Institute, Washington, D.C., 1971). 


Key words: Adherence; instrumental evaluation; porcelain enamel 
direct-to-steel. 


An instrumental method for rating the quality of direct-on adherence 
is described. The method is appropriate, either for monitoring 
day-to-day production or for making interlaboratory comparisons. 
The rating method is firmly based on present practice, is rapid and 
does not require the somewhat subjectuve judgement of trained 
observers. 


Carter, G. C., Kahan, D. J., Bennett, L. H., Cuthill, J. R., Dobbyn, 
R. C., The NBS alloy data center: permuted materials index, 
Nat. Bur. Stand. (U.S.), Spec. Publ. 324, 683 pages (Mar. 1971) 
$7.00, SD Catalog No. C13.10:324. 


Key words: Alloy data; bibliography; index; information; Knight 
shifts; NMR; soft x-ray. 


This Index contains literature references to ~ 10,000 research 
papers on physical properties of metals and alloys. The Index 
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contains all NMR Knight shift papers and soft x-ray emission papers. 
It also contains many soft x-ray absorption papers and a number 
of papers on generally related topics such as susceptibilities, spe- 
cific heats, hyperfine fields, and band structures. The papers are 
annotated in depth and the coded information put onto a magnetic 
tape. The Permuted Materials Index was created from this tape, 
listing alloys under each of their constituent components (i.e., 
CuNi appears under CuNi and under NiCu alloys). 


Cassidy, M. D., Calibration and test services of the National 
Bureau of Standards, Nat. Bur. Stand. (U.S.), Spec. Publ. 250— 
1970 Edition, 209 pages (Dec. 1970) $2.00, SD Catalog No. C13.10:250. 
Supersedes NBS Miscellaneous Publication 250 — 1968 Edition. 


Key words: Calibration services; fees for testing; test fee schedule; 
test procedures. 


This publication contains a descriptive listing, item by item, of most 
of the test and calibration work done at the National Bureau of Stand- 
ards with the respective fees. The content is fully indexed providing 
easy access to data on specific measurement services. Included also 
is information about the Bureau’s organization and statutory func- 
tions, the procedures for requesting tests, and the reports issued 
thereon. Throughout the text are mentioned various other publi- 
cations which form a bibliography in measurement techniques. This 
publication supersedes Special Publication 250—1968 Edition 
and the supplementary material in Measurement Users Bulletins 
Numbers | through 11. 


Cezairliyan, A., A high speed method of measuring thermal 
expansion of electrical conductors, Rev. Sci. Instr. 42, No. 
4, 540-541 (April 1971). 


Key words: High-speed measurements; high temperature; platinum; 
thermal expansion. 


A transient method for the measurement of thermal expansion of 
electrical conductors is described. The method is based on detecting 
the change in radiance coming from a constant radiation source as a 
result of the expansion of the specimen placed between the radiation 
source and a radiation detecting system. The specimen can be pulse 
heated from room temperature to near its melting point in less than 
one second and pertinent experimental quantities can be measured 
with a time resolution of 0.4 ms and a full-scale signal resolution of 
one part in 8000. To check the method, preliminary experiments 
were performed on platinum in the temperature range 300 to 700 K. 
The estimated inaccuracy of the results is within 5 percent. The 
agreement of the results with those in the literature is within 3 
percent. 


Cezairliyan, A., High-speed methods of measuring thermo- 
physical properties at high temperatures (Proc. Conf. Thermo- 
physical Properties, Manchester, England, April 7-10, 1970), Rev. 
Inst. Hautes Temper. Refract T.7, 215-229 (1970). 


Key words: High-speed measurement methods; high temperatures; 
thermodynamics; thermophysical properties. 


In this paper, “high-speed”’ refers to experiments which are of sub- 
second duration. However, in some cases quasi-dynamic experi- 
ments with durations greater than one second are also included. Ad- 
vantages of high-speed measurement of thermophysical properties of 
substances at high temperatures (above 2000 K) are presented. 
Requirements of high-speed measurements are given. Methods 
used for the generation and measurement of thermal energy in short 
times are described. Various techniques (photoelectric, photographic, 
and others) used for the measurement of the temperature of rapidly 
heating or cooling specimen are presented. Techniques for high- 
speed recording of quantities are described. Particular attention is 
given to millisecond-resolution digital data acquisition systems. 
Application of the high-speed techniques to the measurement of 
specific heat, electrical resistivity, thermal radiation properties, 
melting point, thermal diffusivity and others are described. Accuracy 
of high-speed measurements are discussed and are compared with 
those of conventional methods. Potential applications of high-speed 
methods are presented. 


Cezairliyan, A., Short-time (ms) fluctuations of radiance tem. 
perature of graphite arc, Appi. Opt. 10, No. 5, 1178 (May 1971). 


Key words: Graphite arc; high-speed measurements; high tempera- 
tures; pyrometry; radiance standards. 


Short-time fluctuations of the radiance temperature graphite arcs 
are investigated with a millisecond resolution photoelectric pyrom- 
eter. Standard deviation of the fluctuations around their mean is 
found to be 4 K. The range of low and high temperatures for four ex- 
periments extended from 3799 K to 3808 K, with an average of 3804 K 
on IPTS-68. 


Chamberlain, G. E., Mielezarek, S. R., Kuyatt, C. E., Absolute 
measurement of differential cross sections for electron 
scattering in helium, Phys. Rev. A2, No. 5, 1905-1922 
(November 1970). 


Key words: Absolute differential cross section; Born; elastic; elec- 
tron impact; helium; inelastic. 


Absolute measurements of cross sections for electron-impact 
scattering in helium at 5° angle have been made for elastic scatter- 
ing and excitation of the 2'P and 2'S states and for incident energies 
of 50-400 eV. Cross section values of o(2'P,5°) are found to be lower 
than theoretical Born values by 10% at 400 eV, 30% at 100 eV, and 
60% at 50 eV. Deduced values of total 2'P excitation cross sections 
(E = 100 eV) are in agreement with other experimental values. Our 
measurements for elastic scattering agree well with recent theo- 
retical calculations. 


Chandler, H. H., Rupp, N. W., Paffenbarger, G. C., Poor mercury 
hygiene from ultrasonic amalgam condensation, /. Am. 
Dental Assoc. 82, 553-557 (March 1971). 


Key words: Amalgam 
levels; threshold limit. 


condensation; mercury; mercury vapor 


There is a potential hazard of mercury poisoning in using an ultra- 
sonic device for amalgam condensation. A cloud of mercury droplets 
and alloy particles was emitted from the soft amalgam at the work- 
ing tip of the instrument. Mercury vapor levels as recorded by a 
vapor detector held 30 cm from the working tip were 20% of the 
allowable threshold limit value (0.1 mg Hg/m# of air) and probably do 
not represent unsafe levels. The continued long-time use of the ultra- 
sonic instrument would result in the deposition of a great many 
mercury droplets throughout a dental operatory and could thereby 
cause higher mercury vapor levels especially in poorly ventilated 
spaces. In addition, the inhalation of the emitted material by the 
patient and the dental health personnel cannot be considered 
good hygiene. Therefore, the use of this instrument for amalgam 
condensation is contra-indicated until such time as the safety of the 
instrument for this purpose is firmly established. 


Chappell, S. £., Humphreys, J. C., Silicon detector measure- 
ments of energy deposition in aluminum by monoenergetic 
electrons (Proc. Annual Conf. Nuclear Space Radiation, La Jolla, 
Calif., July 1970), IEEE Trans. Nuclear Science NS-17, No. 6, 
272-277 (Dec. 1970). 


Key words: Absorbed energy vs depth; absorbed-energy distribu- 
tions; incident monoenergetic electrons; Monte Carlo calculations; 
semi-infinite aluminum medium; silicon detector. 


The energy deposited at various depths in aluminum by incident 
monoenergetic electrons has been measured with a silicon semi- 
conductor, transmission detector. Beams of monoénergetic elec- 
trons with incident energies of 0.50, 0.75, and 1.0 MeV were directed 
normally on a semi-infinite slab of aluminum in which a 0.196-mm 
silicon detector was positioned at various depths. The pulse-height 
distributions recorded with the detector were converted to absorbed- 
energy distributions from which the probability of energy absorp- 
tion per incident electron in the specific layer, as well as the absorbed 
energy as a function of depth in the material, could be determined. 
The curves of absorbed energy as a function of depth obtained for 
aluminum at each energy were compared to those calculated 
by Berger and Seltzer, employing a Monte Carlo method. Good 
agreement is shown between calculations and measurements. 
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Cohen, J., Edelman, S., Direct piezoelectric effect in poly- 
vinylchloride films, J. Appi. Phys. 42, No. 2, 893-894 (February 
1971). 


Key words: g Modulus; piezoelectricity; polymers; polyvinyl- 
chloride films. 


The piezoelectric modulus, g3:, of polyvinylchloride films has been 
measured at 20 Hz. Values are in the range of approximately 0.4 to 
3 Vm-'/Nm-?, which are larger than the g values of conventional 
piezoelectric materials. 


DeVoe, J. R., Spijkerman, J. J.. Méssbauer spectrometry, Anal. 
Chem. Annual Reviews 42, No. 5, 366R—388R (April 1970). 


Key words: Chemical applications: literature; Méssbauer spec- 
troscopy; review. 


A review of the literature on chemical applications of Méssbauer 
Spectrometry for 1968 and 1969 are presented. This is done pri- 
marily with a table of pertinent information on compounds and tech- 
niques. New developments in the field are also presented. 


Dikkers, R. D., Marshall, R. D., Thom, H. C. S., Hurricane 
Camille— August 1969, Nat. Bur. Stand. (U.S.), Tech. Note 569, 
71 pages (Mar. 1971) 70 cents, SD Catalog No. C.13.46:569. 


Key words: Buildings; failure; hurricanes; mobile homes; roofs; 
structural engineering: tides: wind. 


One week after Hurricane Camille devastated the Mississippi- 
Louisiana Gulf Coast with 125 mph winds and 20-ft tides on August 
17, 1969, a four-man NBS team investigated the damage to buildings 
and other structures. This report presents photographic survey data 
from this investigation along with additional data on wind speeds 
and storm surge. Based on these data, suggestions are made pertain- 
ing to the improvement of building design and construction practices. 


DiMarzio, E. A., Guttman, D. M., Separation by flow and its 
application to gel permeation chromatography, J. Chromatog. 
55, 83-97 (1971). 


Key words: Chromatography; elution volume; flow; gel permeation 
chromatography; hydrodynamic volume; polymer; polymer separa- 
tion. 


When dilute solutions of finite size particles undergoing Brownian 
motion flow through a capillary, the larger particles have higher 
average velocities than the smaller particles. Thus one can obtain a 
separation of particles of different sizes due to fluid flow. The 
elution volumes of suspended particles or polymer molecules are 
derived for various tube geometries. Following Taylor, the effects 
of diffusional broadening of the volume elution peak for finite size 
particles are discussed and the process is shown to be chromato- 
graphic. 

Models of a GPC column are proposed in which there is fluid and 
particle flow through each of the beads as well as around them. 
Diffusion is allowed within and outside of the beads. Equations for 
the location of the volume elution peaks are computed for such 
models and shown to yield functional dependence on the polymer 
radius and column geometry very much like equations derived by 
previous workers for models of GPC columns in which there was no 
flow allowed within the beads. Explicit formulae are given for the 
second and third moments for the above models. It is shown that for a 
monodisperse species the volume elution peak is always a gaussian 
of a finite width. It is shown that beads with open pores that allow for 
flow always have better separation capabilities than beads with 
pores that do not allow for flow. 


Dziuba, R. F., Dunfee, B. L., Resistive voltage-ratio standard 
and measuring circuit, /EEE Trans. Instr. Meas. IM-19, No. 4, 
266-277 (Nov. 1970). 


Key words: Divider; measurement system; ratio; ratio standard; 
resistive divider; voltage ratio. 


This paper describes a highly stable, guarded d-c voltage ratio stand- 
ard and the measuring network and techniques used to establish 
the values of its ratios to an accuracy of 0.2 ppm. The entire system is 
housed within a dry air enclosure whose temperature is maintained 
at 23 + 0.05°C. Discrete ratios from 1/1, to 1000/1, are provided 
with maximum rated voltage set at 1000 volts. The design of the 
standard was chosen so that a redundancy of measurment could be 
incorporated in the system. Thus it is possible to measure each 
successive ratio by a substitution or “boot-strap’’ method, and, by 
satisfying the conditions of the series-parallel principle, to measure 
the 10/1, 100/1, and 1000/1 ratios by a second independent method. 
The design also admits additional checks on the validity of the meas- 
urements. An analysis of measurement errors and a discussion of 
their possible origin are included. Since the intent also was to design 
the ratio standard for low frequency operation some preliminary 
data is included on its ac performance. 


Edgerly, D. E., Weights and Measures Labeling Handbook, 
Nat. Bur. Stand. (U.S.), Handb. 108, 406 pages (May 1971) $3.25, SD 
Catalog No. C13.11:108. 


Key words: Consumer package; exemption; fluid measure; labeling; 
net quantity; principal display panel; regulations; weights and 
measures. 


The Weights and Measures Labeling Handbook has been prepared 
for use by State weights and measures officials in the enforcement of 
State laws and regulations pertaining to the packaging and labeling 
of commodities. The Handbook provides a cross indexed compilation 
of the requirements of the Model State Packaging and Labeling Regu- 
lation, U.S. Department of Agriculture regulations, and regulations 
of the Food and Drug Administration and Federal Trade Commission 
issued pursuant to the Fair Packaging and Labeling Act. 


Engen, G. F., A new method of characterizing amplifier noise 
performance, /EEE Trans. Instr. Meas. IM-19, No. 4, 344-349 
(Nov. 1970). 


Key words: Amplifier noise; noise; noise factor; noise temperature; 
Y-factor. 


Although the use of noise figure or noise temperature to character- 
ize amplifier performance is a well-established practice, it is also 
recognized that this parameter provides only a partial description of 
the amplifier noise properties. In general, the noise figure (or 
temperature) depends upon the generator inpedance and is thus a 
function of the signal-source and amplifier combination. 

Typically, the noise figure is measured by the Y-factor method using 
hot and cold noise sources that are nominally matched (reflection- 
less). The result of this measurement is of value as a figure of 
merit; however, if optimum performance is to be realized, the 
applications engineer must know whether to adjust the signal source 
impedance for maximum power transfer, minimum noise figure, 
or some other criterion, and he must know the deterioration in per- 
formance that results if this is not done. It is the purpose of this paper 
to present an alternative method of characterizing amplifier noise 
performance in terms of parameters that provide ready answers to 
these questions. In addition, the measurement of these parameters 
via a simple extension of the Y-factor method will be described. 


Evenson, K. M., Radford, H. E., Moran, M. M., Jr., CH free radicals. 
detected by infrared lasermagnetic resonance, Appl. Phys. 
Letters, 18, No. 10, 426-429 (May 15, 1971). 


Key words: CH; far infrared; free radical; laser; laser magnetic 
resonance. 


CH radicals in an oxyacetylene flame have been detected by a 
magnetic resonance absorption method, with a water-vapor laser as 
the source oscillator at a frequency of 2.5 THz. The absorption spec- 
trum is identified with the pure rotational transition, J, K=(5/2,2)—> 
(7/2,3) at a wavelength of 118.6 wm. 


Evenson, K. M., Wells, J. S., Matarrese, L. M., Jennings, D. A., 
Variable output-coupling far-infrared Michelson laser, 
J. Appl. Phys. 42, No. 3, 1233-1234 (March 1, 1971). 
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Key words: HCN lasers; H2O lasers; infrared lasers; variable laser 
coupling. 


A method for varying the output coupling of far-infrared gas lasers is 
described. It is based on the Michelson interferometer principle and 
has the advantages of increasing the power available from the laser 
and providing a smooth variation of the output coupling. 


Farabaugh, E. N., Mauer, F. A., Application of a modified 
technique for x-ray topography of large crystals, J. Crystal 
Growth 7, 282-284 (1970). 


Key words: Laser rod; orientation; ruby; single crystal; topography 
x-ray. 


A transmission x-ray topography technique has been developed 
which is suitable for examining large, uncut crystals. A relatively 
large volume is represented in each photograph, and regions that 
are misoriented by more than 0.3° can be detected. Examination of 
single crystal ruby laser rods is used to illustrate the capabilities 
and limitations of the method. 


Fatiadi, A. J., Determination of inososes with an alkaline solu- 
tion of copper (II) oxalate-tartrate complex (the Somogyi 
reagent) and reaction mechanisms involved, Carbohydrate Res. 
17, 419-430 (March 1971). 


Key words: Determination; electron-transfer; inosose; oxidation; 
quantitative; radical; reagent. 


Four inososes have been analyzed with the Somogyi reagent and em- 
pirical equations for their quantitative determination were derived 
because reaction of the Somogyi reagent with inososes affords non- 
stoichiometric quantities of cuprous oxide, and each inosose has a 
different reducing power. 

Results from spectrophotometric and e.s.r. studies of the mechanism 
of oxidation of inososes with the Somogyi reagent at 25 to 55° are in 
agreement with a one-electron transfer process; however, at 90 to 
106°, extensive degradation of inososes by the Somogyi reagent oc- 
curs, doubtless caused by generation of transient radicals during the 
oxidation, as evidenced by results of a radical-scavenging experiment. 


Fickett, F. R., Resistivity of polycrystalline aluminum and 
copper in high magnetic fields: The effect of temperature 
and purity, Appl. Phys. Letters 17, No. 12,525-527 (Dec. 15, 1970). 


Key words: Aluminum; copper: magnetoresistance. 


Data are presented on the resistivity of polycrystalline aluminum and 
copper at 40 kOe and at temperatures from 4 K to 30 K. Specimen 
purity varies over three decades of residual resistance ratio. For 
either metal, the actual resistivity measured in the field at a given 
temperature decreases with increasing specimen purity. This result 


is important for proposed high magnetic field applications of these 
metals. 


Franklin, A. D., Bennett, H. S., Editors, ARPA-NBS Program of 
research on high temperature materials and laser materials, 
January | to June 30, 1970, Nat. Bur. Stand. (U.S.), Tech. Note 565, 
84 pages (Jan. 1971) 75 cents, SD Catalog No. C13.46:565. 


Key words: Aluminum oxide; band structure; borides; bulk optical 
properties; chemical analysis; chromium ions; corundum; crystal 
defects; crystal growth; density of states; diffusion; glass; high 
temperature materials; laser damage; laser degradation; laser glass: 
laser materials; mass transport; neodymium glass; optical properties; 
opto-elastic properties; oxides; ruby; rutile; titanium; titanium 
dioxide; transition metal. 


Brief reviews are given of work performed during this reporting 
period on several problems in the areas of High Temperature 
Materials and Laser Materials. Under the High Temperature 
Materials extremely pure AlO; crystals have been prepared and 
analysed chemically with the aim of providing specimens for re- 
search on mass transport, electronic, and similar properties. Oxygen 
diffusion data on TiQs, obtained by the gas uptake method, agree 
with literature values, and work is now proceeding on developing a 





sectioning technique for oxygen diffusion. The optical constants of 
titanium have been determined over the visible portion of the 
spectrum. Soft x-ray spectra, nmr data, Mossbauer spectra, magnetic 
susceptibility, and other probes of the electronic density of states 
have been determined for a group of diborides of the first-series 
transition metals (Sc, Ti, V, Cr, and Mn) and interpreted in terms of 
the d-electron contribution to the density of states near the Fermi 
level. Both this project and the previous one on the optical properties 
of titanium have been completed. Under Laser Materials, the 
measurement of bulk optical properties and chemical analyses of 
ruby and Nd-glass are complete. Several techniques for determining 
the Cr+ content in laser ruby have been developed and shown to 
agree; with these, it has been shown that very little loss of Cr+ 
occurs when ruby turns orange under optical (x-ray or y-ray) irradia- 
tion. A theory of laser pulse-induced damage at foreign particles in 
Nd-glass has been worked out; it suggests that detection before 
damage might be easiest with pulses of microsecond or longer 
duration. A high-power Nd-glass laser, Q-switched, has been 
purchased for damage studies. The beam characteristics, which 
must be properly controlled for meaningful damage threshold 
measurements, are currently being studied. 


Furcolow, W. H., Technical Standards Coordinator, Clinical 
thermometers (Maximum-self-registering, mercury-in-glass), 
Nat. Bur. Stand. (U.S.), Prod. Stand. 39-70, 12 pages (May 1971) 
15 cents, SD Catalog No. C13.20/2:39-70. 


Key words: Clinical thermometers: glass thermometers, clinical: 
mercury-in-glass thermometers; thermometers,  self-registering, 
clinical. 


This Voluntary Product Standard covers the requirements and 
methods of testing maximum-self-registering, mercury-in-glass 
thermometers of the types commonly used for measuring body 
temperatures, such as oral and rectal types in both regular and 
basal temperature scales. It is intended to serve as a nationally 
recognized basis for certification of compliance by manufacturers 
and for procurement purposes by consumers. The standard includes 
requirements for bulb and stem glasses, mercury, dimensions, 
temperature scale ranges, and graduations, and performance criteria 
for thermometer aging, hard shaking determination, and accuracy 
of scale reading. 


Furcolow, W. H., Technical Standards Coordinator, Hinged 
interior wood door units, Nat. Bur. Stand. (U.S.), Prod. Stand. 
32-70, 15 pages (Mar. 1971) 30 cents, SD Catalog No. C13.20/2:32-70. 


Key words: Door units, wood; hinged door units, interior: pre- 
fabricated door units; wood door units. 


This Voluntary Product Standard covers the designs, grades, sizes, 
materials (including hardware), and construction of prefabricated 
wood door units, as well as a means of identifying units represented 
as conforming to the standard, and definitions of terms used in the 
standard. 


Gadzuk, J. W., Plummer, E. W., Energy distributions for 
thermal field emission, Phys. Rev. B. 3, No. 7, 2125-2129 (April 1, 
1971). 


Key words: Field emission; surfaces and tunneling. 


A sequence of total energy distribution curves for field emission was 
experimentally obtained for a tungsten emitter heated to 1570 K. 
Theoretical curves using the Miller-Good WKB-type approximation 
for tunneling probabilities are in good agreement with the experi- 
mental measurements. A significant feature of both sets of curves 
is a change in slope corresponding to electron emission near the 
top of the surface where the emission mechanism changes from 
tunneling to thermionic emission. This feature is in accord with 
the classical-image force model for the surface potential which 
appears to be valid for distances approaching 3-4 A to the metal 
surface. 


Garner, E.L. Machlan, L. A., Shields, W. R., Standard Reference 
Materials: Uranium isotopic standard reference materials, 
(Certification of uranium isotopic standard reference materials), 
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Nat. Bur. Stand. (U.S.), Spec. Publ. 260-27, 162 pages (April 1971) 
$1.25, SD Catalog No. C13.10:260-—27. 


Key words: Absolute isotopic abundance; ignition procedure; 
isotopic standards, mass spectrometry; stoichiometry; uranium. 


An ignition procedure has been developed that will yield reproducible 
stoichiometry for U;0, The effects of temperature, length of 
ignition, rate of cooling, pressure and type of atmosphere were 
investigated. This ignition procedure has been used for the blending 
of high purity *5U and 2°8U separated isotopes to prepare calibration 
standards for the determination of bias effects in the thermal 
emission mass spectrometry of uranium. Weight aliquoting was 
used to prepare calibration mixes with 7 U/*8U ratios of more than 
10 and less than 0.1 and to add a 2°3U spike for the determination of 
minor isotope abundances in the uranium isotopic standards by 
the isotope dilution technique. 

A description of the unique features of the mass spectrometer in- 
strumentation including the source, NBS collector and expanded 
scale recorder are given. Two specific analytical procedures were 
used for the isotopic analysis of uranium and are adaptable, within 
a general framework, to fit the particular ion current intensity 
requirements of a wide range of isotopic distributions. Mass dis- 
crimination due to evaporation and ionization on the filaments, and 
other parameters such as temperature, time, sample size, sample 
mounting, total sample composition, acidity, filament material, 
pressure, non-ohmic response, R-C response and source memory 
were studied as part of the development effort to establish sound 
analytical procedures. 

The absolute isotopic abundances of 18 uranium SRMs were deter- 
mined by thermal emission mass spectrometry. The general approach 
was to determine absolute 25U/?*U ratios by using calibration mixes 
to correct for filament bias. Then the absolute ***U and #°6U were 
determined by 2*U isotope dilution. For SRM U-—0002, isotope dilu- 
tion was the only practical means of determining the low abundance 
of 5U as well as the **U. The limits given for the isotopic composi- 
tion of the uranium SRMs are at least as large as the 95 percent 
confidence limits for a single determination and include terms for 
inhomogeneities of the material as well as analytical error. 


Hahn, T. A., Thermal expansion of copper from 20 to 800 
K—Standard reference material 736, J. Appl. Phys. 41, No. 13, 
5096-5101 (December 1970). 


Key words: Components of error (within and between samples); 
copper; interferometer; standard reference material; thermal 
expansion. 


Copper is the first of a series of materials that will be certified as 
thermal expansion standards by the National Bureau of Standards. 
The results of tests on five specimens indicate the stock is of con- 
sistent quality so that it may be certified as Standard Reference 
Material 736. A Fizeau interferometer was used for the expansion 
measurements. Above room temperature a controlled atmosphere 
furnace using a calibrated Pt vs. Pt-10% Rh thermocouple was 
used. Below room temperature a cryostat capable of operation 
with both liquid nitrogen and helium was used with a calibrated 
platinum resistance thermometer. Values of expansivity were cal- 
culated between equilibrium temperatures. The expansivity was 
used in the analysis of the data. Third order polynomials were fit 
to the data for each of the five specimens in the overlapping tem- 
perature ranges from 0 to 70 K, 50 to 270 K, and 210 to 800 K to 
test for variations between the specimens. The deviations between 
the five equations were well within the standard deviations of the 
data for each of the specimens in the respective temperature 
intervals. All the expansivity data was then pooled and used to 
obtain an equation for each of the temperature ranges given above. 
These equations and their integrals were used to calculate the 
final values of expansivity and expansion respectively. The results 
of the statistical analysis of the expansion and expansivity data is 
presented. A comparison is made with the data in the literature. 


Hamer, W. J., Standard cells, the primary battery, George W. 
Heise and N. Corey Cahoon, Ed., 1, No. 12, 433-477 (John Wiley 
& Sons, Inc., New York, N.Y., January 1971). 


Key words: Electromotive force; saturated standard cells; standard 
cells; unsaturated standard cells. 


This paper gives a survey of saturated and unsaturated standard 
cells including descriptions of their construction and their be- 
havior under diversified conditions. The survey also includes a 
discussion of the role standard cells play in the maintenance of 
the unit of electromotive force. A discussion of the Hulett stand- 
ard battery is also included. 


Harman, G. G., Kessler, H. K., Application of capacitor micro- 
phones and magnetic pickups to the tuning and trouble 
shooting of microelectronic ultrasonic bonding equip- 
ment, Nat. Bur. Stand. (U.S.), Tech. Note 573, 24 pages (May 1971) 
35 cents, SD Catalog No. C13.46:573. 


Key words: Capacitor microphone; flip-chip; magnetic pickup; 
microelectronic interconnections; spider bonding; ultrasonic bond- 
ing; wire bonding. 


Microelectronic ultrasonic wire bonding equipment typically welds 
wires to integrated circuits at frequencies between 50 and 65 kHz. 
Mechanical vibrations at these frequencies are difficult to measure 
directly and malfunctions of the system may not be recognized. 
Two different methods of measuring these vibrations are described. 
The first method involves use of a capacitor microphone and a 
tapered tip, and the second method use of a small magnetic pickup. 
Procedures are given for establishing a specific ultrasonic vibration 
amplitude, tuning the ultrasonic system to resonance, and diagnos- 
ing both mechanical and electrical problems in wire bonding equip- 
ment. Although these techniques and procedures were developed 
for ultrasonic wire bonding equipment, they are applicable to other 
ultrasonic welding systems of lead attachment, such as flip-chip, 
beam lead and spider bonding. 


Harris, F. K., Fowler, H. A., Olsen, P. T., Accurate Hamon-pair 
potentiometer for Josephson frequency-to-voltage measure- 
ments, Metrologia 6, No. 4, 134—142 (October 1970). 


Key words: Frequency-to-voltage ratio; Josephson effect; potentiom- 
eter; voltage standard. 


Accuracy and precision of parts in 107 have been demonstrated in 
a potentiometer which compares two- to ten-millivolt de signals from 
Josephson junctions against the U.S. “legal volt” secondary standard 
cells. The circuit comprises two (10 100 ohm) Hamon boxes in 
series-to-parallel interchange, linked by a calibrated Kelvin-Varley 
current divider. The critical Hamon-circuit contacts are all made 
through mercury-wetted amalgams; the entire circuit operates in 
smoothly stirred oil. Observations indicate the pair of Hamon boxes 
to be initially balanced within 0.2 ppm at the operating temperature. 
Typical standard deviation in the day-to-day, run-to-run scatter of 
Josephson frequency-to-voltage ratio, measured with this instrument, 
is 0.3 ppm or less. Preliminary values for the observed ratio are in 
general accord with the quoted values of Parker et al., Petley and 
Morris, and Finnegan et al. 


Harvey, J. L., Milliken, L. T., Forthofer, R. J., Trends in motor 
vehicle brake fluids and their standards, Proc. Society of 
Automotive Engineering Conf., Paper No. 710253, pp. 1-17 (Jan- 
uary 11-15, 1971). 


Key words: Brake fluids; brake fluid standards; motor vehicle 
brake fluid; passenger car braking systems. 


The development of motor vehicle brake fluids (MVBF’s) and of 
their specification requirements is reviewed and discussed with 
emphasis on the major problems encountered in service. Those 
factors held to be of major importance in establishing performance 
requirements for MVBF’s are considered also, and applied to the 
major changes that have recently been proposed for revision of 
Federal Motor Vehicle Safety Standard No. 116. The needs and 
possibilities for further revision in the MVBF safety standard 
are examined. Finally, potential effects of some current trends in 
the design of passenger car braking systems are considered. 


Hellwig, H., Areas of promise for the development of future 
primary frequency standards, Proc. 24th Annual Symp. on 
Frequency Control, Atlantic City, N.J., April 27-29, 1970, pp. 
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246-258 (Electronic Industries Assn., Washington, D.C., 1970); 
Vetrologia 6, No. 4, 118-126 (October 1970). 


Key words: Accurate frequency standards; accurate length stand- 
ards; accurate time standards; accurate voltage standards; ion 
storage; laser stabilization; maser oscillator; metrology; quantum 
electronics; saturated absorption; slave oscillators. 


In this paper possibilities are discussed which may lead to the 
development of future primary frequency standards of superior 
accuracy capability. A review is given of the various methods and 
techniques which are currently employed in quantum electronic 
frequency standards or which have a potential usefulness. Various 
effects which influence the output frequency of a primary standard 
are associated with these methods. A classification is given: (a) 
Effects associated with the interrogation of particles (atoms or 
molecules); (b) effects related to the method of confining the parti- 
cles; and (c) effects associated with the particles themselves and 
the way in which they. are treated for an effective interrogation by 
electromagnetic radiation. These three classes of effects are dis- 
cussed in detail, and expectation values for the related uncertainties 
are given. For selected particles certain methods of interrogation, 
confinement, and particle preparation can be combined such as to 
minimize the net uncertainty due to all applicable effects. Differ- 
ent technical solutions are the result. A review of existing and 
proposed devices is given, including quantitative data on the stability 
and accuracy capability. Aspects of the most promising devices 
are discussed, and it is concluded that accuracy capabilities of 
10>" should be within reach of today’s research and development. 


Heydemann, P. L. M., A simple and dependable electrical 
feedthrough for high pressures, Rev. Sci. Instr. 41, No. 12, 
1896 (December 1970). 


Key words: Closure; feedthrough: high pressure; leads. 


This note describes a novel type of high pressure feedthrough 
capable of withstanding more than 25 kbar repeatedly. It can be used 
for measurements with DC requiring very high insulation resistance 
as well as for RF up to 30 MHz. 


Hicho, G. E., Yakowitz, H., Rasberry, S. D., Michaelis, R. E.. Stand- 
ard reference materials: A standard reference material 
containing, nominally four percent austenite, Vat. Bur. 
Stand. (U.S.), Spec. Publ. 260-25, 22 pages (Feb. 1971) 30 cents, 
SD Catalog No. C13.10:260-25. 


Key words: Austenite in ferrite; electron microprobe; powder metal- 
lurgy: quantitative microscopy; SRM; x-ray fluorescence analysis. 


This standard was produced by powder metallurgical techniques 
using known amounts of austenite. Using these techniques, 134 
specimens were prepared. Because these standards are expected to 
be used primarily for the calibration of x-ray diffraction equipment, 
only one surface of each standard is certified, and these surfaces 
range from 3.1 percent to 5.2 percent in austenite content. To make 
the specimens, 310 stainless steel powder (austenitic) was blended 
with 430 stainless steel powder (ferritic) to make a mixture of 5 
percent austenite in ferrite. The material was compacted, sintered, 
polished and etched so the austenite appears white and the ferrite, 
a deep brown. Then quantitative microscopy methods were used 
to determine the percentage of austenite near the surface. Further- 
more, the 310 powder contains 20 percent of nickel while the 430 
powder contains virtually no nickel. Therefore, after establishing 
a meaningful calibration curve, x-ray fluorescence analysis for the 
nickel content was also used as a direct measurement of the amount 
of austenite on the surface of the compact. Both procedures were 
carried out on fifteen specimens statistically selected from the total 
number of compacts produced. Agreement, within experimental 
error limits, was obtained between the x-ray fluorescence results 
and quantitative microscopy results. The x-ray fluorescence method 
was used to characterize all additional compacts. X-ray diffraction 
determinations of austenite content are in good agreement with 
the x-ray fluorescence and quantitative microscopy results. The 
compacts may be used as x-ray diffraction standards for austenite 
or in special cases as x-ray fluorescence standards for nickel content. 





Hoer, C. A., Agy, D. L., Broad-band resistive-divider-type 
directional coupler, /EEE Trans. Instr. Meas. IM-19, No. 4, 
336-343 (November 1970). 


Key words: Coupler; directional coupler; directivity; impedance 
measurements; power measurements. 


A new lumped circuit directional coupler is described for the | 
to 30 MHz frequency range. Directivity greater than 50 dB is achieved 
over this frequency range for coupling ratios of 30, 40, and 50 dB. 
Greater directivity can be achieved over a narrower frequency 
range, and the coupler can be tuned for essentially infinite direc- 
tivity at a fixed frequency. The couplers, although resistive, can 
pass 100, 300, or 1000 Watts for the coupling ratios of 30, 40, and 
50 dB respectively with less than 0.01 dB change in coupling ratio. 
The coupling ratios are constant with frequency to +0.25 dB. 
Equations and curves are given for obtaining component values 
needed to design a resistive divider type coupler with any coupling 


from 3 to 60 dB. 


Hord, J., Voth, R. O., Tabulated values of cavitation B-factor 
for helium, H:, No, F., O., refrigerant 114, and H.O, Nat. 
Bur. Stand. (U.5.), Tech. Note 397, 116 pages (Feb. 1971) $1.00, 
SD Catalog No. C13.46:397. 


Key words: Cavitation: cryogenics; pumps. 


A brief history is given on the development of the B-factor concept 
and its application to the design of liquid pumps. Adaptation of 
the “quasi-static” vaporization model to the cavitation process 
is discussed; previous methods of computing B-factor are reviewed 
and a simplified, more precise computation, consistent with the 
“quasi-static’”’ model, is established. Merits of the different com- 
putational techniques are discussed and two of the methods are 
graphically compared. The best available property data are used 
to compute B-factors for several fluids over a wide range of tempera- 
tures. The results are tabulated as reference data; they are useful 
in the application of the B-factor concept to the prediction of per- 
formance in cavitating liquid pumps. 


Hudson, R. D., Kieffer, L. J.. Compilation of atomic ultraviolet 
photoabsorption cross sections for wavelengths between 
3000 and 10 A, Atomic Data 2, No. 3, 205-262 (May 1971). 


Key words: Atoms; compilation; cross section data; photoabsorption: 
photoionization; photon. 


Graphical displays of selected experimental data on atomic photo- 
absorption cross sections as functions of wavelength are presented 
for 24 atoms. Criteria used in the data selection are discussed. 
References for all data chosen and an author index are included. 
The literature has been covered through October 1969. 


Ives, L. K., Ruff, A. W., Jr., Extended dislocation configura- 
tions in HCP silver-tin alloys of low stacking-fault energy, 


Metal Sci. J. 4, 201-209 (May 1970). 


Key words: Dislocation multipoles; dislocation networks; dislocation 
nodes; dislocations; electron microscopy; hexagonal alloys; silver-tin 
alloys. 


Configurations resulting from the interaction of extended disloca- 
tions in low stacking fault energy hcp silver-tin alloys are studied by 
means of transmission electron microscopy. The interactions of all 
possible sets of two different extended dislocations are considered. 
Since there are two different modes by which dislocations extend 
in the hep structure, these interactions fall into two classes: Class | 
interactions occur between two dislocations extended in the same 
way, and Class II interactions involve two dislocations extended in 
different ways. The resulting configurations consist of several 
different kinds of nodes, node pairs, double ribbons, and dipoles. 
Interactions between arrays of extended dislocations lead to the 
formation of networks of these configurations and to multipoles. 
Many of the observed configurations are similar to those reported 
in hexagonal graphite. The Class | configurations are similar to 
those found in low stacking fault energy fec metals. 
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Jennings, D. A., Power and energy measurement of repeti- 
tively pulsed lasers, Nat. Bur. Stand. (U.S.), Tech. Note 398, 
10 pages (March 1971) 25 cents, SD Catalog No. C13.46:398. 


Key words: Average power; energy per pulse; laser; peak power; 
repetitively pulsed lasers. 


The problem of measuring average power, energy per pulse, and 
peak power of some of the more common repetitively pulsed lasers 
is discussed. The techniques which have been used at the National 
Bureau of Standards are mentioned along with some of the ac- 
curacies obtained. Accuracies of 3 to 10 percent can be achieved 
depending on the laser source and the parameter of interest. 


Kamas, G., Hanson, D. W., Recent VHF/UHF satellite timing 
experiments at the National Bureau of Standards, Proc. 
Precise Time and Time Interval (PTTI) Strategic Planning Meeting 
(U), Washington, D.C., December 10-11, 1970, 1, 79-82 (1970). 


Key words: ATS satellite; LES—6 satellite; NBS satellite timing; 
satellites: TA CSAT satellite; time; timing accuracy. 


A summary discussion is presented of the current work accomplished 
by NBS in the field of time dissemination by geostationary satellites. 
The results of that program and its implication for the future are 
given. 


Kamper, R. A., Zimmerman, J. E., Noise thermometry with the 
Josephson effect, J. Appi. Phys. 42, No. 1, 132-136 (Jan. 1971). 


Key words: Josephson effect; noise; superconductivity; thermometry. 


Thermal noise causes a random frequency modulation of the self- 
oscillation of a Josephson junction, and the temperature of the 
noise source can be determined by analysis of the generated signal. 
We discuss the theoretical limitations of thermometry with this 
principle, and describe a prototype thermometer which has re- 
corded noise temperatures down to 0.075 K. 


Keysar, B. C., Specifications and tolerances for reference 
standards and field standard weights and measures. 2. 
Specifications and tolerances for field standard measuring 
flasks, Nat. Bur. Stand. (U.S.), Handb. 105-2, 6 pages (Jan. 1971) 
25 cents, SD Catalog No. C13.11:105-2. 


Key words: Accurate measurements of volume of liquids; field stand- 
ard measuring flasks; specifications; tolerances; weights and meas- 
ures inspection. 


These specifications and tolerances are recommended as minimum 
requirements for standards used in the field by State and local 
weights and measures officials in quantity determinations of liquid 
commodities. 


Keysar, B. C., Specifications and tolerances for Reference 
Standards and Field Standard Weights and Measures—3. 
Specifications and tolerances for Metal Volumetric Field 
Standards, Nat. Bur. Stand. (U.S.), Handb. 105-3, 8 pages (May 
1971) 25 cents, SD Catalog No. C13.11:105-3. 


Key words: Accurate measurement of volumes: field standard; metal 
volumetric field standards; provers: specifications; test measures; 
tolerances: weights and measures inspection. 


These specifications and tolerances are recommended as minimum 
requirements for standards used in the field by State and local 
weights and measures officials in quantity determinations of liquid 
commodities. 


Kruger, J., Frankenthal, R. P., B. Corrosion of metals, Chap- 
ter in Volume IV, Techniques in Metals Research, Chapter 10B 
oxidation and corrosion, R. A. Rapp, Ed., R. F. Bunshah, Series 
sy 14, Part 2, 571-667 John Wiley and Sons Inc., New York, 
.Y. 1970). 


Key words: Corrosion; electron diffraction: electron microscopy: 
research techniques. 


The application of non-electrochemical techniques to corrosion 
research are discussed. These include optical techniques, x-ray 
techniques, electron microscopy and diffraction, tracer techniques 
and mechanical techniques. 


Kulin, G. .urewitz, P. H., Editors, Hydraulie Research in the 
United States 1970--Including Contributions from Cana- 
dian Laboratories, Nat. Bur. Stand. (U.S.), Spec. Publ. 346, 
354 pages (Mar. 1971) $2.50, SD Catalog No. C13.10:346. 


Key words: Fluid mechanics; hydraulic engineering; hydraulic 
research; hydraulics; hydrodynamics; model studies; research 
summaries. 


Current and recently concluded research projects in hydraulics 
and hydrodynamics for the years 1969-1970 are summarized. 
Projects from more than 250 university, industrial, state and fed- 
eral government laboratories in the United States and Canada 
are reported. 


Kurylo, M. J., Peterson, N. C., Braun, W., Temperature and 
pressure effects in the addition of H atoms to propylene, 
J. Chem. Phys. 54, No. 11, 4662-4666 (June 1, 1971). 


Key words: Absolute rate constant; activation energy; fluorescence; 
H atom; propylene; reaction kinetics. 


The effect of He pressure on the rate of addition of hydrogen atoms 
to propylene was measured at 298 °K. H atom signal was recorded 
by the method of resonance of Lyman a@ radiation at 121.6 nm. 
The results are consistent with the 1/p's dependence of kp; recently 
observed in the H atom-ethylene system. The limiting high pressure 
rate constant obtained in this work is H+C3Hg— C3H; k:= 1.61 
+ 0.04 x 10>"? em’ mol c>! s*". 

The temperature dependence of the H atom addition was investi- 
gated at 50 torr total pressure over the temperature range 177-473 
°K. Analysis of the data below room temperature gives the following 
Arrhenius parameters for addition to the terminal olefinic position: 
K%= (10.18 + 0.26) x 10" exp [—(121] +11)/1.987] em? mol c™! s~'. 
Deviations from this dependence above 298 °K can be assigned to 
non-terminal addition, H atom abstraction, and increased decomposi- 
tion of the excited propyl radicals. 


Lamb, V. A., Johnson, C. E., Valentine, D. R., Physical and mech- 
anical properties of electrodeposited copper, III. Deposits 
from sulfate, fluoborate, pyrophosphate, cyanide, and 
amine baths. J. Electrochem. Soc. 117, Nos. 9, 10, 11, 291C—318C, 
341C-352C, 381C—404C (Sept., Oct., Nov. 1970). 


Key words: Annealing-electroylic copper: electrodeposited copper; 
impurity content; properties of mechanical structure. 


This paper is the third published report on the results of a broad 
program on properties of electrodeposited copper, sponsored 
jointly by the American Electroplaters’ Society, The International 
Copper Research Association, Inc., The Copper Development 
Association, and the National Bureau of Standards. The first paper 
is a comprehensive review of the literature to 1965. The second 
paper is an interim report on experimental results. The present 
and final paper, which incorporates the data in the second, includes 
data from the baths designated in the title, for deposits prepared 
under a wide range of operating conditions and with use of a variety 
of addition agents. Properties measured include tensile strength, 
yield strength, elongation, modulus of elasticity, fatigue limit, hard- 
ness, internal stress, density, electrical resistivity, and thermal 
expansivity. Properties were measured for as-plated deposits, de- 
posits after annealing at several temperatures, after cold-working, 
and at low and high ambient temperatures. Structure of deposits 
was examined by optical, electron micrograph, and x-ray methods. 
Content of impurities in deposits was determined. Correlations are 
developed among properties, structure, impurity content, type of 
bath, and operating conditions. 


Larsen, N. T., Clague, F. R., The NBS type II power measure- 
ment system, Proc. 25th Annual 1970 Instrument Society of 
America Conf. and Exhibit, Philadelphia, Pa., October 26-29, 
1970, 11 pages. 
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Key words: Bridge; leveler; microwave; power; self balancing; 
stabilizer. 


The design, construction, and performance of a new solid state 
microwave power measuring system are described. The instrument 
consists of two modules—a self-balancing bolometer bridge and a 
reference generator. The reference generator is a multi-function 
module which may be used as either a 0.0] percent differential 
voltmeter with a 0 to 10 volt range, a 0 to 10 volt reference voltage 
source of 0.01 percent absolute accuracy, or a precision microwave 
power leveler. A six decade divider provides high precision in all 
three functions. 


Laufer, A. H., Low-temperature chromatographic determina- 
tion of ketene and methyl ketene, /. Chromatogr. Sci. 8, 
677-678 (Nov. 1970). 


Key words: Analysis; gas chromatography; ketene; methyl ketene. 


This note describes the rapid, simultaneous analysis by conven- 
tional gas chromatography of both ketenes in small quantities. 


Liu, Y. M., Coleman, J. A., Electron radiation damage effects 
in silicon surface-barrier detectors, JEEE Trans. Nucl. Sci. 
NS-18, No. 1, 192-199 (February 1971). 


Key words: Alpha particles counting response: capacitance; detec- 
tor noise; electron fluence; electrons radiation damage: front and 
rear contact; leakage current; silicon surface-barrier detector. 


Silicon surface-barrier detectors have been irradiated at room 
temperature with monoenergetic electrons in the energy range of 
200 keV to 1 MeV. The changes of detector reverse leakage current, 
noise, capacitance and alpha-particle counting response were 
determined. In general, detector current and noise increased with 
electron fluence and energy for electron energies of 400 keV and 
above. Detector capacitances tended to decrease slightly for elec- 
tron fluences up to 10" cm~? and increase at higher fluences. No 
significant degradation of detector performance was observed as 
a result of irradiation with 200-keV electrons for fluences up to 
10'*/cm*. The effects of damage on detector performance were 
reduced when the rear, aluminum contact was irradiated rather 
than the front, gold contact. 


Lloyd, E. C., Editor, Accurate characterization of the high- 
pressure ervironment, Proceedings of a Symposium held at the 
National Bureau of Standards, Gaithersburg, Maryland, October 14- 
18, 1968, Nat. Bur. Stand. (U.S.,), Spec. Publ. 326, 343 pages (March 
1971) $4.50, SD Catalog No. C13.10:326. 


Key words: Accurate measurement; equation-of-state; fixed points: 
high pressure; high-pressure equipment; instrumentation; pressure 
scale; shock wave technique; temperature. 


The volume contains 38 papers prepared for the Symposium on Ac- 
curate Characterization of the High-Pressure Environment held on 
October 14-18, 1968, at Gaithersburg, Maryland, under the spon- 
sorship of the National Bureau of Standards and the Geophysical 
Laboratory of the Carnegie Institution of Washington. The papers 
are presented with the discussions that occurred during the sessions. 
The book also includes reports of several informal committees of 
the conferees on choices of reference pressure materials and on 
other matters relevant to improved measurement and calibration. 
The Symposium was intended to provide an authoritative survey 
of problems and techniques presently in use or proposed for precise 
high-pressure measurement and for temperature measurement at 
high pressure. 


Lutz, G. J., Phetoen activation analysis—A Review, Anal. 
Chem. 43, No. 1, 93-104 (January 1971). 


Key words: Bremsstrahlung; carbon; electron accelerators; nitrogen; 
oxygen; photon activation analysis; review: sensitivities. 


Prior to the availability of high energy electron accelerators, photon 
activation analysis was limited to the determination of deuterium 
and beryllium with isotope sources. During the past fifteen years, 





however, the development of high current cyclic and linear electron 
accelerators has made possible many analyses not readily accom. 
plished by other methods. These include the determination of beryl- 
lium and deuterium with isotope sources, accelerator produced 
photons and nuclear considerations in photon activation analysis, 
determination of carbon, nitrogen and oxygen, determination of the 
heavier elements with emphasis on biological, geochemical, oceano- 
graphic and forensic matrices, and errors and corrections in photon 
activation analysis. 


McCaa, W. D., Jr., Nahman, N. S., Generation of reference wave- 
forms by uniform lossy transmission lines, /EEE Trans. Instr. 
Meas. IM-19, No. 4, 382-390 (November 1970). 


Key words: Lossy uniform transmission lines; oscilloscope errors; 
pulse generators; standard waveforms. 


This paper describes a method for generating fractional nanosecond 
pulse waveforms of known shape thereby providing an a priori 
means for establishing reference waveform generators or standards 
to be used in pulse measurements and other applications. The 
method employs the bandlimiting properties of a lossy uniform trans- 
mission line to produce a known waveform and generator impedance. 
A theory for generator characterization and application is presented 
which defines a given generator in terms of time domain or frequency 
domain functions. A new term called the “available” waveform is 
defined. RMS error is employed as the error criteria. 

An example of generator design is given which employes planar skin 
effect loss and provides an output response having 0.2% RMS error 
over the interval 0 < ¢ <600 picoseconds referred to the inherent 
step response of the waveshaping line. The input transition genera- 
tor is assumed to be a 30 picosecond unit ramp generator. Also, an 
application example is given which employs a reference waveform 
generator to evaluate a sampling oscilloscope. The reference 
waveform generator uses a Debye Dielectric lossy transmission line. 


McConnell, P. M., Daney, D. E., Kirgis, J. B., Thermoelastic ex- 
pansion and creep of polyethylene terephthalate and poly- 
pyromelitimide film and polyethylene terephthalate fibers 
from 20 to 295 K, J. Appl. Phys. 41, No. 13,5066—5070 (Dec. 1970) 


Key words: Longitudinal; polyethylene terephthalate; polypyromelit- 
imide; relative creep; thermoelastic; transverse. 


A quartz tube dilatometer was used to measure the lineal thermal 
expansion and creep of single lengths of polyethylene terephthalate 
(PETP) film, Mylar, polypyromelitimide (PPMI) film, Kapton, and 
PETP multi-fiber yarn, Dacron, while stressed under constant 
tension. Tensions below and above the conventionally defined yield 
strength were used and the sample temperature ranged from 20 to 
295 K. Relative creep strain measurements, taken at the constant 
temperatures 77, 195 and 295 K, were found to obey the equation 


€= exp [—2.3 exp (A’y)] 


where y is a function of stress, time, and temperature and 4’ is a 
constant depending on the material. This equation was used to 
correct the thermoelastic expansion measurements for creep at 
the higher stresses. PETP multi-fiber yarn subjected to a slight 
tension was found to elongate during cooldown from 293 to 20 K. 
Higher stresses caused less elongation; i.e., the coefficient of ex- 
pansion increased with stress. This result is believed to be due to 
changes in crystallanity at the higher stresses. A similar stress effect 
was found with PETP film but not with PPMI film. The thermoelastic 
expansivity of the film samples was also found to be sensitive to the 
thickness. 


McDonald, D. G., Evenson, K. M., Wells, J. S., Cupp, J. D., High- 


frequency limit of the Josephson effect, J. Appl. Phys. 42, 
No. 1, 179-181 (Jan. 1971). 


Key words: Detectors; infrared; Josephson effect; lasers; super- 
conductivity. 


Coherent detection of applied radiation via a Josephson current-step 
has been achieved with an applied frequency approximately a 
factor of 3 higher than reported heretofore. The detected signal is at 
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2.5 THz (118 ym). The inadequacy of existing theory to predict 
the frequency dependence of the Josephson effect is discussed. 


McLaughlin, W. L., Hussmann, E. K., Eisenlohr, H. H., Chalkley, 
L.. A chemical dosimeter for monitoring gamma-radiation 
doses of 1-100 krad, Intern. J. Appl. Radiation Isotopes 22, 
135-140 (1971). 


Key words: Chemical dosimeter; disinfestation; dosimeter; dyes; 
gamma rays; insect sterilization; shelf-life extension: sprouting 
inhibition. 


A simple chemical dosimeter is described for measuring gamma-ray 
doses useful for insect sterilization, seed-sprouting inhibition, and 
food shelf-lifs extension. The solutions, colorless before irradiation, 
assume a stable blue-violet color when irradiated to absorbed doses 
from 1 to 100 kilorads. The readout may be made either visually, 
colorimetrically, or spectrophotometrically. The optical density is 
linear with dose, and the response does not vary with dose rate. 


Mackay, D. R., Technical Standards Coordinator, Editorial 
format for Voluntary Product Standards. Part I. Contents 
of a Voluntary Product Standard, Nat. Bur. Stand. (U.S.), Prod. 
Stand. 0-70, 13 pages (March 1971) 15 cents, SD Catalog No. 
C13.20/2:0—70. Supersedes PS 0-67. 


Key words: Editorial format; format; product standards; standard; 
voluntary product standards. 


This publication is intended to assist individuals, organizations, and 
agencies which are concerned with the writing of standards. It 
establishes and describes the material to be contained in a standard 
and provides specific instructions for writing standards. It is pri- 
marily intended for those developing initial drafts of standards to 
be submitted to the Bureau under the Procedures for the Development 
of Voluntary Product Standards published by the Department of 
Commerce. 


Mandel, J., A new analysis of variance model for non-additive 
data- Technumetrics 13, No. 1, 1-18 (February 1971). 


Key words: Factorial experiments; interaction; non-additivity; 
principal components; surface fitting; two-way arrays. 

A method is presented for the analysis of data representing functions 
of two variables, when the response can be tabulated in a rectangular 
array. The procedure is based on a partitioning of the row by column 
interaction effects into a sum of terms, each of which is the product 
of a row factor by a column factor. The factors in each term are 
estimated by a method involving the extraction of characteristic 
roots. 

The method contains as special cases a number of procedures used 
for the handling of non-additivity in two-way arrays. It is very useful 
for the fitting of empirical surfaces, but is also applicable to cases 
in which the data depend on qualitative rather than quantitative 
factors. 

Comparisons with other techniques are made and an illustrative 
example is given. 


Mangum, B. W., Thornton, D. D., Vibrating sample magnetom- 
eter for use at very low temperatures and in high magnetic 
fields, Rev. Sci. Instr. 41, No. 12, 1764-1766 (December 1970). 


Key words: Bimorph; cryogenic; low-temperature; magnetometer; 
piezoelectric; vibrating-sample magnetometer. 


A high-sensitivity vibrating-sample magnetometer, for use at the 
very low temperatures produced by a *He refrigerator or by a *He- 
‘He dilution refrigerator is described. It can be made extremely 
compact by the use of a piezoelectric element as the vibrator. 
The complete magnetometer can be placed in the cryostat between 
the poles of an electromagnet, and can be operated at high applied 
magnetic fields. 


Manning, J. R., Correlation effects and activation energies for 
diffusion in alloys (Proc. Conf. Atomic Transport in Solids and 
Liquids, Marstand, Sweden, June 1970), Z. Naturforsch. 26A, No. 1, 
69-76 (Jan. 1971). 


Key words: Activation analysis; concentrated alloys; correlation 
factor; diffusion; random alloy; vacancies. 


The problems involved in calculating correlation factors for diffusion 
in dilute alloys can be contrasted to those arising in concentrated 
solid solutions. As one moves from the pure element to the dilute 
alloy to the concentrated alloy, the calculation becomes progres- 
sively more difficult. Because of the complex atom configurations 
which can occur’in concentrated alloys, it usually is not possible to 
calculate correlation factors in these alloys exactly. 

Several important simplifications are available in nondilute random 
alloys. A large reduction in complexity can be secured by using a 
random alloy model where each atom is treated as diffusing in a 
uniform matrix, with the matrix properties being determined by 
the composition and jump frequencies in the alloy. Resulting equa- 
tions in this random alloy model can be expressed directly in terms 
of the experimentally measurable tracer diffusion coefficients with 
no unknown vacancy jump frequencies appearing. Also these equa- 
tions have the advantage of being in simple analytic form and not 
requiring numerical methods to evaluate the correlation factors. 
These two features make possible the direct expression of the 
temperature dependence of the correlation factor in terms of the 
experimental activation energies. 

Equations are found for AH/AQ in random binary cubic alloys, 
where AH is the difference between the activation enthalpies for 
diffusion of the two species and AQ is the difference between the 
experimentally measured activation energies of the two species. 
This ratio is never less than unity and can be much larger than unity. 
Values are plotted for diamond, body-centered cubic and face- 
centered cubic structures. From the magnitude and composition 
dependence of AH/AQ, it is concluded that the temperature de- 
pendence of the correlation factor cannot by itself explain the dif- 
ference between the activation energies measured from tracer 
diffusion and from internal friction in the non-dilute range. 


Marsden, C. P., Editor, Silicon device processing. Proceedings 
of a Symposium held at Gaithersburg, Maryland, June 2-3, 1970, 
Nat. Bur. Stand. (U.S.), Spec. Publ. 337, 467 pages (Nov. 1970) 
$5.50, SD Catalog No. C13.10:337. 


Key words: Analysis; device processing; diffusion; epitaxy; junctions; 
resistivity; surface preparation. 


This Proceedings contains the information presented at the Sym- 
posium on Silicon Device Processing held at the National Bureau 
of Standards on June 2-3, 1970. 

This Symposium covered the measurement field and was purpose- 
fully restricted to the generic topics of diffusion, epitaxy, surface 
preparation and interdependence of unit processing operations. 
This emphasis on measurement during the processing operation or 
on the characteristics of the processed material, showed the neces- 
sity and more important, the ambiguities of current methods of 
measurement. The application of some new techniques to measure- 
ment were also discussed. 


Marshak, H., Dove, R. B., A right angle *He cryostat incorporat- 
ing a high field superconducting solenoid, Nat. Bur. Stand. 
(U.S.), Tech. Note. 562, 22 pages (Dec. 1970) 30 cents, SD Catalog 
No. C13.46:562. 


Key words: *He cryostat; *He refrigerator; nuclear orientation: 
superconducting solenoid. 


Construction and operation of a novel *He cryostat, incorporating a 
large superconducting solenoid mounted at right angles to the 
cryostat’s vertical axis, is described. This new cryostat which is part 
of the transportable National Bureau of Standards *He refrigerator, 
has been used successfully for nuclear orientation studies at the 
Atomic Energy Research Establishment, Harwell, England. 


Martin, J. F., Standard Reference Materials: National Bureau 
of Standards-U.S. Steel corporation joint program for de- 
termining oxygen and nitrogen in steel, Nat. Bur. Stand. (U.S.), 
Spec. Publ. 260-26, 40 pages (February 1971) 50 cents, SD Catalog 
No. C13.10:260-26. 
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Key words: Inert gas fusion; one ppm of oxygen and nitrogen in 
steel; research associate program; simultaneous determination of 
oxygen and nitrogen; standard reference material for oxygen. 


Because of a need for a method for the rapid, simultaneous determi- 
nation of small concentrations of oxygen and nitrogen in steel, a 
joint project under the Industrial Fellowship Program of U.S. Steel 
and the Research Associate Program of the National Bureau of Stand- 
ards (NBS) was established. After investigation of various techniques, 
an analytical apparatus was constructed with excellent sensitivity 
for both elements. This apparatus consists of an inert gas fusion sys- 
tem coupled to a gas chromatograph. Recommended modifications 
have been completed on this equipment, and simultaneous deter- 
minations of oxygen and nitrogen in steel are presently being made. 
The detection limit is <1 ppm for both oxygen and nitrogen, and the 
analysis time is about 5 minutes. 

A direct result of this program has been the issuance by NBS of a new 
Standard Reference Material (SRM), a maraging steel, for oxygen. 
This SRM has the lowest oxygen concentration (4.5 ppm) of any stee! 
standard yet issued by NBS. 


Mason, H. L., Editor, Innovative metrology—key to progress. 
Proceedings of the 1970 Standards Laboratory Conference, Nat. Bur. 
Stand. (U.S.), Spec. Publ. 335, 132 pages (Mar. 1971) $1.50, SD Cata 
log No. C13.10:335. 


Key words: Metrology management; National Conference of Stand- 
ards Laboratories; physical measurement; proceedings NCSL. 


The biennial Standards Laboratory Conference of the National Con- 
ference of Standards Laboratories convened at the Gaithersburg 
facilities of the National Bureau of Standards June 15-17, 1970. The 
theme of the meeting, Innovative Metrology— Key to Progress, was 
amplified by 23 papers presented at technical sessions devoted to 
new technologies and applications, laboratory management and op- 
erations, new methods of optimizing calibration intervals, new ways 
of managing, and new international developments. 


Meinke, W. W., Standard reference materials for clinical 
measurements, Anal. Chem. 43, 28A-47A (May 1971). 


Key words: Clinical chemistry; organic SRMs; standard reference 
materials; spectrophotometry SRMs. 


The NBS Program in Standard Reference Materials for Clinical 
Chemistry measurements is described. Each of the SRMs issued 
in this area in the last five years is discussed. Future directions 
of the program are also mentioned. 


Menis, O., Editor, Status of thermal analysis. Proceedings of a 
Symposium on the Current Status of Thermal Analysis Held at 
Gaithersburg, Maryland, April 21-22, 1970, Nat. Bur. Stand. 
(U.S.), Spec. Publ. 338, 189 pages (Oct. 1970) $1.00, SD Catalog 
No. C13.10:338. 


Key words: Atherosclerotic plaque; DSC; elastomers and vulcani- 
zates; explosives; high temperature DTA; kinetics by TGA; selec- 
tion of DTA parameters; temperature standards; TGA; theory of 
isoperibol and adiabatic shield calorimeters. 


The symposium papers offer contributions in differential thermal 
analysis (DTA), differential scanning calorimetry (DSC), develop- 
ment of standards for DTA temperature. scales, and applications 
in high temperature biochemical, polymer and explosive materials. 
The selection and the effect of experimental parameters on the 
types of information on thermal curves is described. The measure- 
ment theory of isoperibol and adiabatic shield calorimeters by the 
method of intermittent heating is described and the magnitude of 
calorimetric instrumental errors is outlined. In an analysis of a 
differential scanning calorimeter, three instrumental time con- 
stants are described and the necessary corrections for the various 
instrumental and thermal time constants are recommended. The 
need for and the status of temperature scale standards for DTA 
is discussed. An evaluation of 12 materials for use as standards 
by cooperating laboratories is presented. High temperature (> 1900 
K) DTA, problems of high temperature calorimetry to 1300 K are 
discussed. The implications concerning the nature of biological 


membranes and decomposition of atherosclerotic plaques are 
derived from differential scanning calorimetry of three-component 
systems of phospholipid, cholesterol or one of its esters and water. 
A description is given of an apparatus developed for measuring 
the rate at which vapors are evolved during the thermal degradation 
of material in a modified thermogravimetric apparatus. The modifi- 
cations of a DTA cell to minimize explosions is described. Data 
derived from the use of lead azide was used to evaluate the tech- 
nique. Thermogravimetric analysis for establishing basic composi- 
tion of elastomer compounds and vulcanizates were reviewed. A 
careful study of oxidation characteristics of carbon black in the 
formulation is also included. 


Menis, O., Rains, T. C., Sensitivity, detection limit, precision 
and accuracy in flame emission and atomic absorption 
spectrometry, Chapter 2, Analytical Flame Spectroscopy, 
Selected Topics, R. Mavrodineanu, Ed., pp. 47-77 (Springer- 
Verlag Publishing Co., New York, N.Y., 1970). 


Key words: Accuracy; atomic absorption; detection limits; differen- 
tial method sensitivity; flame emission; nonaqueous media; preci- 
sion; tables. 


This chapter on sensitivity, detection limits, precision and accuracy 
in Flame Emission and Atomic Absorption spectrometry discusses 
the principles, definitions and methodology in establishing and eval- 
uating these criteria. The various instrumental components and 
their optimum operation are described to enable the analyst to 
achieve these goals. Tables of data for all reported elements are 
presented for the current results on detection limits and sensitivity, 
The problem associated with the environment of the analyte and 
the effect of nonaqueous media are discussed from the standpoint 
of enhanced sensitivity. Finally the criteria for attaining improved 
precision and accuracy are discussed and the role of the differential 
method, ionization buffers and standard additions is discussed. 


Moore, C. E., Selected tables of atomic spectra. A Atomic 
energy levels-second edition. B Multiplet tables. C 1, C7, C 
m, Civ,C Vv, C vi. Data derived from the analyses of optical 
spectra, Nat. Stand. Ref. Data Ser., Nat. Bur. Stand. (U.S.), 3, 
Sec. 3, 71 pages (Nov. 1970) $1.00, SD Catalog No. C 13.48:3/sec. 3. 


Key words: Atomic energy levels, carbon spectra; atomic spectra 
of carbon; carbon spectra; multiplet tables, carbon spectra; spectra, 
carbon; wavelengths, carbon spectra. 


The present publication is the third Section of a series being pre- 
pared in response to the persistent need for a current revision of two 
sets of tables containing data on atomic spectra as derived from 
analyses of optical spectra. As in the the first two Sections, Part A 
contains the atomic energy levels and Part B the multiplet tables. 
All six spectra of carbon, C 1 through C vi are included. The form 
of presentation is described in detail in the text to Section 1, and need 
not be repeated here. 


Motz, J. W., Dick, C. E., Lucas, A. C., Placious, R. C., Sparrow, 
J. H., Production of high-intensity K x-ray beams, J. Appl. 
Phys. 42, No. 5, 2131-2133 (April 1971). 


Key words: Beryllium to gold; electron produced; high intensity; 
monochromatic yields; purity; x-ray beams. 


K x-ray beams are produced by direct electron excitation in beryl- 
lium, carbon, aluminum, titanium, copper, silver, and gold targets. 
With excitation currents of 10 milliamperes from d.c. electron 
accelerators or 1000 amperes from single pulse machines, the K 
x-ray yields are respectively 4 to 9 orders of magnitude larger than 
the yields produced by the conventional fluorescence method. The 
bremsstrahlung contamination in the K x-ray beam is minimized 
by observing the K x-rays at 180 degrees with respect to the incident 
electron beam, in order to exploit the fact that the bremsstrahlung 
intensities become more peaked in the forward direction as the 
electron energy increases. Values are given for the maximum yields 
and purities of the K x-ray beams with the corresponding electron 
energies, and a comparison is made with values obtained by other 
methods involving atomic ion or photon excitation of K x-rays. 
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Mukerjee, P., Mysels, K. J., Critical micelle concentrations of 
aqueous surfactant systems, Nat. Stand. Ref. Data Ser., Nat. 
Bur. Stand. (US.), 36, 227 pages (Feb. 1971) $3.75, SD Catalog 
No. C13.48:36. 


Key words: Association colloid; bibliography; CMC; colloid; colloidal 
electrolyte; critical concentration; critical micelle concentration; 
detergent; hydrophobic bonding; Krafft point; long chain compounds; 
micelle; paraffin chain salts; selected values; soap; solubilization; 
standard values; surface active agents; surface chemistry; surface 
tension; surfactant. 


Critical micelle concentrations (CMC’s), have been collected, orga- 
nized and evaluated. The literature has been scanned for numerical 
values from 1926 up to and including 1966. In addition, over 800 
values, hitherto available only in graphical form or implied in experi- 
mental data, have been extracted from the publications and are 
included. Close to 5,000 entries, based on 333 references, dealing 
with 720 compounds are tabulated in the main tables. Whenever 
available, the temperature, any additives present, the method of 
determination and the literature source are given for each CMC 
value and an indication of the apparent quality of the preparation 
and method used are included. A shorter table gives selected values 
which are believed to be particularly reliable, including highly 
accurate ones. Among these, concordant values from at least two 
independent laboratories are emphasized. 

Included in the Introduction is a general discussion of the importance 
and significance of CMC values and of methods for their determina- 
tion, as well as a summary of the procedures used in the collection, 
evaluation and presentation of these values in the present work. 
Extensive indexes are provided. 


Murphey, W. M., Caswell, R. S., Analysis of results of the Bureau 
International des Poids et Mesures thermal neutron flux 
density intercomparison, Metrologia 6, No. 4, 111-115 (October 
1970). 


Key words; Data analysis; flux density; intercomparisons; least 
squares; standards; thermal neutrons. 


The results of an international comparison of thermal neutron 
flux density standards sponsored by the Neutron Working Group of 
the Bureau International des Poids et Mesures have been analyzed 
using a method conceptually different from a companion analysis 
by Axton. Both methods, however, yield closely the same results. 
All individual values are within +2.6% of the adjusted value, the 
average absolute deviation being about 1%. 


Napolitano, A., Hawkins, E. G., Standard Reference Materials: 
Viscosity of a standard borosilicate glass (Certification of 
Standard Reference Material 717), Nat. Bur. Stand. (U.S.), 
Spec. Publ. 260-23, 10 pages (Dec. 1970) 25 cents, SD Catalog No. 
C13.10:260-23. 


Key words: Beam-bending; borosilicate glass; glass viscosity; paral- 
lel-plate; rotating cylinder; standard reference material; viscosity; 
viscosity standard. 


The viscosity of a borosilicate glass has been measured at the Na- 
tional Bureau of Standards and four other laboratories. Determina- 
tions were made in the range 10 to 105 poises (1525 to 470 °C). 
Measurements were made by the rotating cylinder, fiber elongation, 
beam bending, and parallel-plate methods. The results have been 
evaluated and the glass has been issued as Standard Reference 
Material No. 717. 


Neill, A. H., Jr., Wet traction of tractionized tires, Nat. Bur. 
Stand. (U.S.) Tech. Note 566, 14 pages (Feb. 1971) 20 cents, SD 
Catalog No. C13.46:566. 


Key words: Siping; stability; stopping distance tires; traction; 
tractionizing. 


A series of dynamic vehicle tests was performed at NBS to eval- 
uate the performance of tractionized or siped tires. Stopping dis- 
tance and lateral breakaway data is presented from a two wheel 
diagonally braked automobile which clearly shows that siped tires 


do not represent any improvement in the lateral stability or stopping 
distance characteristics of a typical passenger automobile. 


Nimeroff, I., Color-match classification by variable param- 
eters, Color. Engr. 9, No. 2, 13-17 (March-April 1971). 
color; color-match; isomerism; 


Key words: Classification; 


metamerism. 


There is an increasing need for standard procedures by which 
to assess the nature of achieved color matches. It has been found 
that there are seven different classes of color matches in which 
the colorimetric parameters, the spectral characteristics of ob- 
server, source and specimen, are variable. The different classes 
result from the seven possible combinations of these three vari- 
able parameters, taken one, two and three at a time. Those matches 
in which one parameter varies are called, in accord with the vari- 
able parameter, observer-variable, source-variable or specimen- 
variable matches; those in which two parameters vary are called 
source-specimen-variable, source-observer-variable, or specimen- 
observer-variable matches; those in which all three vary are called 
source-specimen-observer-variable matches. Three of these classes 
of color match are currently designated metameric matches, one 
class is currently designated spectral match, while the other three 
classes have no current designation. The colorimetric effect of spec- 
tral differences in all of these classes of color matches may be charac- 
developed by Nimeroff and Yurow for matches that are currently 
designated metameric matches. A generalization of this index and 
its applications are presented. 


Ogburn, F., Shives, T. R., Brady, C. H., Calibration of a filar 
micrometer, J. Mater. 6, No. 1, 60-66 (March 1971). 


Key words: Metal coatings; micrometers; microscopes; optical 
equipment; optical microscopes; thickness; thickness gages. 


The precision with which a filar micrometer can be calibrated against 
a fixed interval of a stage micrometer was determined at magnifica- 
tions of about 175 and 1900 for a group of 15 individuals familiar 
with microscopes. At each magnification the spread of measurements 
among the 15 operators was about 0.8 percent, whereas the spreads 
for the individual operators were substantially less. 


Osinsky, V. I., Winogradoff, N. N., Excitation and temperature 
dependence of band-edge photolumi ence in gallium 
arsenide, Phys. Rev. B. 3, No. 10, 3341-3346 (May 15, 1971). 





Key words: Gallium arsenide lasers; impurity effects; intensity 
dependence; radiative recombination. 


Photoluminescent spectra of n-type GaAs were studied as a function 
of the excitation intensity, temperature, and doping level. The 
spectra consisted of two major bands representing radiative band-to- 
band recombination and radiative transitions through impurity 
centers, respectively. The intensity of the peak of the former 
went through a minimum and a maximum as the temperature was 
increased from 175 to 500 °K. The temperature corresponding to 
the above maximum increased as the excitation intensity was de- 
creased or the doping level increased. These results suggest that the 
temperature dependence of the peak intensity in band-to-band 
transitions is primarily due to the thermal distribution of the carriers 
over the available energy states. The peak intensity would therefore 
normally be expected to decrease monotonically with an increase 
in temperature, while the above maxima and minima represent per- 
turbations imposed by the presence of temperature dependent tran- 
sitions through radiative or nonradiative impurity recombination 
centers. 


Page, C. H., Comments on “Nonlinear resistors that generate 
subharmonics,”’ Proc. [IEEE Letters 59, No. 2, 301 (Feb. 1971). 


Key words: Nonlinear; subharmonics. 
The analysis reported by Erdey is shown to be incorrect. 
Page, C. H., Vigoureux, P., Editors, The international system of 


units (SI), Nat. Bur. Stand. (U.S.), Spec. Publ. 330, 42 pages (Jan. 
1971), 50 cents, SD Catalog No. C13.10:330. 
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Key words: General conference of weights and measures; inter- 
national committee of wéights and measures; international system 
of units; SI; systémé international des unités; units of measure- 
ment. 


This translation from the French “Le Systémé International d‘ 
Unités”, (SI) published originally by the International Bureau of 
Weights and Measures (BIPM) has been prepared jointly by the 
National Physical Laboratory, UK, and the National Bureau of 
Standards, USA. Included are Resolutions and Recommendations 
of the General Conference of Weights and Measures (CGPM) on 
the International System of Units, together with relevant extracts 
from the International Organization for Standardization (ISO) 
for the practical use of the system. 

Appendix I gives in chronological order the decisions promulgated 
since 1889 by CGPM and the International Committee of Weights 
and Measures (CIPM) on units of measurement and on SI. Appendix 
II outlines the measurements, consistent with the theoretical 
definitions given in this document, which metrological laboratories 
can make to realize the units and to calibrate precision material 
standards. 


Parks, E. J., Thermal analysis of modified cellulose, Tappi 
54, No. 4, 537-544 (April 1971). 


Key words: Alum; aluminum; apparent activation energy; calcium; 
cellulose; crosslinking; differential; ion exchange; organic acid; oxy- 
cellulose; stability; sulfuric acid; thermal analysis; thermogravi- 
metric analysis. 


Differential thermal analysis (DTA) and thermogravimetric analysis 
(TGA) techniques have been applied to characterize thermal proper- 
ties of cellulose and modifications of the type that might reasonably 
occur in papermaking pulps. Periodate, periodite-chlorite, chromate, 
chromate-chlorite, hypochlorous acid, and hypochlorous acid- 
chlorite oxycelluloses all were examined, as were their calcium and 
aluminum exchange products. In addition, the effects on thermal 
properties of cellulose exerted by alum, sulfuric acid, and two organic 
acids that cause color reversion all were investigated. The peak 
temperature of a massive decomposition endotherm (T3), the 
apparent activation energy of this endothermic reaction (E *) and 
the integral procedural decomposition temperature (ipdt) seem to 
be best suited of thermal variables to denote stability. Oxidation 
usually destabilizes cellulose slightly. An exception is exhaustive 
oxidation of periodate oxycellulose with chlorous acid. This causes 
T; and E* to increase, although the ipdt drops. Exchange with 
calcium greatly stabilizes cellulose containing carboxyl groups, 
whereas aluminum destabilizes the same material. Each effect is 
concentration dependent. The importance of pH in these relation- 
ships is not known. Sulfuric acid has a destabilizing effect qual- 
itatively different from that of aluminum covering carboxyls and 
that of alum added to cellulose. Destabilization by alum is probably 
a result of contributions both by sulfate and by aluminum. Organic 
acids have little effect. 

Each treatment has resulted in a higher percentage of residue than 
that of cellulose, which indicates the promotion of crosslinking 
by oxygen, metals, and acid. 


Paule, R. C., Mande!, J., Analysis of interlaboratory measure- 
ments on the vapor pressures of cadmium and silver (Certifi- 
cation of Standard Reference Materials 746 and 748), Nat. 
Bur. Stand. (U.S.), Spec. Publ. 260-21, 30 pages (Jan. 1971), 35 cents, 
SD Catalog No. C13.10:260-21. 


Key words: Cadmium; components of error (within- and between- 
laboratory); heats of sublimation (second and third law); interlabora- 
tory measurements; silver; standard errors; standard reference 
materials; vapor pressure. 


Detailed statistical analyses have been made of results obtained from 
a series of interlaboratory measurements on the vapor pressures of 
cadmium and silver. Standard Reference Materials 746 (cadmium) 
and 748 (silver) which were used for the measurements have been 
certified over the respective pressure ranges 10-''—10-* atm and 
10-"2—10-3 atm. The temperature ranges corresponding to these 
pressures are 350-594 K for cadmium and 800-1600 K for silver. The 
heats of sublimation at 298 K and the associated two standard error 
limits for cadmium and silver are 26660 + 150 cal/mol and 68010 + 300 


cal/mol, respectively. Estimates of uncertainty have been calculated 
for the certified temperature-pressure values as well as for the un- 
certainties expected from a typical single laboratory’s measure. 
ments. The statistical analysis has also been made for both the sec. 
ond and third law methods, and for the within- and between-labora- 
tory components of error. The uncertainty limits are observed as 
functions of both the heat of sublimation and the temperature. 


Peavy, S. T., Varner, R. N., Hogben, D., Source listing of 
OMNITAB II program, Nat. Bur. Stand. (U.S.), Spec. Publ. 339, 
371 pages (Dec. 1970) $4.75, SD Catalog No. C13.10:339. 


Key words: Accuracy; algorithms; ANSI] FORTRAN; documenta- 
tion; error checking; machine independent; OMNITAB II operating 
system subprograms; OMNITAB II source listing; programming 
techniques; transportable; user-oriented computing system. 


OMNITAB II is a general-purpose interpretive computing system 
designed to allow a nonprogrammer to use a high-speed computer 
easily, accurately and effectively. The system permits the user to per- 
form arithmetic operations including complex arithmetic, trigono- 
metric calculations, miscellaneous function calculations, statistical 
analysis, Bessel function calculations, and operations on matrices 
and arrays. 

The OMNITAB II system contains 177 subprograms written in the 
ANSI FORTRAN language. Every effort has been made to make the 
system transportable. This publication contains a complete listing 
of all these subprograms. The listing is preceded by a brief introduc- 
tion which describes the programming techniques used; the use of 
system library functions; and the subprograms used to control the 
flow of operations in the OMNITAB system. 


Phucas, C. B., Technical Standards Coordinator, Package 
quantities of green olives, Nat. Bur. Stand. (U.S.), Prod. Stand. 
40-70, 9 pages (May 1971) 10 cents, SD Catalog No. C13.20/2:40-70. 


Key words: Green olives, package quantities of; olives, green, 
package quantities of; package quantities of green olives. 


This Voluntary Product Standard covers a range of package quan- 
tities that are recommended for green olives and establishes specific 
packaging requirements in terms of net drained weight. Methods of 
labeling products which comply with this Standard are provided. 


Phucas, C. B., Technical Standards Coordinator, Package 
quantities of instant mashed potatoes, Nat. Bur. Stand. (U.S.), 
Prod. Stand. 41-70, 8 pages (April 1971) 10 cents, SD Catalog 
No. C13.20/2:41-70. 


Key words: Instant mashed potatoes, package quantities of; mashed 
potatoes, instant, package quantities of; package quantities of 
instant mashed potatoes: potatoes, instant mashed, package quan- 
tities of. 


This Voluntary Product Standard covers a range of package quan- 
tities based on servings and establishes the definition of a serving 
which is based on the weight of the reconstituted product. 


Phucas, C. B., Technical Standards Coordinator, Package 
quantities of instant nonfat dry milk, Vat. Bur. Stand. (U.S.), 
Prod. Stand. 37-70, 9 pages (April 1971) 10 cents, SD Catalog 
No. C13.20/2:37-70. 


Key words: Dry milk package quantities; instant nonfat dry milk; 
nonfat dry milk, package quantities of; package quantities of 
dry milk. 


This Voluntary Product Standard covers package quantities of 
instant nonfat dry milk based on quart equivalents within the 
range of 3 to 20 quarts. It also establishes the weight of instant 
nonfat dry milk required to make | quart. This Standard covers 
all types of packages of instant nonfat dry milk intended for house- 
hold use including those containing individual quart or multi 
quart equivalent envelopes. 
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Potzik, J.. A Method for determining the dynamic response 
of an elastic load cell element, Rew. Sci. Instr. 41, No. 12, 1726- 
1731 (December 1970). 


Key words: Dynamic force; load cell; viscoelastic damping. 


A method has been developed for determining the dynamic charac- 
teristics of a load cell element by subjecting it to sinusoidally 
time-dependent forces, and an experiment using a piezoelectric 
shaker is described which illustrates the feasibility of the method. 
At present a simple viscoelastic model is assumed for the load cell 
element, and relaxation times on the order of 10-7 seconds have 
been measured. The application of this information to load cells 
subjected to other types of dynamic forces is also discussed. 


Preston, J. D., Forthofer, R. J., Correlation of vehicle, dynamom- 
eter and other laboratory tests for brake friction materials, 
Proc. Society of Automotive Engineers Conf., Detroit, Michigan, 
January 11-15, 1971. Paper No. 710250, pp. 1-10 (1971). 


Key words: Brake friction materials; coefficient of friction; friction 
assessment screening test machine; friction materials test machine, 
girling scale dynamometer; inertia brake dynamometer; vehicle 
road tests. 


The frictional properties of brake lining materials are greatly 
influenced by the manner in which they are used. Test procedures 
designed to measure the coefficient of friction of these materials 
can likewise produce widely differing results depending on the type 
of test conducted. This paper presents data illustrating the per- 
formance correlation of three commercially available brake lining 
materials when subjected to vehicle tests and some of the more 
commonly used laboratory test procedures. 


Radlinski, R. W., Forthofer, R. J., Harvey, J. L., Operating per- 
formance of motor vehicle braking systems as affected by 
fluid water content, Proc. Society of Automotive Engineering 
Conf., Paper No. 710253, pp. 1-11 (January 11-15, 1971). 


Key words: Automotive hydraulic brake fluids; brake fluids per- 
formance; braking systems; low temperature viscosity. 


Automotive hydraulic brake fluids are known to pick up water in 
service. This paper shows that a major portion of the water picked 
up by the brake fluid in a braking system is transmitted through the 
hoses. It discusses two major effects of water pick-up on brake fluids 
performance —reduction in vapor lock temperature and increase in 
low temperature viscosity. Performance data for several ‘“‘wet” 
brake fluids operating in actual braking systems at high and low 
temperatures are presented. 


Ramboz, J. D., A link-compensated ratio transformer bridge, 
Rev. Sci. Instr. 42, No. 4, 522-524 (April 1971). 


Key words: Bridge; impedance; impedance bridge; measurement; 
ratio transformer; ratio transformer bridge. 


A ratio transformer impedance bridge that is capable of measuring 
either capacitive or inductive impedances is described. It is made of 
commercially available laboratory apparatus and can yield measure- 
ments to an uncertainty of less than +0.1 percent. The unique 
feature is an injection transformer which reduces the effect of the 
interconnecting “link-resistance” between the unknown impedance 
and the bridge reference reactance. Five circuit variations of the 
bridge are shown and the equations which yield the solutions for 
the unknown resistance, capacitance, or inductance and phase 
angle are given. 


Reitwiesner, G. W., On computer performance measurement 
programming measuring indexing adroitness by isolating 
complex primes, Vat. Bur. Stand. (U.S.), Tech. Note 572, 25 pages 
(Apr. 1971) 35 cents, SD Catalog No. C13.46:572. 


Key words: Assessment; complex; composite; computer; criteria; 
evaluation; Gaussian primes; indexing; measurement; performance; 
prime; program; test. 


This writing, describing a computer performance test program, is 
concerned not primarily with specific measurements, but rather 
with a procedure for making measurement regarding specific prop- 
erties of computer operation. 

The program is written in a particular problem-oriented program- 
ming language; therefore assessment perforce spans the effects of 
the computer hardware, of the programming language, and of the 
intervening compiler processes. 

The objective of the test is to assess adroitness in certain indexing 
operations. Assessment is accomplished by measuring execution 
time of a recursive programming loop. 

The test problem was chosen as a convenient artifice to use certain 
specific indexing-type operations in the programming employed for 
solution. 

The test program performs a simple computation for which the solu- 
tion is completely definitive, yet for which both the solution and the 
time for achieving it are variable under parameters whose values 
are introduced as program input data. 


Richmond, J. C., Editor, Space simulation. Proceedings of a 
Conference held at NBS Gaithersburg, Maryland, September 14-16, 
1970. Nat. Bur. Stand. (U.S.), Spec. Publ. 336, 984 pages (October 
1970), $5.25, SD Catalog No. C13.10:336. 


Key words: Ablation; particulate radiation; re-entry; space simula- 
tion; thermal radiation; vacuum; weightlessness. 


This volume contains all of the papers presented at the Fifth Space 
Simulation Conference held at the National Bureau of Standards 
September 14-16, 1970, that were available for publication. The 
general scope of the conference was the effect of the space environ- 
ment on materials, components, structures and man. The range of 
topics is too wide to permit a simple classification, but contamina- 
tion, ablation, degradation of materials by the space environment 
and predictive testing account for approximately half of the papers 
presented. Other topics range from purely laboratory problems such 
as radiometry standards, calibration of vacuum gages, cryopumping 
and operation of space simulation facilities to gravity simulation, 
with neutral buoyancy for zero gravity and a man-rated centrifuge 
for high gravity, use of a drifting submarine to study the psychological 
aspects of long-duration missions in a space station, and simulation 
of (1) atmospheric balloon environments, (2) radiation from nuclear 
power sources, (3) solar wind, (4) micro-meteoroid bombardment, 
(5) soil friction on the moon and (6) the Martian atmosphere. 


Risley, A. S., The physical basis of Atomic Frequency Stand- 
ards, Nat. Bur. Stand. (U.S.), Tech. Note 399, 54 pages (April 1971) 
60 cents, SD Catalog No. C13.46:399. 


Key words: Atomic frequency standards; cesium beam; energy 
levels; hydrogen maser; hyperfine interaction; rubidium gas cell; 
transition probability. 


A tutorial discussion of the physical basis of atomic frequency 
standards is given. These principles are then related to the condi- 
tions under which an atom can be used as the working substance of 
a stable and accurate frequency standard. The three primary ex- 
amples of atomic frequency standards—the hydrogen maser, the 
cesium beam, and the rubidium gas cell — are then discussed in terms 
of these principles and conditions. The functions of the fundamental 
parts of each device become apparent through this development. 


Robinson, D. C., Determining of load errors in universal test- 
ing machines due to speed of testing, Mater. Res. Stand. 11, 
No. 1, 14-17 (January 1971). 


Key words; Dynamic loading; load error; testing machine. 


A modification of an American Society for Testing and Materials 
(ASTM) proposed method of determining errors in testing machines 
was used to measure the indicator errors caused by dynamic loading 
for three universal machines of different design. For each of the 
machines, load errors due to the lag of the indicator were found to 
be linearly related to the rate of applied load. Based on the experi- 
mental results, it appears that for testing machines of various design, 
the response of the indicating system to uniformly increasing loads 
resembles that of a proportional-error servomechanism which has a 
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“velocity-lag” steady-state error when following a ramp function 
inpui. A discussion is given of the variables which are found to be 
significant in determining the load errors for the three machines 
investigated. The essential components in the auxiliary load- 
measuring system were a load cell, a fast response digital voltmeter 
and a high speed still camera. The alternate method used was found 
to provide higher quality photographic information from which 
to determine the indicator errors than that obtained when using a 
moving picture camera as proposed in the ASTM procedure. 


Rowland, G. A., The increasing role of states in building codes 
and standards, Proc. 26th. Annual Conf. Reinforced Plastics, 
Composites Division, Washington, D.C., February 9-12, 1971, 
Section 5—B, pp. 1—3 (February 10, 1971). 


Key words: Interstate; intra-state; national conference of states 
on bldg. codes & stds.; uniformity. 


Building codes promulgated, modified, and interpreted on a local 
basis have been deemed to be major deterrents to the develop- 
ment of reasonably large aggregated markets for many new or 
little-used materials needed in meeting the country’s present 
and future building and housing requirements. 

Since its inception in 1967, the National Conference of States 
on Building Codes and Standards has been increasingly active 
in promoting greater involvement of the States in the complex 
building regulatory system. The development of more competent 
technical procedures for accrediting new materials, including 
reinforced plastics, for acceptance statewide, and possibly nation- 
wide, rather than city-by-city will be a very important contribu- 
tion toward progress. 

The responsibilities for health and safety constitutionally rests in 
the several states and the differences in their administration of 
this responsibility is now being seriously reviewed. The Confer- 
ence in its three years of existence has been able to bring together 
delegates from the several states, establish working arrangements, 
and proceed in developing a program of study and recommenda- 
tions for consideration by the Conference in solving many of the 
present building regulatory problems. 

This national effort is a new, refreshing, coordinated approach 
towards what has been a rather chaotic, restrictive, individual 
approach to building regulations. 

A current report of the National Conference of States on Building 
Codes and Standards and the aforementioned organizational 
structure and activities is presented. 


Ruegg, F. C., A cireuit which expands the count storage 
capacity of a multichannel analyzer, Nuci. Instr. Methods 
92, 7-11, (1971). 


Key words: Multichannel scaler; nuclear instrumentation: over- 
flow counter. 


A circuit has been designed which allows the experimenter to store 
up to 10" counts/channel in a 10° count capacity multichannel 
analyzer. The circuit uses the analyzer’s second subgroup to store 
the overflow counts of the first subgroup. The circuit was designed 
to be used with an analyzer operating in the multichannel scaling 
mode. 


Scharf, K., Spectrophotometric measurement of ferric ion 
concentration in the ferrous sulphate (Fricke) dosemeter, 
Phys. Med. Biol. 16, No. 1, 77-86 (1971). 


Key words: Chemical dosimeter; chemical dosimetry; ferrous sul- 
fate; ferrous sulfate dosimeter; ferric ions; Fricke dosimeter; radia- 
tion dosimeter; radiation dosimetry; spectrophotometry; spectro- 
photometric measurements. 


A systematic error in the spectrophotometric measurement of ferric 
ion concentrations in the ferrous sulphate dosimeter may be made 
by an incorrect evaluation of a non-linear spectrophometric calibra- 
tion curve. Methods are discussed for determining the radiation- 
produced change in molarity from the actual calibration curve, and 
a method of normalization of measured absorbances is suggested. 
Normalization factors, converting measured absorbances into nor- 
malized values, can either be calculated by choosing a reference 





value of the molar extinction coefficient, or can be determined by 
comparative absorbance measurements on two spectrophotometers, 
one of them to be a precision instrument. Normalized absorbances 
are proportional to molarity and may be considered to be free of er- 
rors due to instrumental parameters and inaccuracies in acidity and 
temperature of solutions, and of errors in molarity if derived by com. 
parative measurements. 


Schramm, R. E., Corrections and calculations on a x-ray dif- 
fraction line profile: A computer program, Nat. Bur. Stand. 


(U.S.), Tech. Note 600, 34 pages (June 1971) 40 cents, SD Catalog No. 
C13.46:600. 


Key words: Computer programs, Fourier analysis; nickel steels; 
x-ray analysis; x-ray diffraction; statistics. 


This computer program was written to perform corrections and make 
calculations on an x-ray diffraction profile before Fourier analysis. 
The corrections are for background and for variations of intensity 
with the Bragg angle. Also calculated are the separation of the 
Ka; — Ka, doublet and the centroid and position of peak maximum 
with their standard deviations. There is also an option to smooth the 
profile. 


Scott, W. W., Jr., New coaxial RF-DC ammeter, JEEE Trans. 
Instr. Meas. IM-19, No. 4, 318-323 (Nov. 1970). 


Key words: Ammeter; rf-de, coax; thermopile radiation; infrared. 


A unique new type of rf ammeter is described. The ammeter com- 
bines broad frequency coverage with broad current range and is 
designed for use in a 50-ohm coaxial transmission line. Thermal 
compensation is employed to reduce effects of most self-heating 
errors. 

The ammeter consists of a cylindrical ellipsoidal silver reflector of 
infrared energy with a cylindrical thick-film heater along one focus 
and a heat sensing thermopile along the other focus. Advantages 
over older ammeter designs include a shielded construction with 
inherent electrical isolation between the rf and de circuits, a large 
de output (approximately 1 to 100 millivolts) corresponding to a 
decade dynamic range (0.5 to 5 amperes), and a frequency range 
from dc to about 1 GHz. 

This coaxial ammeter is believed to be the first to have very low 
SWR (1.02) in a 50-ohm transmission line. Ammeters of lower current 
range can be built, but at a penalty of increased SWR due to the 
larger heater resistance required. 


Selby, M. C., The Bolovac and its applications, JEEE Trans. 
Instr. Meas. IM-19, No. 4, 324-331 (November 1970). 


Key words: Bolovac; microwave current measurements; microwave 
power measurements; microwave voltage measurements. 


The Bolovac is a novel device for the standardization and measure- 
ment of TEM voltage, current, and power from below 1 MHz through 
18 GHz (potentially through 36 GHz). It can be used for attentuation 
measurement and as a superior impedance termination at all the 
above-mentioned frequencies. To date it is the first and only practical 
means of calibrating microwave current indicators and of evaluating 
oscilloscopes displaying 25-ps or shorter rise-time pulses. 

It is felt that application of the Bolovac wiil radically improve 
power measurements. It reduces the problems arising from mis- 
match-error uncertainties, one of the major present sources of 
concern in power measurement; measurement of complex reflection 
coefficients; complex computations; use of error-limit charts; and 
measurement of effective efficiency and calibration factors. Its 
application should radically increase optimum accuracy and meas- 
urement-time efficiency in calibration facilities outside the National 
Bureau of Standards. 

The Bolovac together with an Rf generator can be used as a source 
of known voltage. It can also be used as an absorption power meter. 
It needs no RF calibration. It has a range as high as 40 dB depending 
on the desired precision and on whether one or more sensors 
(special bolometric disks) are used. Voltages down to 0.1 volt with 
accuracies of 0.5 percent and better (and corresponding power in 
50-ohm systems) have been measured. 

The paper describes the general features and application of the 
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Bolovac and lists references to more comprehensive papers on the 
subject. 


Sieck, L. W., Searles, S., Ausloos, P., High-pressure photoioni- 
zation mass spectrometry. Photoionization of propane at 
11.6-11.8 eV. Formation and reactivity of the (C;Hs); dimer 
ion, J. Chem. Phys. 54, No. 91-95 (Jan. 1, 1971). 


Key words: lon-molecule reaction; kinetics; mass spectrometry; 
photoionization; propane; radiation chemistry. 


The major reaction path of the propane molecular ion with propane 
was found to be the formation of the dimer ion (C;Hg)} via a termolec- 
ular mechanism; 


C3H 
C3Hy + CsHs ——> (C3H,)f + CsHs 


In addition, C;H{ and C,H} were also found as minor reaction 
products at lower pressures. The reactions of the dimeric ions with 
ethylene and NO were also investigated. The charge exchange 
reaction 


(CsHs)f +NO—NO* + 2CsHs 


was found in propane-NO mixtures, suggesting a recombination 
energy in excess of 9.24 eV. The formation of C;HsNO* was also 
detected at higher total pressures. The dimeric ion was also found 
to transfer H, to ethylene without affecting the structural integrity 
of the Carbon skeleton, 


(C3Hs)3 3 CoH, -—> CeH;, - C2H¢ 


indicating that this species exhibits the chemical behavior of a 
saturated hydrocarbon ion. 


Sixsmith, H., Giarratano, P., A miniature centrifugal pump, 
Rev. Sci. Instr. 41, No. 11, 1570-1573 (Nov. 1970). 


Key words; Centrifugal; helium; nitrogen; pump; pump design; 
pump performance. 


A miniature centrifugal pump designed to circulate helium in a flow 
loop is described. It is designed to operate in the liquid helium 
temperature range under ambient pressures up to 100 atmospheres 
(107N/m?), producing a maximum flow rate of 50 cm*/s and a maxi- 
mum pressure rise of 0.5 psi (3400 N/m?). The operating character- 
istics obtained from preliminary tests with cold nitrogen gas and an 
analysis of the performance are included. 


Slattery, W. J., Editor, An index of U.S. voluntary engineering 
standards, Nat. Bur. Stand. (U.S.), Spec. Publ. 329, 1,000 pages 
(Mar. 1971) $9.00, SD Catalog No. C13.10:329. 


Key words: Engineering standards, index of; index of standards, 
recommended practices, specifications, and test methods; Key- 
Word-In-Context index of voluntary standards; standards, volun- 
tary, index of. 


This computer-produced Index contains the permuted titles of 
more than 19,000 voluntary engineering and related standards, 
specifications, test methods, and recommended practices, in 
effect as of December 31, 1969, published by some 360 U.S. tech- 
nical societies, professional organizations, and trade associations. 
The title of each standard can be found under all the significant 
key words which it contains. These key words are arranged alpha- 
betically down the center of each page together with their surround- 
ing context. The date of publication or last revision, the standard 
number, and an acronym designating ihe standards-issuing organi- 
zation appear as part of each entry. A list of these acronyms and 
the names and addresses of the organizations which they represent 
are found at the beginning of the Index. 


Smith, R. W., Jr., Conference Chairman, Editor, Precoordination- 
basis for industrialized building, (Proc. Conf. Gaithersburg, 
Md., Sept. 24-26, 1969), Nat. Bur. Stand. (U.S.), Bldg. Sci. Ser. 
32, 136 pages Van. 1971) $1.50, SD Catalog No. C13.29/2:32. 


Key words: Building; components; precoordination; standards. 


The Conference entitled ‘“Precoordination—Basis for Industrial- 
ized Building’ was held at the National Bureau of Standards, 
Gaithersburg, Md., on September 24-26, 1969. The Conference 
was sponsored by the American National Standards Institute’s 
Committee A62, Precoordination of Building Components and 
Systems, to explore the standards required to establish a basis 
for an industrywide system of building using interchangeable 
components. Coordinated components, conforming to these stand- 
ards, will be compatible and interchangeable in both dimension 
and function and thereby offer unlimited opportunities for product 
and material selection as well as design flexibility. 


Somes, N. F., Dikkers, D. R., Boone, T. H., Lubbock tornado: 
A survey of building damage in an urban area, Nat. Bur. 
Stand. (U.S.), Tech. Note 558, 38 pages (Mar. 1971) 50 cents, SD 
Catalog No. C13.46:558. 


Key words: Anchorage; building performance; glazing; hailstones; 
masonry; mobile homes; roofs; structural engineering; wind load. 


The Building Research Division of the National Bureau of Standards’ 
Institute of Applied Technology sent a three-man team to investi- 
gate the damage to buildings and other structures caused by the 
tornado which struck Lubbock, Texas, on May 11, 1970. The team 
members—the authors of this report—carried out photographic 
surveys on the ground and from a helicopter on the days of May 14, 
15, and 16, 1970. The report is based largely on data gathered 
during this period but includes some data provided by other agen- 
cies and individuals whose assistance is acknowledged in the 
report. The report concludes that current good practice in the 
design and construction of buildings and mobile homes would 
have greatly reduced the damage observed at Lubbock. It also 
notes that natural disasters provide full-scale tests of buildings 
and urges the development of performance criteria with respect 
to wind loads for certain building elements. 


Spiegel, V., Jr., Murphey, W. M., Calculation of thermal neutron 
absorption in cylindrical and spherical neutron sources, 
Metrologia, 7, No. 1, 34-38 (Jan. 1971). 


Key words: Manganous sulfate bath; neutron source absorption; 
neutron source calibration; thermal neutron absorption. 


A calculation of the thermal neutron self-absorption for cylindrical 
or spherical neutron sources has been made. The calculations are 
confirmed by the experimentally-measured difference in manganous 
sulfate bath activity for bare and cadmium-covered Pu-Be and 
Am-Be neutron sources. The calculation is done in single interaction 
approximation and assumes that the incident thermal neutron 
flux is isotropic. The source material may be fissionable and be 
covered by up to three cladding materials. A computer program 
has been written for the numerical calculations. 


Speigel, V., Jr.,. Murphey, W. M., Computer code for the calcu- 
lation of thermal neutron absorption in spherical and cy- 
lindrical neutron sources, Nat. Bur. Stand. (U.S.), Tech. Note 
576, 25 pages (May 1971) 35 cents, SD Catalog No. C13.46:576. 


Key words: Manganous sulfate bath calibration of neutron sources; 
neutron; neutron standards. 


A computer code has been written in FORTRAN IV for the calcu- 
lation of thermal neutron absorption in spherical and cylindrical 
neutron sources. The formalism of the calculation, the structure of 
the computer code, a listing of the code, and some sample results 
are presented. The comparison of the results of this calculation 
to experiment appears elsewhere (1). 


Spijkerman, J. J., Travis, J. C., Pella, P. A., DeVoe, J. R., Pre- 
liminary study on the characteristics and design parameters 
for a Mossbauer resonant detector, Nat. Bur. Stand. (U.S.), 
Tech. Note 541, 65 pages (Jan. 1971), 65 cents, SD Catalog No. 
C13.46:541. 


Key words: Conversion electrons; iron; Méssbauer Spectroscopy; 
resonant detector. 
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Progress in the design and fabrication of a resonant detector for 
Missbauer Spectroscopy is described. This report begins with a 
review of all of the methods of detection for this spectroscopy and 
describes the expected advantages of the resonant detector. If 
one uses conversion electron detection, considerable enhancement 
in signal to noise ratio and decrease in linewidth may be realized. 
Efforts to produce an iron bearing resonant material are described. 


Sullivan, D. B., Peterson, R. L. Kose, V. E., Zimmerman, J. E., 
Generation of harmonics and subharmonics of the 
Josephson oscillation, J. Appl. Phys. 41, No. 12, 4865-4873 
(November 1970). 


Key words: Josephson effect; superconductivity. 


The observation of harmonics and subharmonics of the Josephson 
oscillation is shown to be in agreement with a rather simple model 
of the junction. The generation of harmonics provides an explana- 
tion of induced steps in the current-voltage characteristic which 
occur at submultiples of the usual induced step voltages. The sub- 
harmonic oscillation is seen to be a relaxation-like process which 
can be easily understood in terms of a mechanical analog. 


Taylor, B. N., Repert on the International Conference on 
Precision Measurement and Fundamental Constants, held 
at the National Bureau of Standards, Gaithersburg, Md., Aug. 3-7, 
1970, Metrologia 7, No. 1, 39-42 (Jan. 1971). 


Key words: Fundamental constants; least squares adjustments; 
precision measurements. 


This is a report on the International Conference on Precision 
Measurement and Fundamental Constants held at the National 
Bureau of Standards, Gaithersburg, Maryland, August 3-7, 1970. 
It includes a brief discussion of the purpose and significance of 
the Conference as well as some of the organizational details as 
they pertain to the Conference sponsors, topics covered, publica- 
tion of the Conference Proceedings, etc. The bulk of the report 
is devoted to a summary by topic of the Conference highlights 
and conclusions. 


Thurber, W. R., Determination of deep impurities in silicon 
and germanium by infrared photoconductivity, Nat. Bur. 
Stand. (U.S.), Tech. Note 570, 13 pages (Mar. 1971) 25 cents, SD 
Catalog No. C13.46:570. 


Key words: Deep impurities; germanium; infrared; photoconduc- 
tivity; photoresponse; semiconductors; silicon. 


The feasibility of using infrared photoresponse and photoconduc- 
tivity measurements to study deep impurities in germanium. and 
silicon is examined by reviewing the literature. It is concluded 
that photoconductivity is useful in detecting the presence of specific 
impurities because each impurity has a long wavelength cut off 
in response associated with its ionization energy. However, when 
there are several deep impurities in the same specimen, it is difh- 
cult to be certain of detecting each one because some have broad 
cut offs and many have nearly the same ionization energies. Photo- 
conductivity as a general technique has serious limitations for 
determining the total concentration of deep impurities. The equa- 
tions for determining impurity concentration from the magnitude 
of the photoconductivity signal depend on the relative influence 
of deep and shallow centers. Equations are derived for several 
situations and experimental results from the literature are dis- 
cussed for each one. Only uncompensated centers are available 
for photoionization and therefore the total concentration can not 
be obtained directly. In some situations the response due to a 
deep center is independent of its concentration. Other techniques 
for studying deep impurities are discussed briefly. 


Tighe, N. J., Microstructure of fine-grain ceramics, (Proc. 
15th Sagamore Army Materials Research Conf., Sagamore Confer- 
ence Center, Raquette Lake, N.Y., Aug. 20-23, 1968), Chapter in 
Ultrafine-Grain Ceramics, pp. 109-133 (Syracuse University Press, 
Syracuse, N.Y., 1970). 


Key words: Alumina; AlO;; ceramics; electron microscopy; fine- 
grain ceramics; ion bombardment; magnesia; MgO; microstruc- 
ture; rock; zirconia. 





This describes the use of transmission electron-microscopy to 
characterize the microstructure of fine-grain ceramics. Observa- 
tions have been made on a number of polycrystalline materials 
including alumina, magnesia, zirconia, metal-ceramic composites, 
and rock specimens. 

Thin sections were prepared by ion bombardment. In these sec- 
tions grain boundaries, pores, impurity precipitates and dislo- 
cations could be observed directly. Crystalline second-phase 
material formed as grains and small precipitates could be identi- 
fied by means of electron diffraction. The method of specimen 
preparation and the results obtained from the observation of the 
specimens will be discussed. 


Varner, R. N., Peavy, S. T., Test problems and results for 
OMNITAB II, Nat. Bur. Stand. (U.S.), Tech. Note 551, 190 pages 
(Dec. 1970) $1.50, SD Catalog No. C13.46:551. 


Key words: Accuracy; ANSI] FORTRAN; computer system imple- 
mentation; examples; OMNITAB II; sottware; test problems. 


The lack of test problems and results for many software packages 
is a great hindrance to both the systems programmer and the general 
user. In this publication a set of fifty-two test problems and results for 
the OMNITAB II system is provided to assist individuals in checking 
the implementation of the OMNITAB II program on their particular 
computer. The general user will also find these descriptive examples 
instructive in the use of OMNITAB commands. 


Wachtman, J. B., Jr.. Standard materials for measurements on 
ceramics, Am. Ceram. Soc. Bull. 50, No. 3, 242-247 (March 1971). 


Key words: Ceramics; chemical composition; mechanical properties; 
pressure; properties; standard reference materials; standards; 
temperature; thermal. 


Standard materials having well determined values of specific proper- 
ties are useful in establishing the validity of techniques for character- 
ization and property measurements on ceramics. For some properties 
certified materials are available “off-the-shelf”; for certain other 
properties certification is in progress and a few samples may be 
available on a collaborative basis; while for still other properties 
serious problems remain. Some recent progress and current work 
on standard materials is briefly reviewed for the areas of tempera- 
ture, pressure, radiation, chemical composition, thermal emittance, 
thermal expansion, thermal conductivity, vapor pressure, elasticity 
and viscosity. 


Wagman, D. D., Evans, W. H., Parker, V. B., Halow, I., Bailey, 
S. M., Schumm, R. H., Churney, K. L., Selected values of chem- 
ical thermodynamic properties—Tables for elements 54 
through 61 in the Standard order of arrangement, Nat. 
Bur. Stand. (U.S.), Tech. Note 270-5, 49 pages (Mar. 1971) 55 cents, 
SD Catalog No. C13.46:270-5. 


Key words: Enthalpy; entropy; Gibbs energy of formation; hafnium 
compounds; heat of formation; niobium compounds; scandium 
compounds; tantalum compounds; titanium compounds; vanadium 
compounds; yttrium compounds; zirconium compounds. 


Contains tables of values for the standard heats and Gibbs (free) 
energies of formation, entropies and enthalpies at 298.15 K and 
heats of formation at 0 K for compounds of vanadium, niobium, 
tantalum, titanium, zirconium, hafnium, scandium, and yttrium 
(elements 54-61 in the Standard Order of Arrangement). These 


tables are a continuation of the comprehensive revision of NBS 
Circular 500. 


Wait, D. F., The precision measurement of noise temperature 
of mismatched noise generators, /EEE Trans. Microwave 
Theory Tech. MTT-18, No. 10, 715-724 (Oct. 1970). 


Key words: Compensation method; cryogenic noise generators; 
mismatch error; noise temperature measurement. 


By refining the technique and analysis of an earlier paper, it is 
shown that a compensation method is the basis of the most pre- 
cise measurement of noise temperature, especially for cryogenic 
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noise generators. The measurement system is adjusted so that it 
is in thermal! equilibrium with the generator under test. As is typical 
for the compensation method, the generator under test may be mis- 
matched and the comparison errors and the detector limitations 
can be verified experimentally. For generators with small reflec- 
tion coefficients (| |< 0.1), the comparison error is very close 


to the theoretical minimum, namely twice the system resolution. 
For reflection coefficients as large as 0.5, the comparison error 
increases in atypical system to about four times the system 
resolution. 


Wallace, B., Brown, W. E., Stoichiometric composition of 
whitlockite, J. Dental Res. 50, No. 2, 343-346 (March-April 1971). 


Key words: Ca/P ratio; pyrolysis; pyrophosphate; stoichiometry; 
tricalcium phosphate; whitlockite. 


The mineral whitlockite has been described as having Ca/P ratios 
of 3/2 and 10/7. This study used the pyrolytic formation of pyro- 
phosphate and whitlockite from calcium orthophosphates and 
their mixtures with CaCQOs, to determine this Ca/P ratio. Its mean 
value of 1.506+.015 indicated a formula of Cas(PO,)s. 


Walther, H., Hall, J. L., Tunable dye laser with narrow spectral 
output, Appl. Phys. Letters 17, No. 6, 239-242 (September 15, 
1970). 


Key words: Birefringent filter; dye laser; frequency control. 


A tunable dye laser with narrow-banded spectral output in the order 
of 0.01 A or even smaller is described. This narrow spectral output 
is obtained by means of a birefringent filter (Lyot filter) which was 
inserted into the laser cavity. The properties of this laser setup have 
been investigated and are described. 


Wampier, R. H., A report on the accuracy of some widely 
used least squares computer programs, J. Am. Stat. Assoc. 
65, No. 330, 549-565 (June 1970). 


Key words: Computer programs; curve fitting; Gram-Schmidt orthog- 
onalization; Householder transformations; iterative refinement; least 
squares; linear equations; orthogonalization; orthogonal polynomi- 
als; regression; rounding error; stepwise regression. 


Two linear least squares test problems based on fifth degree poly- 
nomials have been run on more than twenty different computer 
programs in order to assess their numerical accuracy. The programs 
tested, all in present-day use, included representatives from several 
statistical packages as well as some from the SHARE library. 
Essentially four different algorithms were used in the various pro- 
grams to obtain the coefficients of the least squares fits. The tests 
were run on several different computers, in double precision as well 
as single precision. By comparing the coefficients reported, it was 
found that those programs using orthogonal Householder transforma- 
tions or Gram-Schmidt orthonormalization were much more accurate 
than those using elimination algorithms. Programs using orthogonal 
polynomials (suitable only for polynomial fits) also proved to be 
superior to those using elimination algorithms. The most successful 
programs accumulated inner products in double precision and made 
use of iterative refinement procedures. In a number of programs, the 
coefficients reported in one test problem were sometimes completely 
erroneous, containing not even one correct significant digit. 


Waxler, R. M., Laser glass composition and the possibility of 
eliminating electrostrictive effects, /EEE, J. Quantum Electron. 
QE-7, No. 4, 166-167 (April 1971). 


Key words: Electrostriction; glass; lasers; photoelasticity. 


Many authors have stated that electrostrictive self-trapping of the 
light beam initiates laser damage in solid transparent dielectrics. 
By systematically varying composition, it may be possible to find 
a glass in which electrostrictive effects are eliminated. 


Weiss, A. W., A review of theoretical developments in atomic 
f-values, (Proc. 2nd. Intern. Conf. Beam-Foil Spectroscopy in 
Nuclear Instruments and Methods, Lysekil, Sweden, June 8-12, 
1970), Part Il, Lifetimes and transition probabilities, 25-119 in 
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Nucl. Instr. Methods 90, 121—131 (North-Holland Publ. Co., Amster- 
dam, The Netherlands, 1970). 


Key words: Atomic lifetimes; atomic spectra; configuration interac- 
tion; oscillator strengths. 


Theoretical work on atomic f-values are reviewed with particular 
emphasis on the relation to recent beam-foil measurements. One- 
electron models currently in use are reviewed and their limitations 
and range of applicability are discussed. Many-electron, or multi- 
configuration, models are also described and their results compared 
with experiment. Regularities and irregularities along isoelectronic 
sequences will be discussed, and some current problems will be 
presented. 


Weiss, B. Z., Meyerson, M. R., Fatigue crack initiation and 
propagation in chromium diffusion coated Ti bearing steel, 
J. Iron Steel Inst. 208, No. 12, 1069-1077 (December 1970). 


Key words: Chromium diffusion coating; columnar grains; fatigue 
crack initiation; fatigue crack propagation; residual stress; titanium 
bearing steels. 


Chromium diffusion coatings on Ti-bearing steel frequently lead to 
columnar growth beneath the Cr-rich layer. The grains are prefer- 
entially oriented. The depth of the columnar zone does not depend 
on the chromizing time. A high rate of cooling prevents columnar 
growth. No grain boundary diffusion or carbide formation at the 
grain boundaries was observed. The Cr-rich later is structurally 
homogenous and consists of a-solid solution. The hardness of 
the coating does not vary with the chromizing time. 

Compressive residual stresses were found in the Cr-rich later. 
The stresses, whose magnitude does not change with the depth 
of the layer, may be decreased by additional heat treatment. 
Three modes of fatigue crack nucleation were observed all be- 
neath the Cr-rich layer. Cracks were formed after 10% to 18% 
of total fatigue lifetime. Crack propagation in its initial stages 
is primarily dependent on the mode of nucleation. In the later 
stage propagation is dependent on stress and grain size. Fatigue 
properties may be improved considerably by a factor of 3 to 5 by 
additional “normalizing” after chromizing. 


West, E. D., Data analysis for isoperibol laser calorimetry, 
Nat. Bur. Stand. (U.S.), Tech. Note 396, 34 pages (Feb. 1971) 40 
cents, SD Catalog No. C13.46:396. 


Key words: Calorimetry; laser; laser calorimetry; laser energy; 
laser power. 


Isoperibol calorimeters (those operating in a constant-temperature 
environment) are used to measure the power and energy in laser 
beams relative to electrical standards. The derivation of the basic 
formula is reviewed. Two methods are presented for analyzing the 
data taken at equal time intervals: (1) An approximate manual 
method with criteria for avoiding significant errors of approxima- 
tion and (2) A least squares method for use with automatic digital 
computers. 


White, H. J., Jr., Federal information processing standards 
index, January 1, 1971, Nat. Bur. Stand. (U.S.), Fed. Info. 
Process. Stand. Publ. (FIPS Pub) 12, 143 pages, $1.50, SD Catalog 
No. C13.52:12. 


Key words: American National Standards; computers, data ele- 
ments and codes; data processing systems; Federal Information 
Processing Standards; management information systems; Inter- 
national Organization for Standardization; standards; U.S. Gov- 
ernment. 


This publication provides material concerning standardization 
activities in the area of information processing at the Federal, 
National and International levels. Also included are related policy 
and procedural guideline documents. A list of Federal Govern- 
ment participants involved in the development of Federal Infor- 
mation Processing Standards is provided. 


Wiederhorn, S. M., Johnson, H., Effect of pressure on the 
fracture of glass, J. Appl. Phys. 42, No. 2, 681—684 (Feb. 1971). 





Key words: Deep submergence; tracture; fracture energy; glass; 
high pressure; strength. 


The fracture surface energies of three glass compositions were 
measured as a function of ambient pressure and were found to be 
independent of pressure, to 20 kbar, suggesting no change in the 
fracture mechanism. The mechanical behavior of glass thus differs 
from that of plastics or metals which are observed to become stronger 
and more ductile with increasing pressure. The difference in fracture 
behavior is believed due to the fact that fracture of glass is essentially 
a volume conserving process in contrast to metals and plastics for 
which volume expansion occurs. 


Woelfel, J. B., Paffenbarger, G. C., Expanding and shrinking 


77-year-old dentures: report of cases, J. Am. Dental Assoc. 81, 
No. 6, 1342-1348 (December 1970). 


Key words: Acrylic; adaption; contraction; dentures; epoxy; ex- 
pansion; hard rubber; retention; serviceability; stability. 


A set of complete dentures made of hard rubber shrank and warped 
during seven years of observation. The greatest shrinkage occurred 
on the lower molar-to-molar span 1.47% (0.75 mm) and the least 
0.13% (0.09 mm) on the upper flange-to-flange span. A similar set of 
complete dentures made of an epoxy resin expanded but warped 
very little during the seven years. The greatest expansion 1.12% 
(0.56 mm) occurred on the lower molar-to-molar span and the least 
0.86% (0.60 mm) on the upper flange-to-flange span. There was 
little or no difference in the serviceability of these dentures as judged 
by changes in centric and eccentric occlusion, stability, retention, 
and condition of the tissues and by the dentist’s and patient’s 
appraisals. Both sets of dentures deteriorated about the same in 
service. The apposition of the supporting tissues and the dentures 
were about the same as determined by pressure indicator paste 
patterns. The conclusion is, therefore, that the supporting tissues 
followed the dimensional changes in the dentures to a great extent. 


Yakowitz, H., Fiori, C. E., Michaelis, R. E., Standard reference 
materials: Homogeneity characterization of Fe-3Si alloy, 
Nat. Bur. Stand. (U.S.), Spec. Publ. 260-22, 30 pages (Feb. 1971) 
35 cents, SD Catalog No. C13.10:260-22. 


Key words: Fe-3Si alloy; homogeneity testing; metallography; 
microprobe analysis; spectrometric analysis; standard reference 
materials. 


An alloy of iron-3.22 wt. pet. silicon (Fe-3Si) was characterized 
with regard to chemical homogeneity of iron and silicon at the 
micrometer level of spatial resolution. This alloy is satisfactory for 
use as a homogeneous standard for electron probe microanalysis. 
The samples were cut from coarse-grained sheet stock to a final 
size of about 3 mm x 3 mm x 0.28 mm thick. Homogeneity was 
checked by means of quantitative raster scanning in which a square 
matrix (1.1 mmx 1.1 mm) of individual points is analyzed by the 
microprobe. Each matrix represents 400 separate analyses. Usually, 
the same matrix was rerun so that each point was sampled twice. 
The coefficient of variation for both the iron and silicon is less 
than one percent. Quantitative microprobe analysis was also car- 
ried out on this alloy giving a silicon content of 3.14% and an iron 
content of 96.9%. 


Yokel, F. Y., Dikkers, R. D., Strength of load bearing masonry 
walls, J. Struct. Div. Am. Soc. Civil Engineers 97, No. ST5, 1593- 
1609 (May 1971). 


Key words: Bricks: buckling; deflection; loads; forces; masonry; 
moments; slenderness ratio; stability; structural engineering; walls. 


A theory is presented, for the rational analysis of masonry walls 
under eccentric vertical compressive load, or under a combination 
of vertical compressive and transverse load. The theory is based 
on established concepts, which have been modified to account for 
the properties of masonry. Similar methods of rational analysis 
have been adopted for the design of steel columns and are in the 
process of being ad»pted for reinforced concrete columns. An 
investigation of the results of 192 tests of full-scale walls of various 
types of masonry construction established the feasibility of apply- 





ing this theory to closely predict the strength of these walls. The 
analysis accounts for variables not presently considered in the 
design of masonry structures. The correlation between test results 
and theory is illustrated by several typical examples. 


Yokel, F. Y., Mathey, R. G., Dikkers, R. D., Compressive strength 
of slender concrete masonry walls, Nat. Bur. Stand. (U.S,), 
Bldg. Sci. Ser. 33, 32 pages (Dec. 1970) 40 cents, SD Catalog No. 
C13.29/2:33. 


Key words: Buckling; compressive strength; concrete block walls: 
elastic stability; flexural strength; masonry walls; reinforced con- 
crete masonry walls; slenderness effect; structural stability. 


Sixty reinforced and unreinforced concrete masonry walls of 
different slenderness ratios were tested to failure under vertical 
loads applied axially and at various eccentricities. Prism specimens, 
made of similar masonry units and mortars, were also tested under 
the same loading conditions. Analysis of test results indicates that 
wall strength can be conservatively predicted by evaluating cross- 
sectional wall capacity on the basis of prism strength and reducing 
the capacity for slenderness effects by evaluating the added moments 
attributable to wall deflection. Test results were also compared 
with allowable loads computed in accordance with the current 


NCMA standard. 


Yokel, F. Y., Mathey, R. G., Dikkers, R. D., Strength of masonry 
walls under compressive and transverse loads, Nat. Bur. 
Stand. (U.S.), Bldg. Sci. Ser. 34, 74 pages (Mar. 1971) 70 cents, 
SD Catalog No. C.13.29/2:34. 


Key words: Brick; cavity walls; composite walls; compressive 
strength; concrete block; flexural strength; masonry; mortar; slen- 
derness effects: standards; structural stability; walls. 


Ninety walls of 10 different types of masonry construction were 
tested under various combinations of vertical and transverse load. 
It is shown that the effect of vertical load and wall slenderness on 
transverse strength can be predicted by rational analysis. The 
analysis is based on established theory which has been extended 
to account for the properties of masonry. Similar methods of rational 
analysis have been adopted for the design of steel structures and are 
presently being considered for reinforced concrete structures. 


Zimmerman, J. E., Recent Developments in superconducting 
devices, J. Appl. Phys. 42, No. 1, 30-71 (Jan. 1971). 


Key words: Cardiography; instrumentation; 
magnetometry; superconductivity. 


Josephson effect; 


A simple model of the electrical characteristics of thin-film bridge, 
point contact, and tunnel junction Josephson devices is given, along 
with some comments on their relative performance at very low and 
at very high frequencies. A particular example is the de IV character- 
istic of a point contact at the center of a parallel-disc microwave 
cavity. Some recent developments in devices and the application 
of a point-contact loop (SQUID) device to magnetocardiography is 
described. 


Other NBS Publications 


J. Res. Nat. Bur. Stand. (U.S.), 75A (Physics & Chemistry), 
No. 2 (March-April 1971), SD Catalog No. C13.22/sec. 
A:75/2. $1.25. 


Tentative set of key values for thermodynamics—Part I. Report of 
the ICSU-CODATA Task Group on Key Values for Thermody- 
namics. October, 1970. S. Sunner. 

Expanded formulation of thermodynamic scaling in the critical 
region. M. J. Cooper. 

Table of spin-orbit energies for p-electrons in neutral atomic (core) 
np configurations. W. C. Martin. 

Heat capacities of cis-1,4-polyisoprene from 2 to 360 K. S. S. Chang 
and A. B. Bestul. 

Combined low-pressure and high-pressure measurements of density 
and bulk modulus of aviation instrument oil and 2-methylbutane 
and their mixtures. J. C. Houck and P. L. M. Heydemann. 
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Crystal structure of Ca,Na2(COs3)s (shortite). B. Dickens, A. Hyman, 
and W. E. Brown. 


J. Res. Nat. Bur. Stand. (U.S.), 75A (Physics & Chemistry), 
oe oa ar Jane 1971), SD Catalog No. C13.22/sec.A:75/3. 


The solid phase photolysis and radiolysis of ethylene at 20 and 77 
K. R. Gorden, Jr., and P. Ausloos. 

Kinetic mass spectrometric investigation of the ion-molecule reac- 
tions occurring in Cy and C; alkanes following photoionization 
at 106.7 and 104.8 nm. L. W. Sieck, S. K. Searles, and P. Ausloos. 

An alumina standard reference material for resonance frequency 
and dynamic elastic moduli measurement. I. For use at 25 °C. 
R. W. Dickson and J. B. Wachtman, Jr. 

Optical and mechanical properties of some neodymium-doped laser 
glasses. R. M. Waxler, G. W. Cleek, I. H. Malitson, M. J. Doge, 
and T. A. Hahn. 

Effect of environment of viscous flow in inorganic oxide glasses. 
J. H. Simmons and P. B. Macedo. 

Phase relations in the SrO-IrO,-Ir system in air. C. L. McDaniel 
and S. J. Schneider. 

The crystal structure of BaCa(CO;), (barytocalcite). B: Dickens and 
J. S. Bowen. 

Ionization of hydrofluoric acid at 25 °C. P. R. Patel, E. C. Moreno, 
and J. M. Patel. 

Vapor pressure equation for water in the range 0 to 100 °C. A. 
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